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THESIS  ABSTRACT 

To  provide  a  new  perspective  on  the  fate  of  both  natural  organic  matter  and 
hydrophobic  organic  contaminants  (HOCs)  in  marine  sediments,  we  have  investigated  the 
relationship  between  radiocarbon  (14C)  age  and  the  different  modes  of  association  in 
aquatic  sediments  and  soils.  Radiocarbon  is  a  sensitive  tracer  of  OM  provenance,  with 
variations  in  its  natural  abundance  reflecting  the  age  and  origin  of  material.  The  main 
objective  has  been  to  determine  the  significance  of  these  associations,  and  to  assess  how 
they  affect  the  transport,  bioavailability,  preservation  and  residence  times  of  organic 
compounds  in  the  environment. 

Our  results  indicate  that  the  majority  of  HOCs  that  persist  in  marine  sediments  are 
solvent-extractable  and  incorporation  into  insoluble  sediment  residues  is  not 
quantitatively  significant.  For  pristine  sediments,  systematic  variations  in  14C  content  are 
observed  between  different  chemically  defined  sedimentary  organic  fractions.  These 
variations  are  dependent  on  organic  matter  inputs  and/or  the  affects  of  diagenesis.  Our 
observations  also  provide  evidence  for  the  protection  of  labile  marine  carbon  by  chemical 
binding.  Finally,  the  persistence  of  u-alkanes  from  biogenic  sources  compared  to  those 
derived  from  petroleum  indicates  that  protective  matrix  association  can  play  a  crucial  role 
in  determining  the  long-term  fate  of  a  compound.  Overall,  it  is  clear  that  both  natural 
organic  compounds  and  HOCs  can  undergo  very  different  fates  depending  on  their  mode 
of  introduction  to,  and  physical  disposition  in  environmental  matrices. 
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AA14C  %0  for  organic-rich  marine  sediments  (core-top  samples 
only).  Error  bars  represent  0.1  %o  error  for  813C  values  and  10%o 
error  for  A14C. 

Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  1 10 
symbols)  and  saponified  residue  (SA-RES;  open  symbols)  from 
813C  of  bulk  sediment  (TOC;  filled  black  symbol)  expressed  as 
A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA-RES 
from  the  TOC  expressed  as  AAI4C  %o  for  open  shelf  and  river 
dominated  sediments  (core-top  samples  only).  Error  bars  represent 
0.1%o  error  for  813C  values  and  10%o  error  for  A14C. 

Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  111 
symbols)  and  saponified  residue  (SA-RES;  open  symbols)  from 
813C  of  bulk  sediment  (TOC;  filled  black  symbol)  expressed  as 
A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA-RES 
from  the  TOC  expressed  as  AA14C  %o  for  near-coast  sediments. 

Error  bars  represent  0.1  %>  error  for  813C  values  and  10%o  error  for 
A14C. 

Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  113 
symbols)  and  saponified  residue  (SA-RES;  open  symbols)  from 
8,3C  of  bulk  sediment  (TOC;  filled  black  symbol)  expressed  as 
A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA-RES 
from  the  TOC  expressed  as  AA14C  %o  for  Arctic  sediments.  Error 
bars  represent  0.1  %o  error  for  813C  values  and  10%o  error  for  A14C. 
Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  1 14 
symbols)  and  saponified  residue  (SA-RES;  open  symbols)  from 
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Figure  4.18 


Chapter  5 
Figure  5.1 

Figure  5.2 
Figure  5.3 


Figure  5.4 
Figure  5.5 
Figure  5.6 


Figure  5.7 


513C  of  bulk  sediment  (TOC;  filled  black  symbol)  expressed  as 
A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA-RES 
from  the  TOC  expressed  as  AA14C  %o  for  Antarctic  sediments 
(core-top  samples  only).  Error  bars  represent  0.1  %o  error  for 
813C  values  and  10%o  error  for  A14C. 

Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  1 17 
symbols)  and  saponified  residue  (SA-RES;  open  symbols)  from 
813C  of  bulk  sediment  (TOC;  filled  black  symbol)  expressed  as 
AS13C  %o  vs.  the  differences  in  AI4C  of  the  EX-RES  and  SA-RES 
from  the  TOC  expressed  as  AA14C  %o  for  sediment  horizons  below 
the  surface  mixed  layer.  Error  bars  represent  0.1  %o  error  for  S13C 
values  and  10%o  error  for  A14C. 


Structures  of  hydrophobic  organic  contaminants  examined  in  137 
this  study. 

Location  of  sediments  examined  in  this  study.  141 

Gas  chromatograms  of  total  lipid  extracts  from  a)  New  Bedford  148 
Harbor,  b)  Hudson  River,  c)  Chattanooga  Creek  and  d)  Eagle 
Harbor.  The  dashed  line  is  the  GC  baseline.  In  a)  and  d)  PCBs 
and  PAHs  dominate.  In  b)  and  c)  contaminants  are  designated 
with  and  biogenic  hydrocarbons  n-alkanoic  acids,  n-alkanols, 
alkanes  and  sterols  are  labeled  with  “F”,  “O”,  “A”  and  “S” 
respectively.  The  internal  standard  hexatriacontane  is  designated 
with  “is”  and  an  external  contaminant  is  labeled  with  “e” 

Percent  contribution  of  biogenic  hydrocarbons  and  contaminants  149 
(specified  in  the  text)  to  the  total  lipid  extract. 

Structures  of  hydrophobic  organic  contaminants  released  by  151 
alkaline  hydrolysis  from  various  samples  examined  in  this  study. 

Gas  chromatograms  of  saponification  acidic  extracts  from  a)  New  152 
Bedford  Harbor,  b)  Hudson  River,  c)  Chattanooga  Creek  and 
d)  Eagle  Harbor.  The  dashed  line  is  the  GC  baseline.  In  all  samples, 
the  w-alkanoic  acids  are  numbered.  The  «-alkanoic  acids  in  d)  have 
different  retention  times  as  they  were  methylated  in  the 
saponification  process.  The  internal  standard  hexatriacontane  is 
designated  with  “is”. 

Gas  chromatograms  of  saponification  neutral  extracts  from  153 

a)  New  Bedford  Harbor,  b)  Hudson  River,  c)  Chattanooga  Creek 
and  d)  Eagle  Harbor.  The  dashed  line  is  the  GC  baseline.  In  all 
samples,  the  «-alkanols  acids  are  numbered.  Contaminants 
are  designated  in  c)  and  d)  with  The  internal  standard 
hexatriacontane  is  designated  with  “is”. 
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Figure  5.8 
Figure  5.9 

Figure  5.10 

Figure  5.11 
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Figure  6.1 

Figure  6.2 

Figure  6.3 

Figure  6.4 
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Percent  contribution  of  biogenic  hydrocarbons  and  contaminants  154 
(specified  in  the  text)  to  the  saponification  extract. 

Difference  between  A14C  of  unextracted  sediment  (TOC)  and  161 
solvent-extracted  sediment  residue  (EX-RES)  as  a  function  of 
percent  TOC  that  is  comprised  of  solvent-extractable  fossil-fuel- 
derived  contaminants  measured  by  gas  chromatography. 

Difference  between  the  percent  of  total  lipid  extract  (TLE)  that  162 
is  non-GC  amenable  and  the  deviation  from  the  relationship 
between  the  percent  total  organic  carbon  (%TOC)  that  is  fossil 
-fuel  derived  calculated  by  GC-FID  and  isotope  mass  balance. 
Difference  between  AI4C  of  solvent-extracted  (EX-RES)  and  165 
saponified  (SA-RES)  sediment  residues  plotted  against  the 
difference  between  813C  of  solvent-extracted  (EX-RES)  and 
saponified  (SA-RES)  sediment  residues. 


Abundance  of  individual  free  fatty  acids  from  sediment  1 84 

horizons  0-2, 12-14  and  32-34  cm. 

Abundance  of  individual  bound  fatty  acids  from  sediment  1 84 

horizons  0-2, 12-14  and  32-34  cm. 

Abundance  of  individual  n-alkanes  from  sediment  horizons  1 85 

0-2, 12-14  and  32-34  cm. 

A14C  and  813C  compositions  of  individual  free  (filled  symbols)  1 88 

and  bound  (open  symbols)  short  (triangles  and  diamonds)  and 
long  chain  (squares)  fatty  acids  from  0-2, 12-14  and  32-34  cm 
compared  to  the  A14C  and  813C  composition  of  the  solvent 
extracted  sediment  residue  (EX-RES)  from  White  et  al.,  2005a. 

Errors  for  A14C  and  813C  are  smaller  than  the  symbols 
representing  the  data  and  thus  are  not  observed. 

A14C  and  813C  compositions  of  individual  n-alkanes  from  0-2,  191 

12-14  and  32-34  cm  compared  to  the  A14C  and  813C  composition 
of  the  solvent  extracted  sediment  residue  (EX-RES)  from  White 
et  al.,  2005a.  Errors  for  A14C  and  813C  are  smaller  than  the 
symbols  representing  the  data  and  thus  are  not  observed. 
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CHAPTER  1 


GENERAL  INTRODUCTION 


The  burial  of  organic  matter  in  recent  marine  sediments  plays  a  critical  role  in  the 
global  organic  carbon  cycle  (Figure  1.1).  Quantitatively,  there  is  more  organic  carbon 
(OC)  derived  from  non-living  organic  molecules  occurring  in  recent  marine  sediments, 
soil  humus  and  dissolved  in  seawater  than  in  all  living  land  plants  and  marine  organisms 
combined  (Hedges  et  al.,  2000).  Once  deposited  in  marine  sediments,  organic  compounds 
are  subjected  to  various  processes  that  may  influence  their  long-term  fate,  including 
sorption  to  mineral  surfaces,  entrapment  or  partitioning  into  macromolecular  organic 
matter  (MOM),  reaction  with  MOM,  and  degradation  by  microorganisms.  The 
susceptibility  of  an  organic  compound  to  these  processes  is  largely  determined  by  its 
intrinsic  reactivity  as  well  as  its  mode  of  delivery  and  physical  disposition  in  the 
sediment. 

The  main  aim  of  this  thesis  was  to  improve  upon  our  understanding  of  the  nature 
and  distribution  of  organic  compounds  derived  from  both  biogenic  and  anthropogenic 
sources  in  marine  sediments.  By  examining  the  associations  of  organic  compounds  to 
MOM  we  aimed  to  make  predictions  regarding  their  long-term  fate.  In  the  case  of 
hydrophobic  organic  contaminants  (HOCs)  this  is  of  interest  as  previous  studies  have 
shown  that  certain  HOCs  or  their  metabolites  can  bind  to  MOM  in  sediments  (Beller  et 
al.,  1998;  Richnow  et  al.,  1994;  Reddy  et  al.,  2000;  Schwarzbauer  et  al.,  2003),  but  this 
process  is  difficult  to  follow  and  quantitatively  constrain.  Binding  to  sedimentary  MOM 
has  the  potential  to  influence  the  bioavalability,  toxicity  and  transport  of  HOCs,  and  thus 
understanding  these  processes  is  of  great  interest  from  an  environmental  and  human 
health  perspective.  Associations  between  natural  organic  compounds  and  MOM  have 
also  been  described  (Farrington  and  Quinn,  1971;  Lee,  1977;  Zegouagh  et  al.,  1996)  and 
are  of  interest  as  they  may  exert  strong  controls  on  the  preservation  of  labile  organic 
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Figure  1.1  The  global  organic  carbon  cycle  (ca.  1950).  Units  =  1015  gC.  Modified  after 
Hedges,  (1992). 


compounds  (Wakeham  et  al.,  1999)  as  well  as  affect  their  residence  times  in  different 
carbon  reservoirs.  Binding  to  MOM,  for  example,  may  slow  the  delivery  of  a  terrestrial 
derived  organic  compound  to  the  ocean  through  association  with  relatively  immobile 
particulate  phases  (e.g.  soils). 

To  assess  the  associations  between  MOM  and  organic  compounds,  sequential 
chemical  treatments  are  typically  performed.  Samples  are  first  solvent  extracted  to  isolate 
compounds  with  non-covalent  associations  and  these  are  defined  as  “free”  or  “unbound”; 
those  that  remain  associated  to  MOM  and  only  released  by  harsher  chemical  treatments 
are  defined  as  “nonextractable”  or  “bound”.  It  is  important  to  note  that  not  all 
nonextractable  residues  are  covalently  associated  with  MOM  and  some  may  be 
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physically  entrapped  in  MOM  networks  or  mineral  matrices  (Gevao  et  al.,  and  references 
therein;  Knicker  and  Hatcher,  1997). 

In  order  to  determine  the  significance  of  the  associations  of  both  natural  organic 
compounds  and  HOCs  with  MOM,  I  have  examined  variations  in  the  natural  abundance 
of  radiocarbon  (14C)  in  different  sediment  fractions.  Radiocarbon  is  a  sensitive  tracer  of 
organic  matter  and  variations  in  its  natural  abundance  reflect  the  age  and  origin  of 
material.  All  recently  photosynthesized  carbon  contains  significant  amounts  of  14C,  which 
is  produced  in  the  Earth’s  atmosphere  by  cosmic  ray  interactions,  and  makes  its  way  into 
the  biosphere  mainly  via  photosynthetic  fixation  of  CO2.  Additionally,  any  natural 
product  derived  from  carbon  reservoirs  in  communication  with  the  atmosphere  since 
1952  contains  an  excess  quantity  of  14C  due  to  incorporation  of  14C  produced  by  the 
above-ground  testing  of  thermonuclear  weapons.  Due  to  the  slow  equilibrium  between 
the  atmosphere  and  surface  ocean  waters,  compounds  of  oceanic  origin  are  depleted  in 
14C  relative  to  those  of  continental  origin. 

The  5730  yr  half-life  of  14C  makes  this  isotope  ideal  for  the  detection  of  fossil 
fuel-derived  HOCs  (14C-free)  within  different  fractions  of  natural  organic  matter  in 
contemporary  environmental  matrices  (modem  I4C  content).  Using  14C,  it  is  also  possible 
to  make  quantitative  assessments  through  the  principle  of  isotope  mass  balance.  For 
example,  if  we  consider  bulk  total  organic  carbon  (TOC)  in  sediments  to  be  a  mixture  (t) 
comprised  of  exogenous  (x)  and  endogenous  («)  components,  this  can  be  represented  by 

{TOC}  =  {TOC} c  +  {TOC)n  (I) 

and  the  isotopic  composition  of  the  mixture  given  by 

Ft  =  JFx  +  (l  -  f)Fn  (II) 

where  the  Fs  represent  isotopic  abundances  (e.g.,  14C/total  carbon)  and /is  the  fraction  of 
the  mixture  that  is  exogenous.  In  terms  of  radiocarbon  abundances  (A),  we  can  write 

f  =  {At-  An) /{Ax  -  An)  (III) 

where  the  As  represent,  in  terms  of  A14C,  the  abundances  of  14C  in  the  mixture  and  in  its 
endogenous  and  exogenous  components.  This  expression  implies  that / can  be  determined 
precisely  if  there  is  a  big  isotopic  contrast  between  the  endogenous  and  exogenous 
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components  (i.e.,  for  large  values  of  Ax  -  An).  This  is  optimized  in  this  approach 
whenever  the  exogenous  material  is  of  fossil  origin  (and  thus  14C-free,  Ax  =  -1000%o)  and 
the  endogenous  material  is  modem  biological  debris  (An  ~  +200%o).  Stable  carbon 
isotope  (813C)  compositions  of  organic  matter  also  provide  insights  into  the  cycling  of 
natural  organic  compounds  in  the  environment.  Understanding  the  factors  that  control  the 
813C  of  natural  organic  matter  can  place  quantitative  constraints  on  both  source  and 
carbon  flow. 

The  long-term  fate  of  HOCs  in  marine  sediments  is  explored  in  Chapter  2  which 
examines  the  abundance,  composition,  and  mobility  of  polycyclic  aromatic  hydrocarbons 
(PAHs)  in  salt  marsh  sediments  of  Wild  Harbor,  West  Falmouth.  This  site  was  impacted 
by  a  spill  of  No.  2  fuel  oil  in  1969  (Sanders  et  al.,  1980)  and  has  since  been  the  subject  of 
numerous  historical  studies  (Peacock  et  al.,  2005  and  references  therein).  Recent  analysis 
of  total  petroleum  hydrocarbons  in  the  sediments  reveals  that  the  concentrations  of  these 
HOCs  have  not  markedly  changed  over  the  past  30  years  (Reddy  et  al.,  2002  and 
references  therein).  In  this  study,  downcore  profiles  of  PAHs  are  compared  to  a 
numerical  model  to  assess  whether  sediment-water  partitioning  and  molecular  diffusion 
in  the  interstitial  medium  are  the  dominant  processes  controlling  the  vertical  distribution 
of  PAHs  that  we  observe  today. 

Chapter  3  builds  upon  this  previous  study  examining  the  fate  of  PAHs  at  the  West 
Falmouth  site  and  utilizes  variations  in  the  natural  abundance  of  radiocarbon  (A14C)  to 
detect  and  quantify  No.  2  fuel  oil  residues  in  different  fractions  of  sedimentary  organic 
matter.  This  study  examines  the  utility  of  radiocarbon  as  an  inverse  tracer  of  HOCs  in 
marine  sediments  where  the  input  of  HOCs  is  well  constrained  in  time  and  space.  The 
sensitivity  of  this  approach  is  assessed  and  isotope  mass  balance  calculations  of  the  total 
organic  carbon  (TOC)  in  the  sediments  derived  from  solvent-extractable  petroleum  are 
compared  to  values  calculated  from  a  fractional  abundance  of  total  petroleum 
hydrocarbon  concentrations  (measured  by  gas  chromatography)  as  a  function  of  TOC 
content  (determined  by  elemental  analysis). 
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Examination  of  A14C  values  and  molecular  compositions  of  different  sediment 
fractions  and  extracts  from  West  Falmouth  (Chapter  3)  not  only  provides  information 
about  petroleum  residues,  but  also  emphasizes  the  persistence  of  natural  organic 
compounds  (alkanes,  fatty  acids,  alkanols  and  sterols)  in  sediments.  In  some  instances, 
these  natural  organic  compounds  strongly  influence  the  A14C  values  of  the  bulk  residues 
and  extracts  obtained,  and  provide  information  about  the  sources  and  fate  of  these 
compounds.  Chapter  4  therefore,  is  aimed  at  investigating  the  relationships  between  14C 
content  and  the  different  modes  of  association  of  natural  organic  matter  present  in  a 
variety  of  surficial  sediments  from  near-shore  and  continental  margin  sites.  Sediments 
deposited  under  both  oxic  and  anoxic  bottom  waters,  containing  variable  contributions 
from  autochthonous  and  allochthonous  (terrestrial,  ancient  and  petroleum  carbon) 
sources,  subjected  to  different  degrees  of  diagenesis  were  examined  in  terms  of  the  nature 
and  distribution  of  organic  matter. 

With  this  improved  understanding  of  A14C  distribution  in  different  sediment 
fractions  from  uncontaminated  sites,  we  return  to  examine  sediments  impacted  by  various 
fossil  fuel-derived  HOCs  including  petroleum,  creosote,  coal  tar,  polychlorinated 
biphenyls  (PCBs),  benzotriazoles  (BZTs),  2,2-bis(chlorophenyl)-l,l,l-trichlorethane 
(DDT)  and  other  chlorinated  pesticides.  Radiocarbon  and  stable  carbon  isotopic  analyses, 
are  again  utilized,  as  in  the  study  of  the  West  Falmouth  sediments  (Chapter  3),  to 
determine  if  binding  of  these  more  reactive  HOCs  is  significant  in  marine  sediments. 
Structures  of  HOCs  and  their  metabolites  released  from  different  fractions  of  sediment 
are  also  examined  and  potential  binding  mechanisms  discussed. 

The  investigations  in  the  aforementioned  chapters  rely  on  bulk  A,4C 
measurements  of  sediments,  residues  and  extracts.  In  order  to  provide  a  molecular-level 
basis  for  our  findings,  we  return  to  the  West  Falmouth  site  and  perform  compound- 
specific  14C  analysis  of  individual  short  (w-C  14.1s)  and  long  chain  (n-C24-28)  fatty  acids  and 
long  chain  (H-C26-C32)  alkanes  from  contaminated  and  uncontaminated  sediment  horizons 
as  well  as  those  deposited  before  and  after  the  bomb-spike.  Natural  abundance  14C 
measurements  of  specific  biomarkers  provide  a  powerful  tool  to  determine  the  sources  of 
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organic  compounds  in  contemporary  environmental  matrices  (Eglinton  et  al.,  1997; 
Pearson  and  Eglinton,  2000).  Through  this  work,  we  aim  to  not  only  provide  information 
about  the  sources  of  organic  matter  to  saltmarsh  sediments,  but  to  also  conclusively 
determine  whether  any  petroleum-derived  w-alkanes  persist  at  the  site. 

In  Chapter  7,  results  presented  in  this  thesis  are  summarized  and  suggestions  are 
made  for  future  research  directions. 
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ABSTRACT 

Petroleum-derived  hydrocarbons  continue  to  persist  in  Wild  Harbor,  West 
Falmouth,  MA,  following  a  spill  of  No.  2  fuel  oil  in  1969  from  the  barge  Florida.  Recent 
analysis  of  marsh  sediments  revealed  that  residues  of  degraded  oil  are  present  with 
concentrations  of  total  petroleum  hydrocarbons  as  high  as  ~9  mg  g'1 .  Polycyclic  aromatic 
hydrocarbons  (PAHs)  constitute  only  a  minor  fraction  of  these  residues  with  maximum 
concentrations  of  134  (ig  g"1,  but  their  fate  is  of  interest  because  of  their  potential  toxicity 
to  organisms.  Compared  to  typical  unweathered  No.  2  fuel  oil,  the  current  distribution  of 
PAHs  in  the  sediments  reflects  substantial  weathering  by  abiotic  and  biotic  processes; 
specifically  a  preferential  loss  of  naphthalenes  relative  to  phenanthrenes,  as  well  as 
isomer-specific  biodegradation  of  alkylated  PAHs.  Based  on  comparison  to  results  from 
an  earlier  study,  it  appears  that  little  or  no  change  has  occurred  to  the  distribution  of 
PAHs  since  1989,  indicating  that  weathering  at  this  site  has  stalled  or  is  now  proceeding 
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at  a  significantly  slower  rate.  To  assess  whether  sediment-water  partitioning  and 
molecular  diffusion  in  the  interstitial  medium  are  now  the  dominant  processes  controlling 
the  vertical  distribution  of  PAHs,  downcore  profiles  were  compared  to  a  numerical 
model.  Whilst  in  some  cases  the  model  accurately  reproduced  the  measured  data,  there 
were  instances  where  the  distribution  of  PAHs  was  slightly  under  or  overestimated. 
Reasons  for  these  discrepancies  are  discussed  and  are  likely  due  to  either  bioturbation, 
colloid-facilitated  transport  or  both.  Assessment  of  the  influence  of  these  processes  on  the 
spilled  oil  expands  our  understanding  of  the  overall  fate  of  these  compounds  and  their 
potential  long-term  effects  on  the  environment. 

2.1  INTRODUCTION 

A  major  effort  to  understand  the  fate  and  long-term  persistence  of  oil  in  the 
coastal  marine  environment  has  been  directed  at  the  sediments  from  the  Wild  Harbor  salt 
marsh.  West  Falmouth,  MA.  This  site  was  contaminated  when  the  barge  Florida  ran 
aground  and  released  -700,000  L  of  No.  2  fuel  oil  into  Buzzards  Bay  on  September  16, 
1969  (Sanders  et  al.,  1980).  Previous  work  at  the  M-l  site  in  Wild  Harbor,  which  was 
heavily  impacted  by  the  spill  and  the  subject  of  numerous  historical  studies,  revealed  that 
the  concentration  of  total  petroleum  hydrocarbons  (TPH)  in  these  sediments  has  not 
changed  substantially  in  the  past  30  years  (Reddy  et  al.,  2002).  In  addition,  the  use  of 
comprehensive  two-dimensional  gas  chromatography  (GCxGC)  provided  a  detailed 
assessment  of  the  composition  of  petroleum  residues  that  continue  to  persist  (Frysinger  et 
al.,  2003).  These  studies  have  revealed  that  distinct  groups  of  branched  alkanes, 
cycloalkanes,  and  polycyclic  aromatic  hydrocarbons  (PAHs)  remain  whereas  compounds 
typically  abundant  in  No.  2  fuel  oil,  such  as  the  H-alkanes,  are  not  detectable  (Reddy  et 
al.,  2002;  Frysinger  et  al.,  2003).  A  more  recent  survey  of  TPH  in  Wild  Harbor  has  found 
that  only  ~100  kg  continues  to  persist  in  -0.001  km2  of  salt  marsh  (Peacock  et  al.,  2005). 
Examination  of  the  natural  abundance  of  radiocarbon  within  different  pools  of 
sedimentary  organic  matter  has  also  shown  that  the  majority  of  this  oil  has  not  been 
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incorporated  in  the  insoluble  organic  matter  of  the  sediment  and  hence  likely  remains 
available  to  microorganisms  and/or  subject  to  physical  processes  such  as  diffusion  (White 
et  al.,  2005).  Natural  radiocarbon  abundances  within  phospholipids  isolated  from 
contaminated  sediment  intervals  at  site  M-l,  however,  indicate  that  the  oil  is  not  currently 
being  degraded  by  the  resident  microbial  population  (Slater  et  al.,  2005).  The  latter  is 
consistent  with  the  findings  of  Reddy  et  al.  (2002)  who  suggested  that  microbial 
degradation  of  the  branched  alkanes  and  other  constituents  of  the  spilled  oil  may  have 
ceased  over  25  years  ago.  Physical  processes,  therefore,  likely  play  the  most  important 
role  in  determining  the  fate  of  the  residual  oil  at  this  time. 

Marshes  are  considered  to  be  the  shorelines  most  capable  of  long-term  oil 
preservation  (Peacock  et  al.,  2005  and  references  therein)  and  thus  it  is  critical  that  we 
fully  understand  the  processes  that  control  the  fate  of  petroleum  in  these  environments.  In 
addition,  such  information  is  essential  to  damage  assessors  who  must  determine  whether 
to  restore,  rehabilitate,  or  replace  oil-impacted  areas  as  per  the  United  States  Oil  Pollution 
Act  of  1990  (33  U.S.C.  2701-2761).  Although  minor  constituents  of  petroleum,  PAHs  are 
of  particular  interest  due  to  their  carcinogenic,  mutagenic,  and  toxic  properties  (Office  of 
Health  and  Environmental  Assessment,  1985). 

Analysis  of  PAHs  in  surface  sediments  from  Wild  Harbor  following  the  spill 
revealed  a  decrease  in  concentration  of  the  lower  molecular  weight  PAHs  as  the 
dominance  of  the  Ci-naphthalenes  gave  way  to  C2-naphthalenes  by  January  1971,  C3- 
naphthalenes  by  November  1972,  and  C4-naphthalenes  by  September  1975  (Teal  et  al., 
1978).  Twenty  years  after  the  spill,  PAHs  continued  to  be  observed  in  the  sediment  at 
significant  concentrations.  For  example,  2,3,6-trimethylnaphthalene  and  3,6- 
dimethylphenanthrene  were  detected  at  concentrations  greater  than  3  and  5  (ig  g"1, 
respectively  (Teal  et  al.,  1992).  Molecular-specific  differential  weathering  of  PAHs  has 
been  observed  at  other  sites  (Wang  and  Fingas,  1995  and  2000),  and  general  trends  show 
that  degradation  rates  are  proportional  to  the  number  of  aromatic  rings  and  alkyl 
substitutions.  The  concentration  of  the  lower  molecular  weight  and  less  alkylated 
compounds  typically  decreases  over  time.  This  is  in  part  due  to  higher  solubilities  and 
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vapor  pressures  that  increase  the  likelihood  of  removal  by  processes  such  as  evaporation, 
water  solubilization,  diffusion,  and  microbial  degradation.  The  degradation  potential  of 
PAHs  also  differs  within  isomeric  groups,  for  example,  isomers  of  C2-naphthalene  have 
different  susceptibilities  to  microbial  degradation  depending  on  the  substitution  patterns 
(Wang  and  Fingas,  1995;  Farrington  et  al.,  1982).  Examination  of  individual  PAHs 
persisting  in  marine  sediments  at  the  West  Falmouth  site  provides  an  opportunity  to  study 
a  single  pulse  addition  of  oil-derived  PAHs  to  an  otherwise  pristine  system.  While  other 
sites  from  the  salt  marsh  contain  PAHs  with  vertical  distributions  similar  to  those  detailed 
here  (Peacock  et  al.,  2005),  this  study  focuses  on  the  M-l  site  because  of  the  extensive 
historical  data  available  (Teal  et  al.,  1978  and  1992)  as  well  as  the  more  recent 
geochemical  studies  (Reddy  et  al.,  2002;  Frysinger  et  al.,  2003;  Peacock  et  al.,  2005; 
White  et  al.,  2005;  Slater  et  al.,  2005).  In  addition,  it  was  the  most  contaminated  site 
when  Teal  et  al.  (1992)  surveyed  five  sites  in  and  around  Wild  Harbor  in  1989. 
Unfortunately,  the  original  Florida  oil  is  no  longer  available  and  hence  absolute 
differences  between  the  original  oil  and  current  residues  are  not  possible.  Teal  and 
colleagues  also  analyzed  sediment  horizons  at  a  much  coarser  resolution  (5  cm)  and  did 
not  quantify  the  alkylated  PAHs.  For  samples  collected  after  the  spill,  through  the  1970s 
(Teal  et  al.,  1978)  and  in  1989  (Teal  et  al.,  1992),  extracted  ion  chromatograms  for 
alkylated  naphthalenes  and  phenanthrenes  are  only  available  for  qualitative  comparisons. 
Here,  we  present  the  results  from  the  analysis  of  PAHs  from  one  sediment  core  collected 
in  2000  and  discuss  the  abundance,  distribution,  and  vertical  transport  of  these 
contaminants. 

2.2  MATERIALS  AND  METHODS 

2.2.1  Sediment  Collection 

In  August  2000,  a  15 -cm  diameter  core  was  collected  from  the  intertidal  sediment 
at  the  M-l  station  in  Wild  Harbor,  West  Falmouth,  MA.  The  core  was  returned  to  the 
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laboratory,  sectioned  at  2-cm  intervals  and  frozen  until  further  analysis  as  previously 
described  (Reddy  et  al.,  2002). 

2.2.2  Radiochronology 

One  to  two  grams  of  air-dried,  homogenized  sediment  was  counted  for  gamma 

210  214  137 

decay  of  Pb,  Pb  and  Cs  using  a  high  purity  germanium  detector  (Canberra  model 

210  214 

GCW  4023 S)  with  a  closed-end  coaxial  well.  Counting  errors  for  total  Pb  and  Pb 

210 

were  propagated  to  unsupported  Pb  and  ranged  from  5  to  26%  of  the  activity  (±  la), 
with  a  mean  error  of  12%  or  0.58  disintegrations  per  minute  (dpm)  per  gram  of  dry 

.  137 

weight  sediment.  Counting  errors  for  Cs  ranged  from  12  to  30%  of  the  total  activity  (± 

la),  with  a  mean  error  of  7%  or  0.03  dpm  per  gram  of  dry  weight  sediment.  An 

approximate  sedimentation  rate  of  0.35  cm  yr 1  was  estimated  from  the  activity  of 
210 

unsupported  Pb  (Figure  2.SIa)  using  the  constant  initial  concentration  (CIC)  model 

137 

(Eakins,  1983).  This  rate  is  consistent  with  a  maximum  in  Cs  activity  at  a  depth 
corresponding  to  1963  (Figure  2.SIb),  as  well  as  a  coincident  maximum  in  radiocarbon 
(White  et  al.,  2005 both  observed  as  a  result  of  the  peak  in  nuclear  weapons  testing.  The 
oil  spill  itself,  as  well  as  environmental  effects  such  as  erosion,  storm  events,  and 
bioturbation,  have  had  a  significant  effect  on  the  site,  hindering  the  establishment  of  a 
precise  chronology  for  the  core.  Potentially,  the  most  significant  effect  on  sediment 
accumulation  occurred  after  the  oil  spill  when  the  Spartina  alterniflora  was  killed  and  did 
not  grow  back  until  1975  (Reddy  et  al.,  2002). 

2.2.3  Sediment  Characteristics 

Roots  of  the  marsh  grass  Spartina  alterniflora  that  grows  at  the  M-l  site  were 
observed  throughout  the  length  of  the  core.  Moisture  contents  were  determined  as  the 
difference  between  sample  weights  before  and  after  drying.  The  average  solid-to-water 
ratio  of  the  sediment  bed  was  calculated  to  be  0.28  g  ml'1  (ranging  from  0.22  to  0.35  g 
ml'1).  Total  organic  carbon  (TOC)  was  calculated  from  dried  and  homogenized  samples 
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as  described  in  Xu  et  al.  (2001)  and  ranged  from  8-13%  with  an  average  of  1 1%  (Figure 
2.SIc). 

2.2.4  Petroleum  Residue  Analysis 

One  to  five  grams  of  air-dried  sediment  was  spiked  with  the  perdeuterated  parent 
PAHs,  fluorene-dio,  phenanthrene-^o,  pyrene-cfo,  chrysene-<ii2  and  benzo[6]- 
fluoranthene-^12  (0. 1-0.6  pg  of  each),  and  extracted  with  dichloromethane  and  methanol 
(9:1)  by  pressurized  fluid  extraction  (100°C;  1000  psi).  Extracts  were  reduced  in  volume, 
percolated  through  a  copper  column  to  remove  elemental  sulfur,  solvent  exchanged  into 
hexane,  and  charged  onto  a  glass  column  (20  cm  x  0.6  cm)  packed  with  fully  activated 
silica  gel  (100-200  mesh).  The  column  was  eluted  with  ten  milliliters  of  an  equal  mixture 
of  hexane  and  dichloromethane.  Half  of  the  extract  was  removed  and  spiked  with  9,10- 
dihydrophenanthrene  (0. 1-0.6  pg).  The  samples  were  analyzed  on  an  Agilent  6890  series 
gas  chromatograph  with  an  Agilent  5973  mass  selective  detector  (MSD).  The  MSD  was 
operated  in  the  selected-ion  monitoring  (SIM)  mode  for  quantitation  of  target  PAHs 
shown  in  Table  2.1.  Compounds  were  separated  on  a  J&W  DB-XLB  capillary  column 
(60  m,  0.25  mm  ID,  0.25  pm  film)  with  helium  carrier  gas  at  a  constant  flow  of  1  ml 
min  '.  The  GC  oven  had  an  initial  temperature  of  50  °C  (1  minute  hold)  and  was  ramped 
at  10  °C  min"1  until  320  °C  (37  minute  hold).  Parent  PAHs  were  quantified  using  the 
mass  and  area  of  the  corresponding  internal  standard  and  a  linear  term  of  a  six-point 
calibration  curve  generated  with  standard  solutions  prepared  from  the  National  Institute 
of  Standards  and  Technology  (NIST)  Standard  Reference  Material  2260  “ Aromatic 
Hydrocarbons  in  Toluene”.  For  quantification  of  the  alkylated  PAH  series,  the  response 
factors  of  the  parent  PAHs  were  used.  We  recognize  that  extracted  ion  chromatograms 
for  alkylated  phenanthrenes  include  alkylated  anthracenes,  and  that  extracted  ion 
chromatograms  for  Ci-pyrenes  include  Ci-fluoranthenes  and  benzofluorenes,  but  the 
concentrations  of  these  compounds  is  relatively  small  and  thus  for  simplicity  we  only 
refer  to  the  former  in  this  manuscript. 
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Table  2.1  Target  PAHs  analyzed  and  their  concentrations  (pg  g'1)  in  0-2,  14-16  and  28- 
BO  cm  sediment  horizons. 


PAH 

Symbol 

0-2  cm 

14-16  cm 

28-30  cm 

Naphthalene  *  a 

N„ 

nd 

nd 

nd 

2-Methylnaphthalene 

2MN 

nd 

0.013 

nd 

1  -Methylnaphthalene 

1MN 

nd 

nd 

nd 

Ci-naphthalenes  * 

N, 

nd 

0.023 

nd 

C2-naphthalenes  * §i 

n2 

nd 

0.318 

0.011 

C3-naphthalenes  * § 

n3 

0.050 

18.9 

0.081 

C4-naphthalenes  * § 

n4 

0.102 

51.3 

0.058 

Biphenyl  * 

BIP 

nd 

nd 

nd 

2,6-Dimethylnaphthalene 

DMN 

nd 

0.027 

nd 

Acenaphthylene  * 

ANY 

nd 

nd 

nd 

Acenaphthene  * 

ANE 

nd 

0.221 

nd 

2,3,5-Trimethylnaphthalene 

TMN 

nd 

nd 

nd 

Fluorene  * 

F0 

nd 

nd 

nd 

Dibenzothiophene  * 

Do 

nd 

0.438 

nd 

Phenanthrene  * 

Po 

0.066 

0.492 

0.027 

Ci-phenanthrenes  * § 

Pi 

0.056 

8.76 

0.022 

C2-phenanthrenes  * § 

P2 

0.076 

27.9 

0.035 

C3-phenanthrenes  * § 

P3 

0.035 

17.5 

0.025 

C4-phenanthrenes  * 

P4 

nd 

4.68 

0.018 

1  -Methylphenanthrene 

IMP 

nd 

2.17 

nd 

1 ,7-Dimethylphenanthrene 

DMP 

nd 

nd 

nd 

Anthracene  * 

ANT 

nd 

nd 

nd 

Fluoranthene  * 

FLT 

0.178 

0.376 

0.068 

Pyrene  * 

PYRo 

0.124 

0.856 

0.047 

Ci -pyrene  * 

PYRi 

0.050 

1.69 

0.030 

Benzo[a]anthracene  * 

BAA 

0.025 

0.035 

0.015 

Chrysene  * 

CHR 

0.066 

0.089 

0.018 

Benzo  [  b]  fluoranthene  * 

BBF 

0.098 

0.064 

0.030 

Benzo  [A:]  fluoranthene  * 

BKF 

0.028 

0.019 

nd 

Benzo[e]pyrene  * 

BEP 

0.058 

0.036 

0.015 

Benzo[a]pyrene  * 

BAP 

0.023 

0.020 

nd 

Perylene  * 

PER 

nd 

nd 

nd 

Dibenzo[a,  h]  anthracene  * 

DANT 

nd 

nd 

nd 

Indeno[l,23-cd]pyrene  * 

INP 

nd 

nd 

nd 

Benzo[g/n']perylene  * 

BPER 

nd 

nd 

nd 

Coronene  * 

COR 

nd 

nd 

nd 

EPAHs 

1.05 

134 

0.539 

a  XPAHs  is  the  sum  of  the  PAHs  marked  with  a 

b  Concentrations  above  which  adverse  effects  frequently  occur,  Effects  Range-Median  (ERM) 
(Long  et  al.,  1995,)  for  low  molecular  weight  PAHs  is  3.16  pg  g'1.  Here,  the  sum  of  C2-,  C3-,  C4- 
naphthalenes  and  Cr,  C2-,  C3-phenanthrenes  (marked  with  a  §)  is  125  pg  g'1  in  the  most 
contaminated  horizon  (14  to  16  cm). 
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2.2.5  Quality  Control 

Recoveries  of  the  deuterated  internal  standards  ranged  from  50-70%  with  an 
average  of  60%.  The  accuracy  of  the  method  was  assessed  by  analysis  of  the  NIST 
Standard  Reference  Material  1944  “New  York/New  Jersey  Waterway  Sediment’ ’  and 
values  within  70-130%  of  the  certified  values  were  obtained.  The  method  detection  limit 
(MDL)  for  PAHs  was  10  ng  g'1  and  procedural  blanks  were  always  at  or  below  the  MDL. 
The  precision  for  multiple  analyses,  as  expressed  by  the  relative  standard  deviation,  was 
less  than  10%. 

2.3  RESULTS  AND  DISCUSSION 

2.3.1  General  Trends  of  PAHs  in  Sediments 

PAHs  from  the  Florida  spill  continue  to  persist  in  salt  marsh  sediments  in  Wild 
Harbor  with  concentrations  of  EPAHs  as  high  as  134  pg  g'1  at  14  to  16  cm  (Table  2.1), 
coincident  with  the  maximum  concentration  of  TPHs  (~9  mg  g1,  Figure  2.1a).  PAHs 
were  detected  in  each  horizon  of  the  sediment  core  but  levels  were  most  elevated  between 
8  and  20  cm  (Figure  2.1b).  The  small  sub-surface  maximum  observed  at  24  to  28  cm 
most  likely  represents  spilled  oil  that  penetrated  through  channels  in  the  sediment. 
Sediments  at  the  M-l  site  are  very  heterogeneous  due  to  the  presence  of  these  channels 
and  the  extensive  Spartina  alterniflora  root  systems.  Analysis  of  sediments  a  few  years 
after  the  spill  noted  significant  penetration  of  oil  to  depths  of  70  cm  and,  in  some  cases, 
to  depths  greater  than  1  meter  (Bums  and  Teal,  1979).  Concentrations  of  individual  or 
groups  of  alkylated  PAHs  that  are  detectable  in  the  sediments  are  relatively  low  (0.01  - 
0.2  pg  g'1)  in  the  uppermost  (0-2  cm;  Figure  2.2a)  and  deepest  (28-30  cm;  Figure  2.2c) 
intervals  compared  to  significantly  higher  concentrations  (10  -  50  pg  g"1)  at  14-16  cm 
(Figure  2.2b).  The  distribution  of  PAHs  in  the  14-16  cm  horizon  is  typical  of  weathered 
No.  2  fuel  oil.  Specifically,  parent  PAHs,  naphthalene  and  phenanthrene  are  absent,  with 
a  predominance  of  C3-  and  C4-naphthalenes  and  C2-  and  C3-phenanthrenes  compared  to 
an  average  composition  of  four  different  No.  2  fuel  oils  (Figure  2.2d)  (Reddy,  1997).  This 
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Figure  2.2.  PAH  signatures  of  sediment  extracts  from  (a)  0-2  cm,  (b)  14-16  cm,  (c)  28-30 
cm  compared  to  (d)  the  relative  abundance  of  PAHs  present  in  an  average  composition  of 
four  No.  2  fuel  oils  (Reddy,  1997).  Abbreviations  for  PAHs  are  detailed  in  Table  2.1. 
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depletion  of  the  less  alkylated  compounds  and  continuing  dominance  of  the  latter 
compounds  is  supported  by  qualitative  observations  made  seven  (Teal  et  al.,  1978)  and 
20  years  (Teal  et  al.,  1992)  after  the  spill  and  in  other  studies  (Wang  and  Fingas,  1995; 
Wang  et  al.,  2001;  Ho  et  al.,  1999).  Sources  of  these  PAHs  can  be  confirmed  by  ratios  of 
specific  PAHs  such  as  fluoranthene  to  pyrene  and  Ci-phenanthrene  to  phenanthrene. 
Fluoranthene  is  thermodynamically  more  stable  than  pyrene  and  thus  it  is  possible  to 
differentiate  between  pyrogenic  (fluoranthene  :  pyrene  >1)  and  petrogenic  (fluoranthene: 
pyrene  <  1)  sources  (Sicre  et  al.,  1987).  At  the  M-l  site,  the  top  8  cm  and  bottom  18 
through  30  cm  (excluding  the  sub-surface  maximum  at  24  to  28  cm)  have  fluoranthene  to 
pyrene  ratios  ranging  from  1.1  to  1.5,  indicative  of  a  pyrogenic  source  most  likely  from 
atmospheric  deposition  of  combustion-derived  material.  The  middle  horizons  (8  to  18 
cm)  and  the  sub-surface  maximum  (24  to  28  cm),  however,  exhibit  a  petrogenic  signature 
with  fluoranthene  to  pyrene  ratios  ranging  from  0.4  to  0.9.  This  is  consistent  with  No.  2 
fuel  oil  that  has  a  fluoranthene  to  pyrene  ratio  of  0.6  to  0.9  (Gschwend  and  Hites,  1981). 
Because  unsubstituted  PAHs  are  the  most  abundant  compounds  in  combustion-derived 
sources  and  petroleum  has  an  increased  abundance  of  alkylated  compounds  (Gustafsson 
and  Gschwend,  1997),  the  Ci -phenanthrene  to  phenanthrene  ratio  can  also  be  examined. 
The  downcore  profile  (Figure  2.1c)  displays  Ci -phenanthrene  to  phenanthrene  ratios 
ranging  from  10-30  in  the  sediments  heavily  impacted  by  oil  (10  to  18  cm  and  24  to  28 
cm)  and  ~1  elsewhere. 

The  concentrations  of  PAHs  persisting  at  the  M-l  site  were  compared  to 
guidelines  for  contaminant  levels  below  or  above  which  adverse  effects  rarely  and 
frequently  occur,  the  Effects  Range-Low  (ERL),  and  Effects  Range-Median  (ERM), 
respectively  (Long  et  al.,  1995).  The  sum  of  C2-,  C3-,  C4-naphthalenes  andCi-,  C2-,  C3- 
phenanthrenes  approaches  125  |ig  g'1  in  the  most  contaminated  horizon  (14  to  16  cm), 
which  is  40  times  greater  the  ERM  value  for  low  molecular  weight  PAHs  (3.16  pg  g'1). 
In  fact,  all  sediment  horizons  with  a  clear  petrogenic  source  (Ci -phenanthrene  to 
phenanthrene  ratios  >10)  in  this  core  exceed  the  ERM,  thus  indicating  the  potential  to 
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adversely  affect  benthic  infauna.  Efforts  are  underway  to  more  accurately  gauge  the 
biological  effects  of  the  petroleum  residues  that  persist  at  this  site. 

Previous  studies  have  revealed  that  both  aerobic  and  anaerobic  biodegradation  of 
alkylated  PAHs  leaves  distinct  patterns  in  their  relative  abundances  due  to  the  preferential 
loss  of  certain  isomers.  For  example,  2-methylphenanthrene  is  more  easily  biodegraded 
than  other  methylphenanthrenes  (Wang  and  Fingas,  1995;  Huang  et  al.,  2004).  Extracted 
ion  chromatograms  for  the  Ci-,  C2-,  C3-naphthalenes  and  Ci-phenanthrenes  from  oil- 
contaminated  sediment  horizons  (8  to  18  cm),  which  are  not  shown,  exhibit  isomer 
specific  biodegradation  similar  to  trends  observed  by  Wang  and  Fingas  (1995).  Again, 
we  do  not  have  a  sample  of  the  original  Florida  oil  and  cannot  make  absolute 
comparisons.  It  is  also  not  possible  to  determine  exactly  when  biodegradation  occurred, 
but  the  distribution  of  Ci-phenanthrenes  observed  in  this  study  (collected  in  2000)  appear 
to  be  similar  to  those  shown  for  sediments  collected  in  1989  (Teal  et  al.,  1992).  This 
suggests  that  biodegradation  of  PAHs  is  no  longer  a  dominant  process,  thus  leading  us  to 
more  closely  investigate  the  role  of  vertical  transport  by  examining  downcore  profiles  of 
individual  PAHs. 

2.3.2  Downcore  PAH  Profiles 

The  downcore  profile  of  EPAHs  at  the  M-l  site  has  a  significantly  broader  peak 
than  that  of  the  TPH  (Figure  2.1).  The  peak  width,  which  is  determined  by  measuring  the 
width  of  the  peak  base,  is  16  cm  for  EPAHs  and  12  cm  for  the  TPH.  Further  examination 
of  downcore  profiles  of  alkylated  naphthalenes  and  phenanthrenes  (the  most  abundant 
PAHs),  indicate  that  variations  in  vertical  distribution  are  also  evident  at  the  molecular 
level  (Figure  2.3),  for  example,  the  peak  maximum  for  Ci -naphthalene  (19  cm)  is  shifted 
downward  compared  to  the  TPH  maximum  (15  cm)  and  is  considerably  broader  with  a 
width  of  16  cm.  The  C2-  and  C3-naphthalenes,  however,  have  peak  maxima  concurrent 
with  the  TPH,  and  the  C3-naphthalenes  exhibit  a  more  symmetric  profile.  The 
phenanthrenes  show  broader  peaks  (~14  cm)  than  the  TPH,  but  narrower  than  the 
naphthalenes  and  have  peak  maxima  that  coincide  with  the  TPH.  Peak  widths  of  the 
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Figure  2.3  Downcore  profiles  of  TPH  (Reddy  et  al.,  2002)  (a  and  e),  alkylated  naphthalenes, 
Nl5  N2,  N3  (b,  c,  d,  respectively),  and  akylated  phenanthrenes,  Pl9  P2,  P3  (f,  g,  h,  respectively). 
Abbreviations  for  PAHs  are  explained  in  Table  1. 


phenanthrenes  broaden  slightly,  but  consistently  as  alkyl  substituents  are  added,  which  is 
the  reverse  of  that  observed  for  the  naphthalenes. 

The  general  observations  of  the  peak  widths  for  the  TPH,  alkylated  naphthalenes, 
and  phenanthrenes  can  be  explained  by  the  diffusivities  of  these  compounds  in  the  bulk 
sediment.  The  effective  diffusivity  coefficient  (Dieff)  is  calculated  for  each  compound 
assuming  a  two-phase  system  of  sediment  and  interstitial  water,  including  the  effects  of 
porosity,  tortuosity,  and  the  solid-water  partition  coefficient  (K<j)  (Table  2.2).  The 
naphthalenes  have  the  greatest  Djeff  and  hence  diffuse  the  farthest  and  exhibit  the  broadest 
peaks  (Figure  2.3).  The  narrowing  of  the  peak  widths  as  alkyl  substituents  are  added  to 
the  naphthalene  structure  is  reflected  in  a  similar  decrease  in  the  values  of  Djeff  (Table 
2.2).  The  Djeff  of  the  phenanthrenes,  is  an  order  of  magnitude  lower  than  the  naphthalenes 
(Table  2.2)  and  this  is  evidenced  by  their  narrower  peak  widths  (Figure  2.3).  The  Dieff  for 
the  phenanthrenes  also  decreases  as  alkyl  substituents  are  added,  but  unlike  the 
naphthalenes,  this  is  not  reflected  in  the  measured  data  and  a  slight  broadening  of  the 
peak  width  is  observed  as  alkyl  substituents  are  added.  To  increase  the  diffusion  of  these 
more  hydrophobic  compounds,  a  third  phase  is  considered.  The  presence  of  porewater 
colloids  has  previously  been  observed  to  enhance  “dissolved”  concentrations  of 
compounds  beyond  that  predicted  by  their  solubilities  (Chiou  et  al.,  1998;  Brownawell, 
1986)  and  concentrations  of  porewater  dissolved  organic  carbon  (DOC)  in  a  Spartina 
alterniflora  salt  marsh  have  been  observed  to  range  from  0.4  to  1.9  mM  (Gribsholt  et  al., 
2003).  A  significant  portion  (5 1  ±  26%)  of  DOC  in  freshwater  samples  is  known  to  be  of 
high  molecular  weight  (Repeta  et  al.,  2002)  and  would  thus  serve  as  a  good  sorbent  that 
could  facilitate  diffusion.  The  more  alkylated  and  therefore  hydrophobic  the  compound, 
the  greater  this  contribution  would  be,  potentially  explaining  the  increase  in  observed 
peak  width  with  increase  in  alkylation  for  the  phenanthrenes. 

The  overall  increased  dispersion  of  the  PAHs  compared  to  the  TPH  highlights 
their  increased  partitioning  into  water  and  ability  to  diffuse  over  a  wider  area.  The  TPH  at 
this  site  is  represented  by  an  unresolved  complex  mixture  (UCM)  that  is  composed 
mainly  of  saturated  hydrocarbons  (85%)  (Frysinger  et  al.,  2003).  The  latter  have 
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considerably  greater  log  KoW  values  (e.g.  dodecane:  log  KoW  ~  7  (MacKay  et  al.,  1992) 
than  the  PAHs  (Table  2.2)  and  as  a  result  of  their  increased  hydrophobicity,  are  unlikely 
to  have  undergone  significant  diffusion  in  the  past  30  years.  The  relatively  narrow  peak 
width  likely  represents  the  initial  oil-covered  sediment  distribution.  At  its  highest 
concentration  (11-17  cm,  Figure  2.1a),  the  peak  width  of  the  TPH  is  ~6  cm  suggesting 
that  when  the  oil  was  spilled,  it  penetrated  to  approximately  this  depth  at  this  site.  In 
general,  the  peak  widths  of  the  PAHs  can  be  explained  by  their  diffusivities,  but  to 
determine  whether  or  not  this  component  controls  their  distribution,  the  calculated  Djeff 
for  each  compound  and  a  simple  model  was  constructed. 

2.3.3  Mathematical  Modeling  of  Postdepositionai  Transport  of  PAHs 

The  processes  controlling  the  vertical  distribution  of  PAHs  in  the  marsh  sediment 
were  examined  by  comparison  of  the  measured  data  to  a  one-dimensional  model.  This 
approach  was  developed  previously  by  Reiser  et  al.  (1997)  to  study 
alkylbenzenesulfonates  in  lake  sediments  and  is  based  on  a  simplified  version  of  the 
diagenetic  equation  (Berner  et  al.,  1989)  having  the  form: 

(dC/dt)  =  Dieff  (d2C/dx2)  +  co  -(dC/dx)  (I) 

where  C  is  the  concentration  of  the  solute  in  mass  per  unit  volume  of  pore  water  solution, 
Dieff  is  the  effective  diffusivity  coefficient  in  the  bulk  sediment  of  the  compound  of 
interest  (Table  2.2),  co  is  the  sediment  accumulation  rate,  i.e.  the  velocity  of  sediment 
flow  relative  to  the  sediment-water  interface,  t  is  time,  and  x  is  depth  in  the  sedimentary 
column.  The  sediment-water  interface  was  chosen  as  the  origin  of  the  depth  coordinate. 
The  initial  input  of  oil  was  set  at  a  depth  of  0  through  6  cm  allowing  for  settling  of  oil  in 
the  top  layers  of  sediment  as  previously  discussed.  The  velocity  of  burial  of  pore  water 
was  set  equal  to  the  velocity  of  particle  burial  thereby  neglecting  the  advective  flow  of 
pore  water  relative  to  sediment  particles.  The  velocity  of  burial  was  set  constant  at  0.35 
cm  y'!  based  on  sediment  deposited  between  -1910  and  2000  and  as  there  was  no 
significant  evidence  for  mixing  processes  from  the  210Pb  profile  (Figure  2.SIa),  only 
molecular  diffusion  was  assumed  for  Diefr,  thus  excluding  bioturbation  affects.  Any  loss 
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Table  2.2  Properties  of  PAHs  calculated  to  determine  the  effective  diffusivity  (Diefr)  of 
each  PAH. 


PAH 

log  K,)W° 

logK„c0 

Kac 

fj 

Ci -Naphthalene 

3.94 

3.54 

350 

1.0E-02 

5.7E-06 

8.1E-01 

C2-Naphthalene 

4.15 

3.75 

560 

6.4E-03 

5.3E-06 

4.7E-01 

C3-Naphthalene 

4.36 

3.95 

900 

4.0E-03 

5.0E-06 

2.8E-01 

Ci-Phenanthrene 

4.99 

4.57 

3720 

9.6E-04 

4.9E-06 

6.5E-02 

C2-Phenanthrene 

5.20 

4.78 

5970 

6.0E-04 

4.6E-06 

3.8E-02 

C3-Phenanthrene 

5.40 

4.97 

9380 

3.8E-04 

4.4E-06 

2.3E-02 

a  Octanol-water  partition  coefficients  (log  K<,w)  were  calculated  with  an  equation  from  Eastcott  et 
al.  (1988)  based  on  jthe  number  of  carbons  in  each  compound.  These  values  compare  well  with 
expected  and  measured  values  previously  described  (MacKay  et  al.,  1992;  Ruelle,  2000). 
Measured  values  are  not  available  for  all  compounds  examined  in  this  study  and  thus  the 
calculated  values  were  used  to  maintain  consistency  between  compounds. 

*  Log  values  of  the  organic  matter-water  partition  coefficient  (KoC)  were  calculated  from  Log  KoW 
using  a  linear  free-energy  relationship  for  aromatic  hydrocarbons  (Schwarzenbach  et  al.,  2003 
and  references  therein),  where  log  =  0.98  *  log  Kow  -  0.32. 

c  Solid-water  partition  coefficient  (K<j)  was  calculated  by  multiplying  KoCby  the  weight  fraction  of 
solid  which  is  organic  carbon  (foc)  (Schwarzenbach  et  al.,  2003),  which  at  the  M-l  site  is  an 
average  of  0.1  kg  organic  carbon  kg'1  solid  (Figure  2.SIc). 

d  By  combining  Kj,  and  the  solid-to-water  ratio  of  the  sediment  (rsw)  discussed  in  the  text,  the 
fraction  of  the  PAH  dissolved  is  determined,  fiw  =  (1  +  r^Kd)'1. 

e  The  diffusivity  in  water  (DjW)  is  calculated  by  estimating  the  size  (Vj)  of  the  PAHs  using  the 
approach  of  Fuller  et  al.,  (1966)  as  well  as  determining  the  solution  viscosity  at  the  temperature  of 
interest  (q),  which  in  this  case  is  1.139.  Combining  this,  we  can  solve  the  following  equation 
from  Hayduck  and  Laudie  (1974)  where,  Diw  =  (13.26  *  10'5)  /  (q114  *  V;"0'58). 
f  The  effective  diffusivity  (Dieff  )  is  calculated  by  multiplying  Diw  by  fiw  and  dividing  by  the 
tortuosity  (2.27) 


due  to  biodegradation  was  assumed  negligible  when  compared  to  the  concentrations 
present  in  the  sediments  and  considering  that  significant  biodegradation  is  not  occurring 
today  (Slater  et  al.,  2005).  The  validity  of  this  assumption  and  the  role  that 
biodegradation  may  have  played  since  the  oil  was  spilled  is  revisited  in  the  discussion. 
This  model  considers  a  two  phase  system  of  sediment  and  interstitial  water,  but  it  has 
been  shown  that  the  oil  itself  should  be  considered  as  a  third  phase  as  it  may  affect  the 
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partitioning  of  PAHs  depending  on  its  type,  concentration,  and  age  or  weathering  state 
(Jonker  et  al.,  2003).  Sorption  enhancing  effects,  however,  are  only  observed  at  high  oil 
concentrations  above  15%  (w/w)  on  a  sedimentary  organic  carbon  basis.  At  the  M-l  site, 
the  oil  comprises  less  than  10%  (w/w)  of  the  total  organic  carbon  at  all  depth  horizons 
(Figure  2.SIc)  and  thus  this  third  phase  is  not  incorporated  into  the  model.  The  Crank- 
Nicholson  implicit  method  was  employed  to  determine  concentrations  of  Ci-,  C2-,  03- 
naphthalene  and  Cr,  C2-,  C3-phenanthrene  based  on  equation  (I)  as  a  function  of  depth. 
This  method  estimates  spatial  derivatives  at  both  the  present  and  a  future  time  step  and 
then  uses  these  to  obtain  a  spatial  estimate  that  corresponds  to  the  midpoint  of  the  time 
step  (Crank,  1975).  No  attempts  were  made  to  model  the  sub-surface  maximum  at  24-26 
cm.  It  is  worth  noting  that  the  salt  marsh  sediment  in  this  study  is  considerably  more 
heterogeneous  than  the  lake  sediment  examined  by  Reiser  et  al.  (1997),  and  thus  it  is  of 
interest  to  assess  the  utility  of  this  method  for  a  different  sediment  matrix. 

The  downcore  profiles  of  the  alkylated  PAHs  and  the  fit  to  the  corresponding  model  data 
(Figure  2.4)  can  be  summarized  into  three  main  groups.  The  first  consists  of  the  Ci-  and 
C2  -naphthalenes  that  show  broad,  asymmetric  profiles  that  are  overestimated  by  the 
model  in  the  upper  portion  of  the  profile  (Figure  2.4a  and  b).  These  are  the  only  PAHs  to 
exhibit  concentrations  in  the  upper  5cm  that  are  significantly  greater  than  the  detection 
limits  (Figure  2.3b  and  c).  The  model  predicts  the  presence  of  Cj-  and  C2-  naphthalenes 
in  the  upper  5cm  because  their  diffusivities  (0.81  and  0.47  cm2  y'!  respectively)  exceed 
the  rate  of  burial  (0.35  cm  y"1),  but  the  failure  to  correctly  predict  their  relative 
contributions  is  most  likely  due  to  a  loss  term  in  the  surface  sediments  that  is  not 
considered  by  the  model.  This  loss  term  may  include  diffusion  into  the  overlying  water, 
evaporation  or  biodegradation  that  is  more  likely  to  occur  in  shallow  sediments  where  the 
availabilities  of  electron  acceptors  (such  as  oxygen,  nitrate  or  sulfate)  is  high  and  where 
microbial  populations  are  greater  (Slater  et  al.,  2005).  Since  prior  work  indicates  that 
biodegradation  is  no  longer  significant  (Slater  et  al.,  2005),  the  major  loss  may  be  due  to 
evaporation  or  diffusion.  Biodegradation  may  still  be  important,  however,  if  it  is  a 
seasonal  process  that  was  not  actively  occurring  when  the  work  by  Slater  et  al.,  (2005) 
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was  performed.  The  second  group  consists  of  the  C3-  naphthalenes  and  Ci  - 
phenanthrenes  (Figure  2.4c  and  d)  that  display  average  width,  symmetrical  profiles  well 
represented  by  the  model  that  therefore  appear  to  be  controlled  by  diffusion  only.  Slight 
discrepancies  between  the  model  and  measured  data  are  observed,  but  most  likely  arise 
from  an  overestimation  of  the  calculated  diffusivity.  The  third  group  of  compounds  are 
the  least  diffusive  and  most  hydrophobic,  and  consist  of  the  C2  and  C3  -phenanthrenes 
(Figure  2.4e  and  f).  These  compounds  have  relatively  narrow  peak  widths  that  fit  the 
model  well,  but  have  asymmetric  profiles,  displaying  a  sharp  decrease  in  concentration 
below  the  peak  maximum,  but  a  much  shallower  decrease  above.  This  apparent 
“smearing”  of  the  concentration  profile  may  be  a  result  of  bioturbation  that  only  occurs 
above  the  peak  maximum.  Bioturbation  would  affect  all  compounds  present  in  the  oil,  but 
the  smearing  is  not  observed  for  all.  The  magnitude  of  bioturbation  may  be  too  small  to 
be  noticeable  for  the  more  diffusive  PAHs  and  thus  is  only  observed  for  the  least 
diffusive  compounds,  such  as  the  C2-  and  C3-  phenanthrenes.  This  hypothesis  is 
supported  by  the  TPH  profile,  which  has  a  similar  smearing  and  as  discussed  previously, 
is  considered  to  be  mainly  composed  of  hydrophobic  compounds  with  low  diffUsi vities . 
Understanding  past  variations  in  bioturbation  since  the  spill,  and  possible  reasons  for  its 
occurrence  above  the  peak  maximum,  but  not  below,  are  important  considerations.  In  the 
years  following  the  spill,  bioturbation  decreased  significantly  due  to  the  eradication  of  the 
macrobenthos  observed  twelve  hours  after  the  spill,  with  the  only  slight  recovery  in 
infauna  populations  six  years  later  (Krebs  and  Bums,  1977).  Subsequently,  as  the  marsh 
began  to  recover  and  recolonization  took  place,  bioturbation  could  again  become  a  factor. 
The  macrobenthos,  however,  may  not  have  penetrated  into  or  through  the  sediment 
containing  the  highest  concentration  of  the  oil  due  to  its  inherent  toxicity  (as  previously 
discussed),  and  thus  the  bottom  half  of  the  PAH  profile  would  remain  intact,  as  appears 
to  be  the  case  (Figure  2.4e-f).  The  effects  of  pollution  on  the  mixing  intensities  of 
macrobenthos  have  previously  been  examined  (Wheatcroft  and  Martin,  1996).  While  it  is 
difficult  to  constrain  this  factor  in  the  heterogeneous  sediments  examined  here,  it  seems 
likely  that  such  effects  played  a  role  in  determining  the  profiles  of  the  PAHs  that  we 
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P,  normalized  concentration  P2  normalized  concentration 


P3  normalized  concentration 


Figure  2.4  Downcore  profiles  of  alkylated  naphthalenes;  Ni,  N2,  N3  (a,  b,  c  respectively, 
solid  line)  and  alkylated  phenanthrenes;  Pi,  P2,  P3  (d,  e,  f  respectively,  solid  line) 
compared  to  model  data  (dashed  line).  All  concentrations  are  normalized  to  their 
maximum  concentration.  Abbreviations  for  PAHs  are  explained  in  Table  2.1. 


observe  today. 

Considering  these  different  profile  types,  it  is  clear  that  diffusivity  is  important, 
but  that  other  factors  should  be  considered.  The  discrepancies  between  the  measured  and 
model  profiles  highlight  the  complexities  involved  in  understanding  the  post  depositional 


fate  of  organic  compounds.  It  appears  that  a  combination  of  diffusion  in  the  porewater, 
sorption  onto  colloidal  matter,  bioturbation  and  biodegradation  contribute  to  the 
distribution  of  PAHs  observed  in  the  sediments  today.  Given  the  rates  of  diffusion  for  the 
PAHs  (Table  2.2)  and  the  sedimentation  rate  at  the  site  (0.35  cm  y'1),  only  the  Ci-  and  C2 
-naphthalenes  are  likely  to  diffuse  out  of  the  sediment  and  into  the  overlying  water.  In 
summary,  this  study  provides  a  quantitative  evaluation  of  some  of  the  most  toxic  and 
persistent  compounds  present  in  the  environment  today.  The  role  that  these  post- 
depositional  physical  processes  may  exert  on  organic  contaminants  in  marine  sediments 
is  not,  however,  limited  to  PAHs  (Achman  et  al.,  1996)  and  the  findings  of  this  study 
should  be  considered  when  predicting  the  long  term  fate  of  a  variety  of  organic 
contaminants  as  well  as  when  using  sediments  to  reconstruct  histories  of  environmental 
contamination. 
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deviation  from  triplicate  analyses). 
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ABSTRACT 

To  provide  a  new  perspective  on  the  fate  of  petroleum  in  the  marine  environment, 
we  utilized  variations  in  the  natural  abundance  of  radiocarbon  (14C)  to  detect  and  quantify 
petroleum  residues  that  have  persisted  in  Wild  Harbor  sediments,  West  Falmouth,  MA, 
for  more  than  thirty  years.  The  5730  yr  half-life  of  14C  makes  this  isotope  ideal  for  the 
detection  of  fossil  fuel-derived  contaminants  (14C  free)  within  different  fractions  of 
natural  organic  matter  (modem  14C  content)  in  environmental  matrices.  Samples  of  both 
contaminated  and  uncontaminated  sediments  were  sequentially  treated,  first  by  solvent 
extraction,  followed  by  saponification,  and  then  acid  hydrolysis.  Radiocarbon  analysis  of 
the  sediment  residues  and  select  extracts  was  performed  to  probe  for  the  presence  of 
fossil  fuel  contaminants  and/or  their  metabolites  in  different  pools  of  sedimentary  organic 
matter.  Our  results  indicate  that  the  majority  of  fossil  carbon  is  solvent  extractable  and 
has  not  been  incorporated  in  the  insoluble  organic  matter  in  sediment.  Unextracted 
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sediments  contaminated  with  petroleum  contain  significantly  less  14C  than  extracted 
sediments,  and  isotope  mass  balance  calculations  suggest  that  up  to  ~9%  of  the  total 
organic  carbon  (TOC)  in  the  petroleum  contaminated  sediment  horizons  is  derived  from 
solvent-extractable  petroleum.  These  estimates  are  similar  to  values  calculated  when  the 
total  quantities  of  oil  (measured  by  gas  chromatography  with  flame  ionization  detector 
(GC-FED))  are  compared  to  TOC  content  (determined  by  elemental  analysis).  These 
results  pave  the  way  for  applications  of  this  isotopic  approach  to  more  complex 
environmental  systems  where  the  fate  of  contaminants  is  less  certain. 

3.1  INTRODUCTION 

The  fate  of  hydrophobic  organic  contaminants  (HOCs)  in  the  aquatic  environment 
is  largely  controlled  by  their  low  aqueous  solubility  and  high  affinity  for  particles, 
rendering  sediments  their  primary  sink.  Once  incorporated  into  sediments,  HOCs  are 
subject  to  varied  and  complex  processes  that  confound  our  understanding  of  their 
ultimate  fate.  For  example,  laboratory  (Richnow  et  al.,  1997,  1998,  2000;  Hatcher  et  al., 
1993)  and  field  studies  (Beller  et  al.,  1998;  Richnow  et  al.,  1994;  Reddy  et  al.,  2000; 
Schwarzbauer  et  al.,  2003)  have  shown  that  certain  HOCs  (or  their  metabolites)  can  bind 
to  natural  organic  matter,  but  this  process  is  not  easy  to  follow  and  is  especially  difficult 
to  quantify  using  current  analytical  approaches.  These  associations  are,  however,  of 
interest  as  they  can  potentially  affect  the  bioavailability,  toxicity,  transport  and  long-term 
fate  of  HOCs  (Alexander,  1995).  The  aforementioned  studies  define  HOCs  in  sediment  as 
either  “free”  or  “unbound”  and  therefore  solvent  extractable,  or  “non-extractable”  or 
“bound”  and  released  only  by  harsher  chemical  treatments  such  as  base  or  acid 
hydrolysis. 

Defining  and  quantifying  free  and  bound  pools  is  critical  when  assessing  the 
potential  impact  that  contaminants  may  have  on  the  environment  and  human  health.  In 
soils,  it  is  proposed  that  irreversible  formation  of  bound  residues  could  be  a  major  sink 
for  certain  organic  contaminants  (Dec  and  Bollag,  1997).  Promotion  of  this  process  has 


52 


thus  been  considered  as  a  bioremediation  strategy.  It  is  not  clear  how  robust  this  approach 
would  be  when  applied  to  aquatic  sediments  as  opposed  to  soils,  as  there  are  significant 
differences  between  the  two  environments.  Such  differences  include  water  content, 
oxygen  availability,  pH  and  microbial  population  (Hedges  et  al.,  1997)  as  well  as  the 
presence  of  reactive  species  such  as  hydroxyl  radicals  that  may  affect  potential  reactions 
between  HOCs  and  organic  matter.  Processes  affecting  the  overall  fate  of  HOCs  in 
sediments  must  therefore  be  examined  independently. 

The  use  of  natural  abundance  14C  measurements  is  a  novel  approach  to  trace 
HOCs  into  various  pools  of  sedimentary  organic  matter  (Reddy  et  al.,  2002a).  The 
absence  of  14C  in  fossil  fuel-derived  contaminants  serves  as  a  natural  isotopic  label  (i.e. 
an  inverse  tracer)  that  is  intrinsic  to  the  contaminant.  Thus  its  presence  and  provenance  in 
environmental  samples  can  be  studied  directly.  This  is  desirable,  as  model  studies 
involving  additions  of  labeled  compounds  (i.e.,  laboratory  incubations  or  even  additions 
in  situ)  may  not  accurately  reflect  processes  occurring  in  the  field.  Work  by  Steinberg  et 
al.  (1987),  for  example,  examined  the  degradation  of  14C-labelled  ethylene  dibromide 
(EDB)  in  EDB  contaminated  soils  and  noted  that  the  added  label  was  metabolized  faster 
than  the  native  contaminant.  Using  natural  abundance  14C  measurements,  it  is  also 
possible  to  quantitatively  assess  contributions  of  fossil  carbon  through  the  principle  of 
isotope  mass  balance.  Variations  in  14C  have  previously  been  used  to  constrain 
contributions  of  fossil  hydrocarbons  to  marine  sediments  (Zafiriou,  1973;  Farrington  and 
Tripp,  1977;  Giger  et  al.,  1984;  Wakeham  et  al.,  2004).  In  these  studies,  significant 
depletions  in  14C  confirm  the  presence  of  a  fossil  component  to  the  sedimentary  organic 
matter. 

In  the  present  study,  sediments  from  Wild  Harbor,  West  Falmouth,  MA,  a  site 
contaminated  with  No.  2  fuel  oil  from  the  Florida  oil  spill  of  1969,  were  examined  to 
attempt  to  constrain  contributions  from  fossil-fuel  derived  compounds  to  the  sediment. 
Since  the  spill,  this  site  has  been  the  subject  of  extensive  and  long-term  studies  (Reddy  et 
al.,  2002 b  and  references  therein)  and  provides  a  well  characterized  site  at  which  to 
evaluate  the  utility  of  these  molecular  and  isotopic  approaches. 
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3.2  MATERIALS  AND  METHODS 

3.2.1  Samples 

In  August  2000,  a  15-cm  diameter  sediment  core  was  collected  from  the  intertidal 
zone  at  the  M-l  station  in  Wild  Harbor,  West  Falmouth,  MA  as  previously  described 
(Reddy  et  al.,  20026).  The  core  was  extruded  and  sectioned  at  2-cm  intervals  from  0  to  36 
cm.  The  sections  were  placed  in  precombusted  glass  jars  and  frozen  until  further  analysis. 
Ten  to  twenty  grams  of  each  sediment  horizon  was  defrosted  and  wet  sieved  through  a 
1.4  mm  sieve  to  remove  the  intact  roots  of  the  marsh  grass  Spartina  altemiflora  that 
grows  at  the  M-l  site.  The  sieved  sediment  was  left  to  air  dry  before  homogenizing  with 
a  mortar  and  pestle.  The  sequential  chemical  treatment  of  these  sediments  is  summarized 
in  Figure  3.1. 

■a 

3.2.2  Sediment  Chronology 

One  to  two  grams  of  air-dried,  homogenized  sediment  was  counted  for  gamma 

210  214  137 

decay  of  Pb,  Pb  and  Cs  using  a  high  purity  germanium  detector  (Canberra  model 
GCW  4023S)  with  a  closed-end  coaxial  well.  An  approximate  sedimentation  rate  of  0.35 
cm  yr 1  was  estimated  from  the  activity  of  unsupported  21°Pb  using  the  constant  initial 
concentration  (CIC)  model  (Eakins,  1983).  This  rate  is  consistent  with  a  maximum  in 
137Cs  activity  at  a  depth  corresponding  to  1963,  which  was  the  peak  in  nuclear  weapons 
testing  activities  that  occurred  during  the  late  fifties  and  early  sixties.  Obtaining  a  precise 
chronology  of  this  core  is  difficult  because  of  physical  effects  such  as  erosion,  storm 
events  and  bioturbation.  The  oil  spill  also  had  a  significant  effect  on  the  site,  killing  the 
Spartina  altemiflora.  The  latter  did  not  grow  back  until  1975  with  an  attendant  influence 
on  sediment  accumulation  (Reddy  et  al.,  20026). 

3.2.3  Solvent  Extraction 

One  to  five  grams  of  air-dried  and  homogenized  sediment  were  extracted  with 
dichloromethane  (DCM)  and  methanol  (MeOH)  (9:1)  by  pressurized  fluid  extraction  (100 


54 


°C,  1000  psi).  The  resulting  total  lipid  extracts  (TLEs)  were  reduced  in  volume  by  rotary 
evaporation  and  percolated  through  a  copper  column  to  remove  sulfur  prior  to  further 
analysis. 


SEDIMENT  (TOC) 


Pressurized  fluid  extraction 
(100°C,  1000  psi) 


Alkaline 
hydrolysis 
KOH/MeOH 
reflux  2  h 


Ether  cleavage 
57%  HI  reflux  1  h, 
LiAIH4  reflux  1 .5  h 


r 

1 

r 

1 

ACID  HYDROLSIS 
EXTRACT  (AH-EX) 

ACID  HYDROLYZED 
RESIDUE  (AH-RES) 

Figure  3.1  Flow  diagram  to  illustrate  the  sequential  treatment  of  sediments  and  residues. 
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3.2.4  Saponification 

Saponification  was  performed  to  cleave  ester-linked  moieties  from  the  solvent- 
extracted  sediments,  as  in  (Richnow  et  al.,  1994).  A  solution  of  0.5  N  potassium 
hydroxide  (KOH)  (in  MeOH,  100  mL)  and  water  (H2O)  (20  mL)  were  added  to  one  to 
five  grams  of  pre-extracted  sediments  from  horizons  0-2,  8-10,  10-12,  12-14,  14-16,  30- 
32  and  32-34  cm  and  refluxed  for  2  h.  Extracted  sediment  samples  from  30-32  and  32-34 
cm  were  combined  to  obtain  the  desired  quantity  of  sediment.  Approximately  67%  of  the 
total  pooled  sample  with  a  TOC  content  of  9.8%  came  from  the  30-32  cm  horizon  and  the 
remainder  with  a  TOC  content  of  8.7%  came  from  the  32-34  cm  horizon.  Extracted 
sediments  from  12-14  (11.1%  TOC  content)  and  14-16  cm  (10.9%  TOC  content)  were 
processed  separately  as  well  as  a  mixture  comprised  of  equal  quantities  of  the  two 
horizons  to  obtain  enough  material  for  further  chemical  treatment  post  saponification. 
After  cooling,  the  reaction  mixture  was  separated  by  centrifugation  and  the  remaining 
sediment  rinsed  with  MeOH,  followed  by  DCM,  and  then  hexane.  The  alkaline  solution 
was  acidified  with  4N  hydrochloric  acid  (HC1)  to  pH  2  and  extracted  with  DCM.  The 
combined  extracts  were  reduced  in  volume  by  rotary  evaporation  and  dried  over  Na2S04. 

3.2.5  Acid  Hydrolysis 

Hydriodic  acid  (HI)  has  previously  been  used  to  cleave  ether  bound  compounds  in 
sedimentary  samples  (Hoefs  et  al.,  1997;  Hold  et  al.,  1998;  Schouten  et  al.,  2000)  and 
was  used  here  to  release  any  petroleum  residues  bound  by  ether  linkages.  One  to  two 
grams  of  previously  saponified  sediment  residues  from  0-2,  10-12,  12-16  and  30-34  cm 
horizons  were  each  added  to  57%  HI  solution  (in  H20,  50  mL)  and  refluxed  under  a 
nitrogen  atmosphere  for  1  h.  After  cooling,  the  sediment  residue  was  removed  by 
centrifugation  and  rinsed  with  H2O,  followed  by  hexane  and  allowed  to  air  dry.  The 
reaction  solution  was  extracted  with  hexane  and  the  combined  extracts  washed  with  H2O 
saturated  with  sodium  chloride  (NaCl)  to  remove  any  excess  acid,  followed  by  5% 
sodium  thiosulfate  (Na2S2(>3)  solution  (to  remove  iodine).  The  extract  was  reduced  in 
volume  by  rotary  evaporation  before  being  charged  onto  a  glass  column  (20  cm  x  0.6  cm) 
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packed  with  fully  activated  silica  gel  (100-200  mesh).  The  column  was  eluted  with 
hexane  and  DCM  (9:1)  and  the  extract  air-dried  before  being  taken  up  in  1,4-dioxane. 
Lithium  aluminum  hydride  (LiAlFL)  (600  mg)  was  added  to  the  extract,  which  was  then 
refluxed  under  nitrogen  for  1.5  h  to  reduce  any  alkyliodides  to  hydrocarbons.  After 
cooling  to  room  temperature,  ethyl  acetate  was  slowly  added  to  deactivate  the  remaining 
LiAIHU.  The  reaction  mixture  was  then  centrifuged  and  the  recovered  supernatant  washed 
with  H20  saturated  with  NaCl  followed  by  H20  alone,  before  being  reduced  in  volume 
and  dried  over  Na2SC>4. 

3.2.6  Gas  Chromatography-Mass  Spectrometry  (GC/MS)  Analysis 

TLEs,  saponified  extracts  and  extracts  from  acid  hydrolysis  were  derivatized  with 
bis-(trimethylsilyl)-trifluoroacetamide  (BSTFA)  and  pyridine  to  improve  chromatography 
by  silylating  active  hydrogens  present  in  the  samples.  For  compound  identification, 
electron  ionization  (El)  spectra  were  acquired  on  an  Agilent  6890  series  gas 
chromatograph  with  an  Agilent  5973  mass  selective  detector  (MSD).  Compounds  were 
separated  on  a  J&W  DB-5MS  capillary  column  (60  m  x  0.32mm  ID,  0.25  |Xm  film 
thickness)  with  helium  carrier  gas  at  a  constant  flow  of  1  ml  min"1.  The  GC  oven  had  an 
initial  temperature  of  40°C  (1  minute  hold)  and  was  ramped  at  20°C  min"1  until  130°C 
and  then  at  2°C  min"1  until  320°C  (30  minute  hold).  Spectra  were  acquired  between  m/z 
40-650  at  a  scan  rate  of  1  cycle  s'1.  Compounds  present  in  the  extracts  were  either 
tentatively  identified  from  mass  spectra  and  gas  chromatographic  retention  characteristics 
or  identified  by  comparison  to  pure  standards. 

3.2.7  Bulk  Sediment  and  Extract  Analyses 

Aliquots  of  the  dried,  homogenized  untreated,  extracted,  saponified  and  acid 
hydrolyzed  sediment  residues  were  treated  with  10%  HC1  to  remove  any  carbonates.  The 
sediments  were  then  analyzed  for  organic  carbon  content  with  a  Carlo  Erba  1108 
elemental  analyzer.  The  stable  carbon  isotope  ratio  (813C)  and  radiocarbon  abundance 
(A14C)  of  the  sediments  as  well  as  select  TLEs  and  saponified  extracts  (representing  both 
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the  contaminated  region  of  the  core  and  the  uncontaminated  areas  above  and  below)  was 
measured  on  purified  carbon  dioxide  (CO2)  after  their  combustion  in  the  presence  of 
cupric  oxide  (CuO).  Stable  carbon  values  were  determined  by  isotope  ratio  mass 
spectrometry  and  14C  content  by  accelerator  mass  spectrometry  (AMS)  at  the  National 
Ocean  Sciences  Accelerator  Mass  Spectrometry  (NOSAMS)  facility  at  Woods  Hole 
Oceanographic  Institution  (WHOI),  Woods  Hole,  MA,  after  conversion  of  the  CO2  to 
graphite  (McNichol  et  al.,  1994).  All  14C  measurements  are  normalized  to  813C  values  of 
-25%o  and  expressed  as  A14C  values,  which  is  the  per  mille  (%o)  deviation  from  the 
international  standard  for  14C  dating,  Standard  Reference  Material  4990B  “Oxalic  Acid”. 
Routine  precision  for  813C  and  A14C  measurements  are  ~0.1  and  10%o,  respectively.  All 
radiocarbon  abundances  obtained  in  this  study  were  corrected  for  in-situ  14C  decay  using 
A14C  =  [fme  <1950^-1]1000  (Stuiver  and  Polach,  1977),  where  X= 1/8267  (year '),  fm- 
fraction  modem  14C,  corrected  for  isotopic  fractionation  using  813C,  and  “x”  equals  the 
year  of  deposition  (determined  from  the  2,0Pb  chronology.  Table  2).  This  yields  the  A14C 
values  at  the  time  of  sediment  deposition. 

3.3  RESULTS  AND  DISCUSSION 

3.3.1  Vertical  distribution  of  petroleum  hydrocarbon  contaminants  in  salt  marsh 
sediments 

A  recent  study  examining  sediments  from  the  Wild  Harbor  site  shows  that  in 
discrete  locations  within  the  salt  marsh,  oil  residues  persist  more  than  30  years  after  the 
original  oil  was  spilled  (Reddy  et  al.,  20026).  This  study  provided  a  refined  inventory  of 
the  petroleum  hydrocarbons  that  have  been  removed  from  the  sediments,  such  as  the  n- 
alkanes,  as  well  as  those  that  continue  to  persist.  The  latter  include,  branched  alkanes, 
cycloalkanes  and  aromatic  hydrocarbons.  The  persistence  of  these  compounds  is  apparent 
when  examining  the  downcore  profile  of  total  petroleum  hydrocarbons  (TPH)  in  marsh 
sediments  (Figure  3.2a).  A  maximum  of  ~  8-9  mg  g'1  is  observed  in  the  12-16  cm  depth 
interval,  which,  based  on  2,0Pb  and  137Cs  profiles  corresponds  well  with  the  timing  of  the 
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spill.  Minimum  values  of  ~  0-2  mg  g'1  are  observed  throughout  the  rest  of  the  core  (0-12 
cm  and  16-34  cm).  Other  significant  observations  from  the  downcore  profiles  include  the 
A14C  values  of  the  sediments,  residues  and  extracts  from  the  sequential  chemical 
treatments.  These  display  very  similar  profiles  with  a  maximum  in  14C  levels  at  a  depth  of 
8-10  cm,  corresponding  to  the  well  established  “bomb  14C  spike”  resulting  from 
aboveground  thermonuclear  weapons  testing  in  the  late  1950s  and  early  1960s  (Levin  et 
al.,  1985)  (Figure  3.2b,  c,  d,  e).  Significant  offsets  are  observed,  however,  between  the 
different  sediment  residues  and  these  coincide  with  the  maximum  TPH  content.  The 
implications  of  these  isotopic  variations  will  be  discussed  further. 

3.3.2  Radiocarbon  as  an  Inverse  Tracer 

Radiocarbon  abundances  (A14C)  of  the  organic  carbon  from  untreated  sediments, 
residues  and  corresponding  extracts  were  measured  from  horizons  throughout  the 
sediment  core  (Table  3.1  and  Figure  3.2).  In  the  most  heavily  contaminated  horizons  (12- 
lb  cm)  (Figure  3.2a),  a  significant  offset  is  observed  between  the  unextracted  and  solvent 
extracted  sediments  (indicated  by  arrows  in  Figure  3.2b).  This  is  due  to  the  influence  of 
,4C-free  petroleum  from  the  oil  spill.  The  unextracted  sediments  contain  significantly  less 
14C  than  the  extracted  sediments.  Correspondingly,  the  14C  content  of  the  TLE  is 
significantly  lower  than  that  from  the  uncontaminated  horizons  (Figure  3.2c).  The  TLEs 
from  the  contaminated  horizons  of  the  core  are  comprised  of  a  mixture  of  fossil  carbon 
(A14C  =  -1000%o)  and  biogenic  hydrocarbons  (14C  ~  +100%o)  and  this  is  reflected  in  the 
A14C  values  of  these  extracts.  This  is  supported  by  GC/MS  analysis  of  TLEs  throughout 
the  entire  core,  which  indicate  the  presence  of  biogenic  hydrocarbons  including  higher 
plant  wax-derived  w-alkanes  from  n-C 25  to  /1-C33  with  a  strong  odd/even  dominance, 
alkanoic  acids  from  n-C\e  to  n-C 34  with  an  even/odd  dominance,  n-alkanols  from  H-C22  to 
H-C36  with  even/odd  dominance,  and  the  C29  sterols  P-sitosterol  and  stigmasterol,  which 
are  major  sterols  in  higher  plants.  The  higher  plant  triterpenoid  lupeol  was  also  observed. 
The  contaminated  sections  (8-16  cm)  are,  however,  dominated  by  an  unresolved  complex 
mixture  (UCM)  of  hydrocarbons  in  the  boiling  range  of  n- C13  to  M-C24  alkanes,  consistent 
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Table  1.  Stable  carbon  and  radiocarbon  abundances  for  sediment  residues  and  extracts 


U  ** 
2  si 
<  w 


U 

<fr 

< 

U 

<] 

u 

< 


tO  w 


U  ^ 
3  ^ 

tO  w 


60 

X 

H 

t 


U  ^ 

r<]  £ 


u  ^ 
p.  si 
o  w 


4 

03 


1/3 

M 

c~~ 

*n 

VO 

tq 

vq 

tN 

§ 

vri 

00 

© 

cn 

in 

CN 

CO 

X 

w 

N" 

vo 

oo 

ov 

i 

U 

8 


bo 

X 

% 

CA 


|JH  3  '"g 

r-»  o\  H 

CO  w  H 


u  ^ 

5=4 


3 


X 

w 


U  * 
3  sl  8 
to  W  H 


fr  s 

<U  <J 
TJ  w 


o 

00 

r-; 

00 

vq 

co 

_  VO 

Ov 

cn 

o 

00 

oo 

vo 

CN 

OV 

N" 

ov 

cn 

f-H 

l-H 

tN 

r— 1 

Ot) 

vo 

cn 

i 

3  f' 
3  T 

«> 

■ 

irJ 

i 

cd 

o' 

VO 

1 

co 

CO 


VO 

On 


CN 

r-H 

VO 


CN 

•4 

<N 


CN 


t>  in 

CnI  I —  O  OV 


CN 

o 

co 


CO  vo  >o 
(N  VO  On  —I 
t"  tN  O  t- 
m  tn  in  vo 

■  i  i  i 


CN 

r-' 

CO 

to 

■ 


iq 

tN 


VO 

>o 

tN 

i 


•vt 

o 

fH 

Ov 

VO 

VO 

»— t 

On 

CO 

»H 

tN 

O 

o' 

CN 

OV 

rr 

vo 

vo 


co 

wi 

«o 


co  vq  tN  to  t— i 

O  •— i  >d  co  co 

vo  in  co  ^  ^ 

i  i  i  l  i 


00 


<N 

I 


«q 

CN 

tN 


UO 

rH 

tN 


Ov 

td 


q  vq  ©  vq 

O  On  CO  tN 
tN  1  CN  CN 


w> 

Ov 


oo  vq 

co  tN 
tN  tN 


to 

ON 


00 

■4 

tN 


I 

tN 


■'t 


co  cn  wo 

vd  vo  vd 

■  i  i  i  i 


i- 

r- 

00 


vo 


©  On  CO  t";  CO  00  00 

odr^odr^r^r^vdvd 

»-H  T-H  i—l  1-H  t-H 

I  I  I  I  I  I  I  ■ 


00  t-;  CO  — < 

vo  vo  vd  VO 

T— H  *— I  I  ^ 

I  I  I  I 


oo 

wo 


CN 

CN 

o 

1-H 

q 

Ov 

O) 

q 

-t: 

O 

rH 

r- 

N" 

o\ 

00 

00 

oo 

00 

00 

td 

td 

td 

6 

vd 

vd 

in 

c 

r—i 

1 

r-H 

1 

1 

1 

1 

1 

1 

r-H 

1 

t—H 

1 

^H 

1 

i 

r— H 

1 

■ 

»— H 

1 

CN 

vo 

VO 

00 

o 

CN 

CN 

Tf 

CN 

Tf 

vo 

oo 

'n- ' 
vr"H 

r-H 

vhH 

t-H 

r-H 

CN 

CN 

cn 

m 

cn 

o 

CN 

vo 

1 

OO 

O 

CN 

-4 

CN 

VO 

00 

O 

o 

CN 

6 

i— i 

r-H 

r-H 

t-H 

CN 

cn 

cn 

cn 

X 

<u 


S  .2 
“  pq 

T3  t 

tS  H 

,  c5^ 

■'$  ©“ 

<L>  -O 

©  8  *g 

W  CA  4_» 

^  | 

1  s  I 

O  <L> 

•Is  a  £ 

CA  Q  r-* 

±  .2  o  +j 

•s  it  § 

«S  -g  ^ 

x)  iu  3 

a>  +-»  u  Q 

tj  a  rj  o 

b  £  &  o 
8  §  .2  S 

2  2  w  o 

j.g3  i 

I?  I 

.a  §*  .2 

ll  1  S 

.2  “  03  <* 

w  g  ©  a 
2  -5  0  o 

>  ^  *8  I 

U  2  -d  ■§ 

2u-§  g 


61 


h  nm  refers  to  samples  that  were  not  measured. 

1  This  samples  is  a  mixture  of  sediment  horizons  12-14  and  14- 16cm. 
7  This  samples  is  a  mixture  of  sediment  horizons  30-32  and  32-34cm. 


with  a  biodegraded  No.2  fuel  oil  source  (Reddy  et  al.,  2002c;  Frysinger  et  al.,  2003). 

Comparison  of  the  14C  contents  of  the  solvent  extracted  and  saponified  sediments 
shows  little  difference,  with  the  exception  of  the  12-14  cm  horizon  (Figure  3.2d).  This 
shift  and  smaller  changes  downcore  indicate  less  14C  in  the  saponified  residue  than  the 
extracted  residue  from  which  it  originated.  This  suggests  that  the  compounds  liberated  by 
the  process  of  saponification  in  these  sediments  are  more  enriched  in  14C  than  the 
remaining  saponified  sediment  residue.  This  is  supported  by  the  14C  contents  of  the 
saponified  extracts  (Figure  3.2e),  which  display  a  modem  signature  reminiscent  of  the 
bomb  spike  unperturbed  by  petroleum  contamination.  Moreover,  based  on  GC/MS 
analysis,  the  saponified  extracts  were  similar  in  composition  throughout  the  sediment 
core  and  indicate  the  presence  of  n-alkanols  from  n-C^i  to  n- Cn  with  even/odd 
dominance,  alkanoic  acids  from  n-Cio  to  n-C2&  with  even/odd  dominance  and  sterols  (as 
before).  The  extracts  from  the  contaminated  zone  contained  no  evidence  of  additional 
components  that  could  derive  from  hydrocarbons. 

There  are  observable  offsets  in  A14C  between  the  saponified  and  hydriodic  acid 
hydrolyzed  sediment  residues  throughout  the  contaminated  and  uncontaminated  horizons 
of  the  core  (Figure  3.2f).  In  all  cases,  the  acid  hydrolyzed  residue  is  more  depleted  in  14C 
than  the  saponified  residue.  This  suggests  that  more  recently  synthesized  material  has 
been  removed  by  this  treatment.  The  remaining  residue  is  significantly  older  and  may  be 
partially  comprised  of  aliphatic  biopolymers  that  are  highly  refractory  and  not 
hydrolysable  by  these  chemical  treatments  (Hedges  et  al.,  2000).  In  the  contaminated 
horizon  (12-16  cm)  the  effect  of  the  bomb  spike  is  evident  by  a  more  14C  enriched  acid 
hydrolysis  sediment  residue  relative  to  the  pre-bomb  sediment  residue  in  the  30-34  cm 
horizon  (approximate  date  1909)  and  the  most  recently  deposited  post-bomb  material 
deposited  in  the  0-2  cm  horizon.  GC/MS  analyses  of  the  extracts  indicate  the  presence  of 
benzo[Z>]thiophene  and  its  alkylated  derivatives  as  well  as  alkylated  benzenes  in  all  of  the 
depth  horizons  examined.  These  compounds  are  unlikely  to  be  natural  products  and 
instead  may  arise  from  alteration  of  natural  compounds  during  the  hydriodic  acid 
hydrolysis  treatment.  The  extract  obtained  from  acid  hydrolysis  of  the  contaminated 


62 


section  of  the  core  (12-16  cm),  however,  also  displays  a  few  compounds  that  are 
exclusive  to  this  sample.  These  compounds  may  be  related  to  the  contamination  in  some 
way  either  directly  or  indirectly  as  a  result  of  the  oil  spill.  They  could,  for  example,  be 
due  to  a  different  natural  source  of  organic  matter  related  to  a  shift  in  the  biological 
community  structure  in  response  to  the  spill.  The  structure  of  these  compounds  is 
currently  unknown  and  their  identification  and  source  will  be  the  focus  of  future  work. 
For  the  present  investigation,  it  is  important  to  note  that  the  abundance  of  these 
compounds  is  insufficient  to  influence  the  14C  content  of  the  sediment  residues. 

3.3.3  Stable  Carbon  Isotopic  Compositions  of  Sediments,  Residues  and  Extracts 

Stable  carbon  isotopic  compositions  of  the  organic  carbon  from  untreated 
sediments,  as  well  as  corresponding  residues  and  extracts  were  measured  from  horizons 
throughout  the  sediment  core  (Table  3.1  and  Figure  3.3).  Patterns  observed  after  each 
sequential  chemical  treatment  are  similar  to  those  observed  for  A14C.  As  before,  an  offset 
is  observed  between  the  unextracted  and  extracted  sediments  (Figure  3.3b)  within  the 
contaminated  horizons  (Figure  3.3a),  which  we  attribute  to  the  presence  of  13C-depleted 
petroleum.  The  813C  value  of  the  remaining  oil  at  this  site  is  ~-27%o  (Reddy  et  al.,  2002c) 
and  thus  contaminated  sediments  are  shifted  from  the  813C  of  the  natural  organic  matter 
towards  this  end  member.  Organic  carbon  in  salt  marsh  sediments  containing  Spartina 
alterniflora  have  previously  been  shown  to  have  813C  values  between  -17  and  -24%o 
(Boshker  et  al.,  1999).  These  values  are  proposed  to  reflect  a  variety  of  carbon  inputs 

ii 

including  those  depleted  in  C  such  as  phytoplankton,  benthic  algae  or  C3  vascular  plant 
material  as  well  as  the  13C  enriched  Spartina  alterniflora  carbon.  The  813C  of  the 
uncontaminated  TOC  in  the  core  top  is  more  depleted  (~  -18%o)  than  the  uncontaminated 
sediments  lower  in  the  core  (-16%o).  This  2%o  shift  may  reflect  selective  preservation  of  a 
13C  enriched  biomass  source  such  as  Spartina  over  the  more  13C  depleted  phytoplankton 
biomass.  The  observed  shift  in  813C  could  alternatively  indicate  a  non-steady  state  of 
carbon  inputs  to  the  salt  marsh  over  time.  Middleburg  et  al.,  (1997)  suggest  that  an 
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S,  unfilled  diamonds).  Dashed  horizontal  lines  highlight  the  zone  of  maximum  TPH  contamination. 


important  mechanism  affecting  the  813C  of  sedimentary  TOC  in  salt  marshes  is  the  input 
of  allochthonous  organic  matter,  which  includes  organic  matter  sorbed  on  mineral  matter 
as  well  as  marine  phytoplankton,  microphytobenthos  and  non-local  macrophytes.  The 
TLE  is  also  quite  enriched  (Figure  3.3c),  and  in  the  contaminated  horizons  shows  more 
depleted  (-24%o)  values  closer  to  that  of  petroleum  (-27%o,  as  before),  which  is  to  be 
expected.  Very  little  difference  is  observed  between  the  extracted  and  saponified 
sediment  residues  (Figure  3.3d)  and  saponified  extract  (Figure  3.3e).  Differences  in  813C 
are,  however,  observed  between  the  saponified  and  the  acid  hydrolyzed  sediment  residues 
(Figure  3.3f).  In  all  cases  the  acid  hydrolyzed  residue  is  more  C  depleted  than  the 
saponified  residue,  and  this  difference  is  most  significant  (3.5%o)  in  the  0-2  cm  horizon. 
The  smallest  difference  of  ~0.5%o  is  observed  in  the  contaminated  horizon  12-16  cm.  The 
offsets  observed  between  these  two  residues  in  the  different  sediment  horizons  could 
arise  from  isotopic  fractionation  associated  with  slightly  different  chemical  conditions 
imposed  on  these  sediments  throughout  the  multi-step  chemical  treatment.  Alternatively, 
these  differences  in  8,3C  may  represent  different  sources  of  carbon  to  the  marsh  as 
previously  discussed.  As  mentioned,  the  acid  hydrolyzed  sediment  residues  may  contain 
refractory  aliphatic  biopolymers.  The  presence  of  these  biopolymers,  which  are 
synthesized  from  lipids  that  are  typically  depleted  relative  to  bulk  biomass,  may  help  to 
explain  the  more  depleted  813C  signal  of  the  acid  hydrolyzed  residues.  It  is  evident  from 
this  data,  however,  that  there  is  no  significant  manifestation  of  any  remaining  oil  on  the 
813C  values  of  the  residues. 

3.3.4  Isotope  Mass  Balance 

The  fraction  of  TOC  that  is  derived  from  petroleum  in  these  sediments  can  be 
calculated  based  on  isotope  mass  balance.  TOC  from  the  bulk,  unextracted  sediment  is 
considered  to  be  a  mixture  ( t )  comprised  of  exogenous  (x)  and  endogenous  («) 
components  that  can  be  represented  by 

(roc)  =  (: toc)x + (: mc)n  (i) 

and  the  isotopic  composition  of  the  mixture  is  given  by 
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Ft  =  fFx  +  {\-f)Fn  (II) 

where  the  Fs  represent  isotopic  abundances  and  /  is  the  fraction  of  the  mixture  that  is 
exogenous,  in  this  case,  derived  from  petroleum.  In  terms  of  radiocarbon  abundances  (A), 
we  can  write 

/  =  (At  —  An)/(Ax  —  An)  (III) 

The  corrected  A14C  values  for  the  untreated  sediments  were  used  as  At  and  the 
corrected  A14C  values  for  the  extracted  sediment  residues  were  used  for  An.  Assuming  Ax 
to  reflect  an  absence  of  14C  (i.e.,  A14C  =  -1000%o),  solutions  of /  were  obtained  for  the 
sediment  horizons  examined  in  this  study.  These  are  expressed  as  the  percentage  of  TOC 
derived  from  solvent-extractable  petroleum  and  are  shown  in  Table  3.2  alongside  values 
calculated  by  comparing  TPH  concentrations  (measured  by  GC-FID)  to  TOC  content 
(determined  by  elemental  analysis).  The  mass  of  TPH  was  converted  to  carbon 
equivalents  considering  the  mass  contribution  of  hydrogen  to  the  hydrocarbon,  which  in 
an  alkane  typical  of  the  UCM  (w-Cn  to  H-C24  as  before)  is  approximately  15%  of  the 
total.  Potential  errors  could  arise  depending  on  the  specific  contribution  of  hydrogen  to 
the  TPH,  which  could  range  from  15%  for  the  most  saturated  compounds  such  as  the 
alkanes  to  5%  for  the  most  unsaturated  compounds  such  as  the  3-ring  aromatic 
hydrocarbons.  If  a  5%  contribution  of  hydrogen  was  used  instead  of  15%,  the  %TOC 
derived  from  petroleum  values  in  Table  3.2  would  increase  by  no  greater  than  1%.  The 
reported  values  (Table  3.2)  show  good  agreement  throughout  the  core  indicating  a 
maximum  of  6-9%  TOC  derived  from  petroleum  in  the  contaminated  horizons  (12-16 
cm).  In  the  relatively  uncontaminated  horizons,  above  and  below  this,  ~  0-2%  of  TOC  is 
derived  from  the  petroleum.  There  is  a  difference  of  ~l-2%  between  the  two  methods  and 
this  is  likely  due  to  the  error  associated  with  the  14C  abundance  measurements  and/or  the 
mass  contribution  of  hydrogen  to  the  TPH  as  previously  discussed.  It  is  clear,  however, 
that  within  this  error,  the  fraction  of  petroleum  that  is  measurable  by  GC-FID  is 
comparable  with  the  isotope  mass  balance. 

These  results  indicate  that  at  the  bulk  level,  the  vast  majority  of  petroleum  appears 
to  be  solvent-extractable  and  not  bound  to  the  insoluble  organic  matter  in  the  sediment.  It 
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is  important  to  note  that  the  TLE  could  contain  soluble  esterified  contaminants.  However, 
since  the  free  contaminants  quantified  by  GC-FID  correspond  well  to  the  mass  balance, 
the  proportion  that  is  bound  must  be  extremely  small.  Given  that  the  contaminant  in  this 
study  consists  of  relatively  inert  hydrocarbons,  it  is  not  surprising  that  significant  binding 
is  not  observed.  For  binding  to  occur,  compounds  must  possess  functional  groups,  and 


Table  3.2  Elemental  and  Petroleum  concentrations  in  West  Falmouth  sediments,  2000 


depth 

(cm) 

approx 

age 

TOC 

(%) 

TPH 

(mgg'V 

%  TOC  from  petroleum  measured  by 
mass  balance  gas  chromatography  b 

0-2 

1997 

8.42 

0.00 

0.16 

0.00 

2-4 

1991 

9.78 

0.00 

-0.03 

0.00 

4-6 

1986 

12.0 

0.00 

-1.12 

0.00 

6-8 

1980 

12.8 

0.04 

-0.10 

0.03 

8-10 

1974 

10.8 

2.08 

2.71 

1.63 

10-12 

1969 

12.8 

2.51 

0.95 

1.67 

12-14 

1963 

11.6 

7.57 

8.25 

5.56 

14-16 

1957 

11.7 

8.95 

8.62 

6.52 

16-18 

1951 

11.2 

1.41 

-1.36 

1.07 

18-20 

1946 

10.5 

0.24 

2.77 

0.19 

20-22 

1940 

11.1 

0.44 

2.09 

0.34 

30-32 

1911 

10.0 

0.00 

2.06 

0.00 

32-34 

1906 

8.09 

0.00 

-0.25 

0.00 

0  Data  from  Reddy  et  al.,  (20026)  where  TPH  is  total  petroleum  hydrocarbons. 

*  To  calculate  the  percent  total  organic  carbon  (%  TOC)  derived  from  petroleum,  the  mass  of 
TPH  was  converted  to  carbon  equivalents  considering  a  mass  contribution  of  hydrogen  to  the 
hydrocarbons,  which  in  an  alkane  typical  of  the  UCM  (n-CI3  to  n-C2 4  as  before)  is  approximately 
15%  of  the  total. 
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whilst  petroleum  metabolites  resulting  from  microbial  degradation  are  chemically  more 
reactive,  we  found  no  detectable  evidence  of  any  binding  of  degradation  products.  This 
may  be  because  degradation  at  this  site  is  slow  or  has  ceased.  Recent  work  by  Slater  et 
al.,  (2005),  who  have  examined  the  14C  content  of  lipid  markers  of  microbial  activity 
present  in  the  sediments,  supports  this  claim.  The  study  by  Slater  et  al.  indicates  that  there 
is  currently  no  measurable  flow  of  petroleum  carbon  into  the  biomass  and  that  the 
microbial  population  are  preferentially  utilizing  the  natural  organic  matter  in  the  sediment 
as  a  carbon  source.  This  conclusion,  coupled  to  the  findings  of  this  study  indicating  the 
absence  of  significant  binding  or  incorporation  of  the  petroleum,  leaves  us  to  speculate  on 
the  reason(s)  why  hydrocarbon  contamination  persists  in  these  sediments  more  than  30 
years  after  the  spill.  It  has  been  proposed  that  this  persistence  is  due  to  the  lack  of  sulfate, 
a  necessary  electron  acceptor  for  anaerobic  degradation  of  petroleum  or  that  the  residual 
compounds  are  structurally  too  complex  for  anaerobic  degradation  or  are  inaccessible  to 
microbial  attack  (Reddy  et  al.,  2002 b).  Further  studies  will  explore  these  scenarios. 

3.3.5  Sensitivity  of  Isotopic  Approach 

While  it  is  clear  from  the  present  study  that  14C  can  be  used  as  an  intrinsic  tracer 
of  contaminants  in  environmental  matrices;  there  are  limitations  to  this  approach.  For 
example,  it  may  be  that  contaminants  are  binding  to  sediments  at  this  site,  but  not  at 
concentrations  observable  using  14C  as  an  inverse  tracer.  Figure  3.4  details  the  expected 
change  in  A14Ctoc  as  a  function  of  the  percentage  of  petroleum  present  in  the  sediment. 
Sediment  containing  0.1%  petroleum  represents  a  relatively  non-impacted  A14Ctoc  of 
+100%o,  a  value  chosen  to  represent  modem  sediments.  Previous  work  measuring  14C 
content  of  discrete  carbon  pools  in  Buzzards  Bay,  MA,  1991  yielded  A14C  values  of  eel 
grass  particulate  organic  carbon  (POC)  of  62%o,  and  Spartina  alterniflora  as  198%o 
(McNichol  et  al.,  1994).  We  assume  that  the  bulk  of  the  organic  matter  at  our  site  is 
Spartina  alterniflora  as  opposed  to  phytoplankton  and  therefore  choose  a  post-bomb, 
heavily  Spartina  influenced  value  of  100%o.  Potential  errors  inherent  in  assuming  this 
isotopic  end  member  value  are  shown  in  Figure  3.4.  This  figure  clearly  demonstrates  that 
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Figure  3.4.  The  sensitivity  of  A14C  of  total  organic  carbon  (A14CToc)  to  additions  of 
petroleum  (A14C=  -1000%o)  to  sediments  with  varying  A  14Cjoc  values.  Modem  sediments 
are  represented  by  a  A14Ctoc  of  +100%o.  Two  other  potential  end  members  were  chosen 
to  represent  the  particulate  organic  carbon  (POC)  of  die  salt  marsh.  These  were  Eel  grass 
and  Spartina  and  are  represented  by  A14Ctoc  values  of  +62%o  and  +198%o  respectively, 
as  in  McNichol  et  al.,  (1994). 


given  an  analytical  precision  of  10%o,  at  least  1%  of  petroleum  (A14C  =  -1000%o)  must  be 
present  to  significantly  affect  A14Ctoc-  Further  analysis  of  sediments  on  a  molecular 
rather  than  bulk  level  would  undoubtedly  resolve  petroleum-derived  compounds  that  are 
bound,  although  it  is  clear  from  the  isotope  mass  balance  that  such  findings  would  not  be 
quantitatively  significant.  It  is  also  possible  that  bound  compounds  of  petroleum  origin 
could  be  present  in  the  extracts  obtained  from  each  of  the  sequential  chemical  treatments, 
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as  mentioned  previously,  but  this  is  likely  an  extremely  small  quantity  as  the  percentage 
of  fossil  carbon  determined  by  mass  balance  compares  well  to  that  measured  by  GC-FED. 
Future  work  will  be  directed  at  applying  this  approach  to  sediments  and  soils  containing 
more  reactive  contaminants  to  determine  the  significance  of  binding,  as  well  as  molecular 
level  14C  studies  to  develop  a  more  detailed  understanding  of  the  provenance  and  fate  of 
contaminants  in  environmental  samples. 
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x  j 

prior  to  solvent  extraction, c  EX-RES  is  the  extracted  sediment  residue,  SA-RES  is  the  saponified  sediment  residue, e  AH- 
RES  is  the  acid  hydrolyzed  sediment  residue,  ^TLE  is  the  total  lipid  extract  obtained  from  solvent  extraction, 8  SA-EX  is  the 
saponification  extract. 


CHAPTER  4 


RELATIONSHIPS  BETWEEN  RADIOCARBON  CONTENT  AND  MODES  OF 
ASSOCIATION  OF  NATURAL  ORGANIC  MATTER  IN  COASTAL  AND 
CONTINENTAL  MARGIN  SEDIMENTS 


ABSTRACT 


We  have  investigated  the  relationships  between  radiocarbon  (14C)  and  stable 
carbon  (13C)  isotopic  composition  and  the  different  modes  of  association  of  organic 
matter  (OM)  present  in  surficial  sediments  from  near-shore  and  continental  margin  sites 
that  range  in  terms  of  inputs  and  depositional  conditions.  Our  aim  is  to  gain  insight  into 
the  sources  and  sinks  of  OM  in  the  marine  environment,  and  to  enhance  our 
understanding  of  the  nature  and  distribution  of  OM  in  sediments.  Radiocarbon  is  a 
particularly  sensitive  tracer  of  OM  provenance,  with  variations  in  its  natural  abundance 
reflecting  age  and  pre-  and  post-depositional  history.  In  order  to  investigate  isotopic 
variability  within  different  sedimentary  OM  pools,  sediments  were  sequentially  treated 
first  by  solvent  extraction  to  examine  non-covalently  associated  compounds  followed  by 
saponification  to  cleave  ester-linked  moieties.  Radiocarbon  analysis  was  performed  on 
the  bulk  sediment  and  resulting  sediment  residues.  The  molecular  composition  of  the 
extracts  was  examined  using  gas  chromatography/mass  spectrometry  (GC/MS),  and  the 
relative  contributions  of  terrestrial  and  marine  biomarkers  released  from  the  sequential 
chemical  treatments  were  assessed.  For  select  sites,  pre-industrial  sediment  horizons  were 
also  studied  to  examine  diagenetic  influences  on  the  relationships  between  age,  chemical 
association  and  reactivity  of  sedimentary  OM.  Radiocarbon  abundances  (A14C)  of  the 
bulk  sediment  reflect  a  mixture  of  modem,  pre-aged  and  fossil  (i.e.  14C-“dead”)  OM  of 
marine  and  terrestrial  origin  with  little  difference  in  A14C  between  the  unextracted  (TOC) 
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and  solvent  extracted  sediments  (EX-RES).  Distinct  offsets  in  A14C  between  EX-RES  and 
saponified  sediment  residues  (SA-RES)  observed  at  most  sites  indicate  the  removal  of 
more  14C  enriched  material  by  this  treatment.  Saponification  extracts  (SAE)  are  mainly 
composed  of  short-chain  (n-C \2  to  H-C24)  alkanoic  acids  with  an  even/odd  dominance 
indicating  a  predominantly  marine  algal  or  microbial  source.  These  observations  provide 
evidence  for  the  protection  of  labile  marine  carbon  by  chemical  binding.  Exceptions  to 
these  overall  patterns  exist  and  are  attributed  to  either  fossil  carbon  inputs  and/or  the 
impact  of  diagenetic  processes. 

4.1  INTRODUCTION 

The  burial  of  organic  matter  (OM)  in  recent  marine  sediments  plays  a  critical  role 
in  the  global  organic  carbon  cycle.  Quantitatively,  there  is  more  organic  carbon  (OC) 
derived  from  non-living  organic  molecules  occurring  in  recent  marine  sediments,  soil 
humus  and  dissolved  in  seawater  than  in  all  living  land  plants  and  marine  organisms 
combined  (Hedges  et  al.,  2000).  Approximately  90%  of  OC  buried  globally  in  the  ocean 
is  sequestered  in  shelf  and  slope  sediments  (Hedges  and  Keil,  1995  and  references 
therein),  and  removal  of  organic  carbon  from  the  biosphere  and  storage  in  sediments 
affects  many  key  processes  such  as  global  carbon  budgets  and  long  term  controls  on 
atmospheric  CO2  and  O2. 

To  understand  the  sources,  transport  and  degradation  of  OC  in  the  oceans,  the 
nature  and  distribution  of  OM  in  marine  sediments  must  be  studied.  In  addition  to  the 
intrinsic  reactivity  of  organic  molecules,  interactions  between  individual  organic 
compounds  and  macromolecular  organic  matter  (MOM)  are  important  as  they  have  the 
potential  to  influence  the  overall  fate  of  specific  organic  compounds,  including 
contaminants  (Farrington  and  Quinn,  1971;  Lee,  1977;  Zegouagh  et  al.,  1996).  For 
example,  compounds  which  are  “free”  and  non-covalently  associated  with  MOM  may  be 
more  easily  degraded  than  those  that  are  “bound”  and  covalently  linked  by  an  ester,  ether 
or  C-C  bond.  Compounds  defined  as  “free”  are  those  that  are  solvent  extractable,  whereas 
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those  released  only  by  harsher  chemical  treatments,  are  considered  “non-extractable”  or 
“bound”.  It  is  important  to  note  that  not  all  “non-extractable”  compounds  are  covalently 
bound  to  OM  and  some  may  be  physically  entrapped  in  OM  or  mineral  matrices  (Knicker 
and  Hatcher,  1997  and  references  therein)  and  only  released  when  the  matrices 
themselves  are  altered  by  chemical  treatments.  Conversely,  the  extractability  of  a 
compound  does  not  necessarily  imply  that  it  is  susceptible  to  degradation.  For  example, 
organic  compounds  sorbed  to  mineral  surfaces  may  be  readily  extractable,  but  physically 
protected  from  degradation  (Keil  and  Hedges,  1993).  Nevertheless  these  different 
associations  may  exert  strong  controls  on  the  preservation  of  labile  organic  compounds 
(Wakeham,  1999),  as  well  as  their  residence  times  in  different  carbon  reservoirs.  Binding 
to  OM,  for  example,  may  slow  the  delivery  of  a  terrestrial  derived  organic  compound  to 
the  ocean  through  association  with  relatively  immobile  particulate  phases  (e.g.  soils). 

Complex  mixing  of  autochthonous  and  allochthonous  materials  occur  on 
continental  margins.  To  assess  the  relative  contributions  of  marine,  terrestrial  and  fossil 
OC  to  marine  sediments,  abundances  and  distributions  of  molecular  markers,  such  as 
hydrocarbons,  n-alkanoic  acids  and  alkanols  are  typically  examined.  Specific 
characteristics  of  these  markers  are  studied  and  include  carbon  preference  and  carbon 
isotopic  compositions.  For  example,  short  chain  (n-Cn  to  M-C22)  H-alkanoic  acids  and 
alkanols  with  an  even-over-odd  predominance  are  considered  to  be  derived 
predominantly  from  marine  algae,  bacteria  or  zooplankton  whereas  their  longer  chain  (n- 
C24  to  /1-C34)  counterparts  have  terrestrial  origins  from  higher  plant  epicuticular  waxes. 
Alkanoic  acids  derived  from  sedimentary  microbial  biomass  may  also  contribute 
significantly  (Hollander  and  Smith,  2001).  For  w-alkanes,  an  odd-over-even 
predominance  in  carbon  chain  length  >n- C23  is  typical  of  higher  plants  and  is  numerically 
represented  by  the  carbon  preference  index  (CPI,  =  Z[odd  numbered  alkane  abundances 
from  M-C23  to  n-C35]/Z[even  numbered  alkane  abundances  from  n-C 22  to  H-C34]).  Initial 
plant  wax  CPI  values  are  high  (CPI  =  4-40;  Collister  et  al.,  1994),  but  tend  towards  1 
with  increasing  degradation,  as  well  as  increasing  thermal  maturity  for  sources  such  as 
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petroleum,  petroleum  source  rocks  and  continental  rocks  of  intermediate  thermal  maturity 
(Brassell  and  Eglinton,  1980). 

Stable  carbon  isotope  (813C)  compositions  of  OM  can  be  used  as  indicators  of 
both  source  and  carbon  flow.  Sedimentary  bulk  OM  is  comprised  of  contributions  from 
marine,  terrestrial  and  bacterial  sources,  and  this  is  reflected  in  the  isotope  ratios 
obtained.  Terrestrial  sources  derived  from  C3  vegetation  typically  exhibit  813C  values  ~  - 
26%o  whereas  marine  sources  are  more  enriched  (~  -21%o;  Fry  and  Sherr,  1984;  Prahl  et 
al.,  1994).  Natural  abundance  radiocarbon  (A14C)  measurements  are  also  informative  as 
14C  serves  as  a  sensitive  tracer  of  OM  inputs.  All  recently  photosynthesized  carbon 
contains  significant  amounts  of  14C,  which  is  produced  in  the  Earth’s  atmosphere  by 
cosmic  ray  interactions,  and  makes  its  way  into  the  biosphere  mainly  via  photosynthetic 
fixation  of  CO2.  Additionally,  any  natural  product  derived  from  carbon  reservoirs  in 
communication  with  the  atmosphere  since  1952  contains  an  excess  quantity  of  14C  due  to 
incorporation  of  14C  produced  by  the  above-ground  testing  of  thermonuclear  weapons. 
The  5730  year  half  life  of  14C  allows  compounds  of  geologic  age  to  be  easily 
distinguished  as  they  contain  no  14C.  Continental  margin  sediments  contain  significant 
quantities  of  both  modem  and  “pre-aged”  terrestrial  organic  carbon  (Eglinton  et  al.,  1997; 
Masiello  and  DrufFel  2001;  Blair  et  al.,  2003;  Goni  et  al.,  2005;  Drenzek  et  al.,  submitted) 
with  potential  sources  of  “pre-aged”  carbon  arising  from  remobilization  of  soil  carbon, 
weathering  of  sedimentary  rock  and  resuspension  and  redistribution  of  marine  organic 
matter  (Mollenhauer  et  al.,  2005). 

In  this  study,  A14C  and  813C  measurements  of  bulk  sediments  and  residues 
resulting  from  sequential  chemical  treatments  are  combined  with  analyses  of  the 
molecular  composition  of  extracts  released  by  these  treatments.  By  examining  the  nature 
and  distribution  of  OM  in  marine  sediments,  we  aim  to  gain  insight  not  only  into  the 
varied  sources  of  OM,  but  also  by  understanding  their  associations  we  hope  to  make 
predictions  concerning  their  overall  fate. 
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4.2  MATERIALS  AND  METHODS 


4.2.1  Study  Area  and  Samples 

Geographic  locations,  dates,  depths  and  details  of  the  sites  examined  in  this  study 
are  detailed  in  Table  4.1  and  Figure  4.1.  These  samples  represent  a  variety  of  coastal  and 
continental  margin  settings,  many  of  which  have  previously  been  subjected  to  molecular 
isotopic  studies  (see  references,  Table  4.1).  The  samples  are  divided  into  five  general 
groups  determined  by  their  depositional  settings.  Group  ‘a’  includes  organic-rich  marine 
sediments  from  Guaymas  basin  (GB),  a  hydrothermal  system  with  active  petroleum 
production  in  the  sediment;  sediments  from  the  upwelling  systems  of  the  Peru  Margin 
(PM);  and  surface  (0-1  cm)  and  pre-industrial  sediments  (13-15  cm)  from  the  Benguela 
upwelling  (BU)  system  off  of  the  coast  of  Namibia,  Africa.  Santa  Monica  Basin  (SMB) 
sediments  were  examined  from  the  periphery  of  the  basin  which  has  oxygenated  bottom 
waters  (Station  1),  and  from  the  less  oxygenated  depocenter  (Station  2)  from  surface  (0-1 
cm),  pre-bomb  (10-11  cm)  and  pre-industrial  (30-32  cm)  sediments.  Group  ‘b’  is 
comprised  of  open  shelf  and  river  dominated  sediments  from  the  New  England  Margin 
(NEM);  the  Washington  Margin  (WM),  including  a  site  from  the  inner  shelf  (Station  1) 
and  surface  (0-1  cm)  and  pre-industrial  (34-36  cm)  sediments  from  further  offshore 
(Station  4),  as  well  as  surface  (0-2  cm)  and  pre-industrial  (100-101  cm)  sediments  from 
the  Eel  River  Margin  (ERM)  which  has  a  narrower  flood  plain  than  WM.  Group  ‘c’ 
encompasses  coastal  sites  where  sediments  were  recovered  from  shallow  (<10  cm)  or 
intertidal  zones.  These  include  anoxic  sediments  from  the  Pettaquamscutt  River  (PR), 
sediments  from  Penobscot  Bay  (PB);  sediments  from  Lowes  Cove  (LC),  which  is  an 
intertidal  mudflat;  marine  sediments  from  Blakely  Harbor  (BH);  and  salt  marsh 
sediments  from  West  Falmouth  (WF).  Group  ‘d’  includes  river-dominated  sediments 
from  the  Mackenzie  shelf  (MS)  and  Beaufort  slope  (BS)  in  the  Canadian  Arctic.  Group 
‘e’  is  comprised  of  surface  (0-2  cm)  sediments  from  the  Ross  Sea  (RS),  Antarctica  and 
one  deeper  sediment  horizon  (12-15  cm).  Significant  inputs  of  pre-aged  OC  to  sediments 
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Figure  4.1  Location  of  sediments  examined  in  this  study. 


Table  4.1  Locations,  dates,  depths  and  sedimentary  environments  of  samples  analyzed  in  this  study. 
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c  MC  -  multi-core. 
rfGGC  -  giant  gravity  core. 
e  FC  -  freeze-core. 

^GS  -  grab  sampler. 


from  groups  ‘d’  and  ‘e’  have  previously  been  observed  (Ohkouchi  et  al.,  2003;  Goni  et 
al.,  2005;  Drenzek  et  al.,  submitted). 

All  groups  include  various  end-members  in  terms  of  depositional  environments,  and 
include  sediments  deposited  under  both  oxic  and  anoxic  bottom  waters  with  variable 
contributions  from  marine,  terrestrial,  ancient  and  petroleum  carbon  sources.  This  enables 
an  assessment  of  the  potential  affects  of  diagenesis  and  different  oxidation  states  on  the 
nature  and  distribution  of  OM  in  marine  sediments  and  enables  predictions  to  be  made 
about  more  complex  sediments  that  are  a  mixture  of  the  end-member  environments 
examined  here. 

4.2.2  Sequential  Lipid  Extraction 

Sediments  were  air-dried  and  homogenized  with  a  mortar  and  pestle  before 
passing  through  a  1mm  sieve.  Depending  on  the  OC  content  and  quantity  of  sediment 
available,  two  to  twenty  grams  of  each  sediment  sample  were  taken  and  sequentially 
treated  as  previously  described  (White  et  al.,  2005).  Briefly,  the  sediments  were  extracted 
with  dichloromethane  (DCM)  and  methanol  (MeOH)  (9:1)  by  pressurized  fluid  extraction 
(100°C,  1000  psi).  The  resulting  total  lipid  extracts  (TLE)  were  reduced  in  volume  by 
rotary  evaporation  and  percolated  through  a  copper  column  plugged  with  combusted 
glass  wool  to  remove  sulfur  and  any  entrained  particulate  matter  prior  to  further  analysis. 
The  solvent  extraction  step  was  not  performed  in  this  study  for  GB,  MS  and  BS  samples 
because  the  sediments  had  already  been  solvent  extracted  using  a  Soxhlet  apparatus  as 
part  of  previous  investigations  (Pearson  et  al.,  2005;  Drenzek  et  al.,  submitted).  Aliquots 
of  the  TLE  from  these  sediments  were  provided  for  GC/MS  analysis.  Solvent  extracted 
sediment  from  all  samples  was  refluxed  at  70°C  for  2  h  with  a  solution  of  0.5  N 
potassium  hydroxide  (KOH)  (in  MeOH,  100  mL)  and  water  (H2O)  (20  mL).  After 
cooling,  the  reaction  mixture  was  separated  by  centrifugation  (1500  rpm),  the  supernatant 
removed,  and  the  remaining  sediment  rinsed  with  MeOH,  followed  by  DCM,  and  then 
hexane.  These  extracts  were  combined  with  the  supernatant  in  a  separatory  funnel.  The 
alkaline  extract  (supernatant)  was  back-extracted  with  hexane  (3  x  50  mL)  to  obtain  a 
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neutral  fraction  before  acidification  with  4N  hydrochloric  acid  (HC1)  to  pH  2  and 
extraction  with  DCM  (3  x  50  mL)  to  obtain  an  acidic  fraction.  Both  fractions  were 
reduced  in  volume  by  rotary  evaporation  and  dried  by  passing  through  a  small  glass 
pipette  column  of  anhydrous  Na2S04  plugged  with  glass  wool  (this  also  served  to  remove 
any  particles  from  the  extracts). 

4.2.3  Bulk  Sediment  Analyses 

Aliquots  of  dried  unextracted  (TOC),  extracted  (EX-RES)  and  saponified 
sediment  residues  (SA-RES)  were  treated  with  10%  HC1  to  remove  any  carbonates.  The 
sediments  and  residues  were  then  analyzed  for  organic  carbon  content  with  a  Carlo  Erba 
1108  elemental  analyzer.  The  stable  carbon  isotope  ratio  (SI3C)  and  radiocarbon 
abundance  (A14C)  of  the  sediments  were  measured  on  purified  carbon  dioxide  (CO2)  after 
their  combustion  in  the  presence  of  cupric  oxide  (CuO).  Stable  carbon  values  were 
determined  by  isotope  ratio  mass  spectrometry  (irMS),  and  14C  content  by  accelerator 
mass  spectrometry  (AMS)  at  the  National  Ocean  Sciences  Accelerator  Mass 
Spectrometry  (NOSAMS)  facility  at  Woods  Hole  Oceanographic  Institution  (WHOI), 
Woods  Hole,  MA,  after  conversion  of  the  CO2  to  graphite  (McNichol  et  al.,  1994).  All 
14C  measurements  are  normalized  to  8*3C  values  of  -25%o  and  expressed  as  A14C  values, 
which  is  the  per  mille  (%o)  deviation  from  the  international  standard  for  14C  dating, 
Standard  Reference  Material  4990B  "Oxalic  Acid’’.  Precision  for  813C  and  A14C 
measurements  are  ~0.1  and  2-5%o  (see  Appendix  for  errors  specific  to  each  sample 
analyzed),  respectively.  For  our  samples,  the  reported  A14C  =  \fme  (195°'*)X-1]1000  (Stuiver 
and  Polach,  1977),  where  X= 1/8267  (year '),  fm-  fraction  modem  14C  (corrected  for 
isotopic  fractionation  using  813C),  and  “x”  equals  the  year  of  collection  (see  Table  4.1). 
This  corrects  for  decay  of  14C  since  time  of  collection  to  time  of  measurement.  Isotopic 
and  elemental  analysis  of  the  TOC,  EX-RES  and  SA-RES  from  WF,  and  the  TOC  from 
the  ERM  (100-101  cm  sample  only),  GB,  BU,  RS,  MS  and  BS  had  previously  been 
performed  (see  Table  4.2),  and  were  not  repeated  for  this  study. 
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4.2.4  Gas  Chromatography-Mass  Spectrometry  (GC/MS)  Analysis 

Aliquots  of  TLEs  and  saponified  neutral  and  acidic  extracts  were  spiked  with  an 
internal  standard  ( n-C alkane,  4  pg),  prior  to  derivatization  with  bis-(trimethylsilyl)- 
trifluoroacetamide  (BSTFA)  and  pyridine  to  improve  chromatography.  For  compound 
identification,  electron  ionization  (El)  spectra  were  acquired  on  an  Agilent  6890  series 
gas  chromatograph  interfaced  to  an  Agilent  5973  mass  selective  detector  (MSD).  A  post 
column  split  to  a  flame  ionization  detector  allowed  for  quantification  of  compounds  of 
interest  relative  to  the  internal  standard.  Compounds  were  separated  on  a  J&W  DB-5MS 
capillary  column  (60  m  x  0.32mm  ED,  0.25  pm  film  thickness)  with  helium  carrier  gas  at 
a  constant  flow  of  1  ml  min'1.  The  initial  GC  oven  temperature  was  40°C  (1  minute  hold) 
and  was  ramped  at  20°C  min'1  until  130°C  and  then  at  2°C  min'1  until  320°C  (30  minute 
hold).  Spectra  were  acquired  between  m/z  40-650  at  a  scan  rate  of  1  cycle  s'1.  Compounds 
present  in  the  extracts  were  identified  from  mass  spectra  and  gas  chromatographic 
retention  characteristics.  The  flame  ionization  detector  response  was  calibrated  regularly 
by  injection  of  n-alkane,  n-alkanoic  and  n-alkanol  standards. 

4.3  RESULTS  AND  DISCUSSION 

4.3.1  General  Trends  in  Surface  Sediments 

For  the  majority  of  surface  sediments  examined  in  this  study,  the  A14C  of  the  bulk 
untreated  sediment  prior  to  solvent  extraction  (TOC)  have  A14C  values  that  range  from  ~ 
-300  to  -50%o  (Table  4.2  and  Figure  4.2  a,  b  and  c),  reflecting  a  mixture  of  modem,  pre¬ 
aged  and  fossil  OM  of  marine  and  terrestrial  origin.  These  A14C  values  compare  well  with 
those  obtained  for  other  core-top  sediments  previously  reported  in  the  literature  (e.g., 
McNichol  et  al.,  1993;  Wang  et  al.,  1996;  Eglinton  et  al.,  1997;  Masiello  and  DrufFel, 
2003;  Hwang  et  al.,  2005;  Komada  et  al.,  2005).  Salt  marsh  sediments  from  WF  are, 
however,  more  enriched  (+44%o,  Figure  4.2c),  likely  due  to  incorporation  of  recently 
synthesized  vascular  plant  biomass  from  the  roots  of  the  marsh  grass  Spartina 
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Table  4.2  Percent  organic  carbon,  stable  carbon  and  radiocarbon  abundances0  for  bulk  samples  and  their  respective  residues. 
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a  All  A14C  values  displayed  in  table  are  corrected  for  14C  decay  since  time  of  collection  to 
time  of  measurement  (see  text). 

b  TOC  is  the  untreated  sediment  prior  to  solvent  extraction. 
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Figure  4.2  A14C  and  813C  of  unextracted  sediment  residues  (TOC)  for  (a)  organic-rich 
marine  sediments;  (b)  open  shelf  and  river  dominated  sites;  (c)  near-coast  sediments;  (d) 
Arctic  sediments;  (e)  Antarctic  sediments  and  (f)  all  core-top  samples  compared.  Error 
bars  for  A14C  and  813C  are  10%o  and  0.1%o  respectively  and  are  smaller  than  the  symbols 
shown  in  the  figure  and  thus  are  not  visible. 
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altemiflora  into  the  sedimentary  matrix.  Other  exceptions  include  sediments  from  ERM, 
RS,  MS  and  BS  which  are  all  significantly  depleted  (A14C  ~  -750  to  -300%o;  Figure  4.2b, 
d  and  e)  due  to  varying  contributions  of  fossil  organic  carbon.  This  fossil  carbon  is 
derived  from  kerogen  from  bedrock  in  the  case  of  ERM  sediments  (Blair  et  al.,  2004); 
pre-aged  soil  material  and  bitumen  or  kerogen  in  MS  and  BS  sediments  (Goni  et  al., 
2005;  Drenzek  et  al.,  2005)  and  from  relict  organic  carbon  from  the  Antarctic  continent  in 
the  case  of  RS  sediments  (Ohkouchi  et  al.,  2003). 

The  813C  values  of  the  bulk  TOC  for  all  sites  range  from  -28  to  -18%o  (Table  4.2, 
Figure  4.2),  but  within  each  group  of  sediments,  it  does  not  vary  by  more  than  6%o.  In 
general,  the  organic-rich  marine  sediments  (Figure  4.2a)  are  the  most  C-enriched  (-23.4 
to  -19%o),  and  the  Arctic  and  Antarctic  sediments  (Figure  4.2d  and  e)  are  the  most 
depleted  (-27.4  to  -23.8%o),  likely  reflecting  significant  contributions  of  enriched  marine 
carbon  in  the  former  and  ,3C-depleted  marine  carbon  and/or  fossil  carbon  in  the  latter. 
Sediments  from  near-coast  and  open  shelf  and  river  dominated  sites  lie  between  these 
two  extremes  with  813C  of  the  bulk  TOC  ranging  from  -25.4  to  -20.8%o  (Figure  4.2b  and 
c).  The  WF  sample  is  the  exception.  It’s  more  enriched  (-18.2%o)  8,3C  value  consistent 
with  those  of  other  salt  marsh  systems  where  the  C4  plant  Spartina  altemiflora  is  the 
dominant  macrophyte  (Middleburg  et  al.,  1997).  A  weak  positive  correlation  between 
isotopic  compositions  (813C  vs.  A14C)  of  the  TOC  of  all  sediments  examined  is  apparent 
(Figure  4.2e;  r^O.5,  p=0.0004)  indicating  that  A14C  as  well  as  813C  can  be  used  to 
indicate  the  source  of  organic  matter.  It  is  worth  noting  that  this  correlation  improves 
(^=0.6)  if  the  sediments  with  significant  fossil  contributions  (MS,  BS  and  RS  Fairy  and 
Emperor)  are  considered  outliers. 

In  general,  there  is  little  difference  between  the  A14C  value  of  the  TOC  and  the 
corresponding  solvent  extracted  sediment  residue  (EX-RES)  for  the  majority  of  sites 
examined  (Table  4.2).  This  concurs  with  previous  work  on  sediments  from  NEM  that 
demonstrated  the  A14C  of  individually  isolated,  extractable  alkenones,  which  are 
biomarkers  derived  from  marine  algae,  are  similar  to  that  of  bulk  OM  (Ohkouchi,  Pers. 
Comm.).  GC/MS  analysis  of  TLE  (Figure  4.3  i)  from  each  sample  provides  insights  into 
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the  composition  and  source  of  materials  removed  by  solvent  extraction.  It  should  be 
borne  in  mind,  however,  that  much  of  the  TLE  is  not  amenable  to  GC,  so  interpretations 
should  be  made  with  caution.  In  each  case,  varying  proportions  of  short  (n-Cn  to  n-C 22) 
and  long  (n- C24  to  M-C34)  chain  alkanoic  acids  and  alkanols  with-even-over-odd 
predominance  (Figure  4.4-4.8),  reflect  inputs  from  marine  and  terrestrial  sources, 
respectively.  There  are  also  significant  contributions  from  n-alkanes,  which  are 
predominantly  odd  carbon  numbered.  In  samples  where  they  contribute  <5%  to  the  GC- 
amenable  portion  the  even  n-alkanes  concentrations  are  below  the  detection  limit. 
Generally,  CPI  values  range  from  1.6  (predominantly  fossil-derived)  to  13.3 
(predominantly  plant-wax  derived).  There  are  also  significant  contributions  from  sterols 
for  which  the  assignment  of  sources  is  less  straightforward  (Volkman  2005),  but  is 
dominated  by  those  with  carbon  numbers  from  C27-C29. 

At  the  majority  of  sites,  systematic  offsets  in  A14C  between  the  EX-RES  and  SA- 
RES  are  observed  (Table  4.2).  With  the  exception  of  WM  St.  1  and  GB  (discussed  later) 
the  SA-RES  is  always  more  depleted  than  that  of  the  EX-RES  indicating  removal  of  more 
14C  enriched  material  by  the  process  of  saponification.  GC/MS  analysis  indicates  that  the 
GC-amenable  fraction  of  the  saponification  extracts  (SAE)  is  mainly  composed  of  short- 
chain  alkanoic  acids  (n-Cn  to  H-C22)  (Figure  4.3  ii )  with  strong  even-over-odd 
predominance  (Figure  4.9-4.13),  indicating  a  marine  algal  or  bacterial  source.  This  is 
supported  in  a  study  by  Ohkouchi  et  al.  (2003)  that  demonstrated  relatively  uniform  A14C 
values  for  ester-bound  short  chain  alkanoic  acids  in  surface  sediments  from  the  Ross  Sea 
that  were  more  enriched  relative  to  the  bulk  sediment  and  reflected  the  post-bomb 
dissolved  inorganic  carbon  (DIC)  reservoir.  Such  observations  provide  evidence  for  the 
protection  of  labile  marine  carbon  by  chemical  binding  with  the  sedimentary  matrix.  Free 
tt-alkanoic  acids  are  especially  susceptible  to  degradation  if  they  are  short-chain  and 
unsaturated  (Sun  and  Wakeham,  1994;  Zegouagh  et  al.,  1996)  and  although  they  are  also 
derived  from  bacterial  communities  responsible  for  this  degradation,  these  would  most 
likely  be  present  in  the  TLE  (as  intact  phospholipids)  as  opposed  to  the  SAE.  Long  chain 
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Figure  4.3  Relative  contribution  of  n-alkanoic  acids,  n-alcohols  and  alkanes  for  a) 
organic-rich  marine  sediments;  b)  open  shelf  and  river  dominated  sediments;  c)  near¬ 
coast  sediments;  d)  Arctic  sediments  and  e)  Antarctic  sediments  to  /)  the  total  lipid 
extract  (TLE)  and  (ii)  the  saponification  extract  (SAE).  See  text  for  description  of ‘other’. 
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Figure  4.4  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids,  H-alkanols  and 
n-alkanes  in  the  total  lipid  extract  (TLE)  in  organic-rich  marine  sediments.  represents  the 
unsaturated  counterpart  of  the  previous  n-alkanoic  acid.  Values  superimposed  on  histogram  bars 
represent  those  that  are  off  scale  and  odd  over  even  predominance  (OEP)  for  n-alkanes  is 
indicated.  Diagonal  lines  are  present  where  compounds  are  below  detection  limits. 
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Figure  4.6  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids,  n 
alkanols  and  n-alkanes  in  the  total  lipid  extract  (TLE)  to  Near-coast  sediments. 
represents  the  unsaturated  counterpart  of  the  previous  n-alkanoic  acid.  Values 
superimposed  on  histogram  bars  represent  those  that  are  off  scale  and  odd  over  even 
predominance  (OEP)  for  n-alkanes  is  indicated. 
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Figure  4.7  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids,  n- 
alkanols  and  n-alkanes  in  the  total  lipid  extract  (TLE)  to  Arctic  sediments.  **’  represents 
the  unsaturated  counterpart  of  the  previous  n-alkanoic  acid.  Values  superimposed  on 
histogram  bars  represent  those  that  are  off  scale  and  odd  over  even  predominance  (OEP) 
for  n-alkanes  is  indicated. 
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Figure  4.8  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids  and  n- 
alkanols  in  the  total  lipid  extract  (TLE)  to  Antarctic  sediments.  **’  represents  the 
unsaturated  counterpart  of  the  previous  n-alkanoic  acid.  Values  superimposed  on 
histogram  bars  represent  those  that  are  off  scale.  Diagonal  lines  are  present  where 
compounds  are  below  detection  limits. 
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Figure  4.10  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids  and 
n-alkanols  in  the  saponification  extract  (SAE)  for  open  shelf  and  river  dominated 
sediments.  **’  represents  the  unsaturated  counterpart  of  the  previous  n-alkanoic  acid. 
Values  superimposed  on  histogram  bars  represent  those  that  are  off  scale. 
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Figure  4.11  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids  and 
n-alkanols  in  the  saponification  extract  (SAE)  for  near-coast  sediments.  **’  represents  the 
unsaturated  counterpart  of  the  previous  n-alkanoic  acid. 
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Figure  4.12  Percent  contribution  to  compound  class  for  individual  n-alkanoic  acids  and 
n-alkanols  in  the  saponification  extract  (SAE)  for  Arctic  sediments.  **’  represents  the 
unsaturated  counterpart  of  the  previous  n-alkanoic  acid.  Values  superimposed  on 
histogram  bars  represent  those  that  are  off  scale.  Diagonal  lines  are  present  where 
compounds  are  below  detection  limits. 
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M-alkanoic  acids  originating  from  higher  plants  are,  however,  more  resistant  to 
degradation  most  likely  due  to  a  protective  matrix  association  (Zegouagh  et  al.,  1996). 

To  determine  the  significance  of  the  contributions  of  different  biomarkers  observed  in  the 
TLE  and  SAE,  the  proportion  of  TOC  that  is  solvent  extractable  and  saponifiable  as  well 
as  the  quantity  of  this  material  that  is  GC  amenable  was  calculated  (Table  4.3).  To 
calculate  the  percentage  of  TOC  that  was  extractable,  portions  of  the  extracts  were 
weighed  and  converted  to  carbon  equivalents  considering  contributions  from  hydrogen 
and  oxygen  (an  average  of  80%  of  biogenic  compounds  examined  here  are  carbon 
derived).  This  value  was  then  divided  by  the  amount  of  organic  carbon  in  the  bulk 
unextracted  sediment  (TOC,  Table  4.2).  The  fact  that  only  GC-amenable  compounds  are 
discussed  is  a  limitation  of  this  study,  but,  this  highlights  the  strengths  of  the  isotopic 
approach  as  differences  in  isotopic  compositions  between  sediments  and  residues 
provides  information  about  the  carbon  pools  that  are  not  GC-amenable.  Examination  of 
A14C  shifts  between  bulk  residues  and  the  A14C  of  the  extracts  removed  that  lead  to  these 
shifts  has  been  performed  in  one  case  (WF)  and  these  values  are  consistent  with  one 
another  (White  et  al.,  2005). 

In  the  context  of  these  overall  aforementioned  trends,  the  isotopic  abundances  of 
sediments  and  their  residues  as  well  as  their  compound  distributions  in  “free”  (TLE)  and 
“bound”  (SAE)  pools,  will  be  discussed  in  more  detail.  The  difference  in  isotope 
abundances  (A14C  and  813C)  of  the  residues  (EX-RES  represented  by  grey  filled  circles; 
SA-RES  represented  by  open  circles)  from  the  bulk  TOC  (black  filled  circle  at  the  origin 
of  the  plot)  are  shown  for  core-top  sediments  in  figures  4.13-4.17.  These  will  also  be 
discussed  for  sediment  samples  below  the  mixed  layer  (Figure  4.18).  Where  present, 
arrows  indicate  significant  offsets  between  the  residues  as  a  result  of  the  chemical 
treatments. 
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Table  4.3  Percent  organic  carbon  for  all  samples  that  is  extractable,  saponifiable  and  non-hydrolyzable  and  the  percent  of  the 
extracts  that  is  amenable  by  gas  chromatography. 
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4.3.2  Organic-Rich  Marine  Sediments 

Isotopic  abundances  of  the  TOC  for  organic-rich  marine  sediments  from  surface 
sediments,  display  relatively  small  offsets  between  bulk  sediments  and  their  respective 
residues  as  a  result  of  the  chemical  treatments.  This  is  demonstrated  by  the  clustering  of 
the  EX-RES  and  SA-RES  in  both  A,4C  and  813C  space  around  the  origin,  which 
represents  the  isotopic  composition  of  the  TOC  (Figure  4.13).  This  indicates  that  the 
organic  carbon  removed  by  the  sequential  treatments  is  isotopically  similar  to  that  in  the 
original  bulk  sediment  as  well  as  that  remaining  in  the  residues.  Analysis  of  organic 
phosphorus  in  sediments  of  the  Peru  Margin  by  sequential  treatment  also  indicated  little 
variance  in  organic  phosphorus  between  different  pools  of  OM  suggesting  that  the 
different  pools  of  OM  have  similar  compositions  (Laarkamp,  2000).  This  is  likely  due  to 
the  fact  that  these  organic-rich  sediments  are  predominantly  derived  from  a  marine 
source.  Although  the  isotopic  offsets  between  the  residues  are  small,  some  general 
patterns  are  evident.  Excluding  the  GB  sample,  the  mean  offset  in  AA14C  and  A813C 
between  the  EX-RES  and  TOC  is  7  and  0.2%o  respectively  indicating  a  slight  enrichment 
in  EX-RES  in  both  isotopes  relative  to  the  TOC.  This  is  consistent  with  previous  work  by 
Eglinton  et  al.  (1997)  which  observed  a  I4C-depleted  TLE  relative  to  the  TOC  in  surface 
sediments  from  the  Black  Sea.  Other  work  by  Hwang  et  al.  (2005),  consistently 
demonstrates  lipid  that  are  depleted  in  A14C  and  813  relative  to  bulk  sediment.  Excluding 
the  GB  sample  again,  the  mean  offset  in  AA14C  and  A813C  between  the  SA-RES  and  TOC 
can  be  calculated  and  is  -9  and  0.3%o  respectively.  The  SA-RES  is  enriched  in  813C 
similarly  to  the  EX-RES,  but  is  slightly  depleted  in  A14C  relative  to  the  bulk  TOC, 
indicating  the  removal  of  14C-enriched  material  by  saponification. 

Comparisons  between  sediments  from  the  SMB  depocenter  (Station  2)  where 
bottom  water  oxygen  contents  are  <0.1  ml  l'1  (Hickey,  1992)  and  the  periphery  (Station 
1)  where  bottom  water  oxygen  contents  are  higher,  have  previously  been  made  to 
examine  the  differential  degradation  of  alkenones  and  differences  in  fatty  acid 
composition  (Gong  and  Hollander,  1997,  1999;  Mollenhauer  and  Eglinton,  in 
preparation).  Alkenone  contents  were  observed  to  be  lower  in  sediments  from  the 


106 


periphery  than  those  from  the  depocenter,  but  at  depths  corresponding  to  deposition 
under  anoxic  conditions  in  the  depocenter,  the  alkenone  contents  were  of  equal 
abundance,  indicating  the  control  of  post-depositional  oxidative  degradation  on  these 
compounds.  In  this  study,  surface  sediment  OC  contents  of  the  periphery  were  lower 
(2.2%,  Table  4.2)  compared  to  the  depocenter  supporting  these  previous  findings  (3.2%, 
Table  4.2).  Several  experiments  examining  surface  water  primary  production,  physical 


Figure  4.13  Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  symbols) 
and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AAl4C  %o  for  organic-rich  marine  sediments  (core-top 
samples  only).  Error  bars  represent  0.1%o  error  for  813C  values  and  10%o  error  for  A14C. 
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and  chemical  properties  of  the  upper  water  column  and  sedimentation  processes  have 
been  performed  to  demonstrate  that  the  primary  difference  between  the  sediments  of  the 
depocenter  and  the  periphery  was  the  oxygen  content  (Gong  and  Hollander,  1999).  Very 
similar  distributions  of  extractable  and  saponifiable  biomarkers  are  observed  in  the 
surface  sediments  from  the  two  stations  (Figures  4.3a,  4.4  and  4.9),  with  the  only 
significant  difference  being  the  absence  of  short  chain  alkanols  (n-C  12-C22)  in  the 
periphery  sediments  (Station  1)  suggesting  enhanced  degradation  of  these  compounds  in 
the  presence  of  oxygen.  Overall,  the  surface  bulk  sediment  and  corresponding  residues 
from  the  depocenter  are  more  enriched  (-85  to  -66%o)  than  those  from  the  periphery  (-146 
to  -129%o;  Table  4.2),  indicating  that  14C  rich  organic  material  that  is  younger  and  likely 
more  labile,  is  preferentially  degraded  at  the  periphery. 

Sediments  from  GB  behave  very  differently  than  the  other  organic-rich  marine 
sediments  (Figure  4.13)  and  the  residues  become  successively  more  enriched  in  14C  as  the 
sequential  treatments  are  performed.  The  TOC  is  the  most  depleted  (-197%o)  followed  by 
the  EX-RES  (-184%o)  and  the  SA-RES  (-163%o).  Guaymas  Basin  is  an  unusual 
environment  in  which  hydrothermal  petroleum  is  produced  in  the  sediments  as  a  result  of 
magmatic  heating.  The  successive  enrichment  of  residues  observed  results  from  the 
removal  of  relatively  !4C  depleted  material  at  each  step.  Solvent-extractable  n-alkanoic 
acids  isolated  from  the  microbial  mat  of  Beggiotoa,  exhibit  depleted  A14C  values  (-418%o 
to  -227%o)  that  are  similar  to  that  of  hydrothermal  petroleum,  indicating  the  consumption 
of  pre-aged  carbon  by  the  entire  bacterial  assemblage  (Pearson  et  al.,  2005).  For  this 
study,  we  have  examined  the  sediment  underlying  the  microbial  mat,  and  our  results  are 
consistent  with  Pearson  et  al.  suggesting  that  solvent  extraction  removes  compounds  that 
are  relatively  depleted  in  14C.  Notably,  the  TLE  in  this  study  is  comprised  predominantly 
of  short  chain  «-alkanoic  acids  (-80%,  Figure  4.3a).  The  SAE  reveals  a  significant 
(-70%)  contribution  of  polycyclic  aromatic  hydrocarbons  (PAHs)  including 
naphthalenes,  phenanthrenes  and  their  alkylated  derivatives,  fluoranthene,  anthracene  and 
fluorene  (shown  as  “other”  in  Figure  4.3a).  These  PAHs  are  apparently  nonextractable, 
but  are  unlikely  to  be  covalently  bonded  to  the  OM  as  they  are  not  functionalized. 
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Instead,  they  may  be  encapsulated  in  the  sedimentary  matrix,  possibly  via  association 
with  organic  matter  or  mineral  grains.  These  compounds  are  directly  derived  from  the 
14C-depleted  hydrothermal  petroleum,  thus  explaining  the  remaining  enriched  residue. 

4.3.3  Open  Shelf  and  River-Dominated  Sediments 

Significant  offsets  in  A14C  are  observed  between  TOC  and  EX-RES  of  surface 
sediments  from  WM  St.  1  and  ERM  and  to  a  lesser  extent  NEM  (Figure  4.14).  These 
offsets  are,  however,  in  opposite  directions.  In  the  case  of  WM  St.l,  the  A14C  EX-RES  is 
significantly  more  depleted  than  that  of  the  TOC.  For  this  latter  sample,  the  offset  may 
reflect  the  weak  association  and  facile  removal  of  more  recently  synthesized  material  by 
solvent  extraction,  leaving  behind  a  more  depleted  residue.  This  interpretation  is 
supported  by  the  molecular-level  composition  of  the  TLE  (Figure  4.3b),  which  has 
significant  contributions  from  long-chain  («-C24  to  /1-C34)  alkanols  (~20%)  and  alkanoic 
acids  (~5%)  with  an  even  over  odd  predominance  as  well  as  long  chain  (/1-C23  to  H-C31) 
alkanes  (~15%)  with  an  odd-over-even  predominance  (CPI  =  10.4;  Figure  4.5).  These 
biomarker  compounds  are  characteristic  of  vascular  plant  biomass  and  are  particularly 
abundant  at  WM  St.  1,  which  is  on  the  inner  shelf.  Absolute  abundances  of  n-alkanes  on 
the  WM  continental  shelf  average  165  (ig  g'1  carbon  (C)  (Prahl  et  al.,  1994).  Radiocarbon 
analysis  of  individual  terrestrial  biomarkers  for  a  site  proximal  to  St.  1  determined  that 
they  are  younger  than  the  TOC  (Eglinton,  Pers.  Comm.)  supporting  these  bulk  level 
observations.  No  significant  offset  in  A14C  between  TOC  and  EX-RES  (Table  4.2)  is 
observed  at  WM  St.  4,  which  is  further  offshore  than  St.  1  (upper  continental  slope)  for 
surface  sediments,  or  those  at  a  depth  of  34-36  cm  corresponding  to  ~318  yr  bp  (M. 
Uchida,  Pers.  Comm.).  This  is  likely  due  to  a  decrease  in  the  contribution  of  terrestrial 
plant  material,  which  is  not  particularly  evident  in  the  composition  of  the  TLE  (Figure 
4.3),  but  previous  research  has  indicated  decreasing  concentrations  in  this  material  with 
increasing  distance  offshore,  for  example,  absolute  abundances  of  «-alkanes  on  the  slope 
are  91  |ig  g'1  C  and  52  |ig  g'1  C  in  basin  sediments  (Prahl  et  al.,  1994). 
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In  sharp  contrast  to  WM1,  the  radiocarbon  content  of  TOC  from  ERM  0-2cm  is 
considerably  more  depleted  (~100%o)  than  the  EX-RES  (Figure  4.14).  Analysis  of  the  n- 
alkanes,  which  comprise  ~20%  of  the  GC-resolvable  TLE  (Figure  4.3b)  reveals  a  low 
CPI  (1.6;  Figure  4.5),  betraying  the  presence  of  petrogenic  hydrocarbons  associated  with 
the  supply  of  OM  from  erosion  of  ancient  sedimentary  rocks  on  the  adjacent  continent 
(Blair  et  al.,  2004).  Large  quantities  of  unaltered  bedrock  containing  predominantly  a 


and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AA14C  %o  for  open  shelf  and  river  dominated  sediments 
(core-top  samples  only).  Error  bars  represent  0.  l%o  error  for  813C  values  and  10%o  error 
for  A14C. 
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type  III  kerogen  and  plant  debris-rich  surface  soils  are  delivered  from  the  continent  as  a 
result  of  large  scale  flood  events.  As  a  consequence,  widespread  and  rapid  (0.2-0.4  cm 
yr'1)  accumulation  of  modem  sediment  occurs  along  the  shelf  break  (Sommerfield  and 
Nittrouer,  1999).  The  sediments  examined  in  this  study  are  from  the  middle  of  the  flood 
depocenter  on  the  shelf. 


Figure  4.15  Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  symbols) 
and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AA14C  %o  for  near-coast  sediments.  Error  bars  represent 
0.1  %o  error  for  813C  values  and  10%o  error  for  A14C. 
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4.3.4  Near-Coast  Sediments 

Near-coast  sediments  (Figure  4.15)  display  similar  isotopic  shifts  to  the  general 
trends  observed  except  for  LC,  which  exhibits  a  more  I4C  enriched  TOC  relative  to  the 
EX-RES  (-82  and  -107%o  respectively;  Table  4.2).  This  is  likely  due  to  the  removal  of 
14C-enriched  terrestrial  biomass  as  indicated  by  the  abundance  of  long  chain  n-alkanoic 
acids  and  n-alkanols  with  even/odd  dominance  as  well  as  n-alkanes  with  a  CPI  of  3.7 
(Figures  4.3c  and  4.6).  A  significant  shift  in  813C  and  A14C  space  is  also  observed  for  BH 
between  the  EX-RES  and  SA-RES  (Figure  4.15),  which  becomes  more  depleted  in  both 
813C  and  A14C  due  to  the  removal  of  relatively  13C  and  14C  enriched  material  that  is 
likely  marine  derived  as  indicated  by  short  chain  n-alkanoic  acids  and  alkanols  that 
dominate  the  SAE  (Figure  4.3c). 

4.3.5  Arctic  Sediments 

Surface  sediments  and  residues  from  MS  and  BS  exhibit  the  most  depleted  A14C 
values  (-823  to  -618%o)  observed  in  this  study  (Table  4.2).  Similar  shifts  in  A14C  are 
observed  in  all  sediments  (Figure  4.16),  with  the  TOC  being  the  most  enriched  in  14C  and 
the  SA-RES  the  most  depleted  resulting  from  the  removal  of  OC  that  is  increasingly 
enriched  in  14C  as  the  sequential  treatments  progress.  This  offset  is  greatest  (~  -100%o) 
between  the  EX-RES  and  SA-RES  (Figure  4.16)  likely  due  to  the  SAE  being  composed 
predominantly  (-60-100%)  of  fresh  marine  biomarkers  such  as  short  chain  n-alkanoic 
acids  (as  in  Zegouagh  et  al.,  1996)  and  n-alkanols  (Figure  4.3d  ii)).  The  TLEs,  are 
dominated  by  short  chain  n-alkanoic  acids,  which  have  previously  been  shown  to  have 
modem  14C  values  (Drenzek  et  al.,  submitted)  with  lesser  contributions  from  long  chain 
n-alkanoic  acids  and  short  and  long  chain  n-alkanols.  There  are  also  varying 
contributions  from  n-alkanes  (5-25%)  with  OEP  ranging  from  1. 9-2.7  previously  shown 
to  be  derived  from  a  mixture  of  higher  plant  leaf  waxes  and  erosion  of  sedimentary  rocks 
in  the  drainage  basin  (Yunker  et  al.,  1993;  Drenzek  et  al.,  submitted).  What  remains  as 
the  non-hydrolyzable  portion  of  the  TOC  is  extremely  depleted  in  14C  (-823  to  -721%o). 
These  values  are  consistent  with  A14C  compositions  of  pyrolysate  alkanes  from  Station  5 
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(as  in  Drenzek  et  al.,  submitted)  that  are  thought  to  be  derived  from  kerogen  and/or 
vascular  plant  cutan  that  is  pre-aged  on  the  continent  prior  to  delivery  to  the  Beaufort 
Sea. 


Figure  4.16  Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  symbols) 
and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AA14C  %o  for  Arctic  sediments.  Error  bars  represent 
0.1  %o  error  for  813C  values  and  10%o  error  for  A14C. 
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4.3.6  Antarctic  Sediments 

Significant  offsets  in  A14C  and  813C  are  observed  between  sediments  and 
corresponding  residues  in  Antarctic  sediments  from  the  Ross  Sea  (Figure  4.17).  These 
sites  are  the  same  as  those  examined  by  Ohkouchi  et  al.  (2003).  The  Emperor  and  Gentoo 
sites  are  in  regions  with  extensive  blooms  of  haptophyte  algae;  Chinstrap  represents  a 
region  where  diatoms  bloom  during  the  austral  summer.  In  contrast,  large  algal  blooms 
have  not  been  observed  at  the  Fairy  site  which  is  the  most  proximal  to  the  Ross  ice  shelf 


Figure  4.17  Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  symbols) 
and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AA14C  %o  for  Antarctic  sediments  (core-top  samples 
only).  Error  bars  represent  0.1%o  error  for  813C  values  and  10%o  error  for  A14C. 
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(Dunbar  et  al.,  1985).  For  all  sites  except  Emperor,  the  TOC  is  similar  to,  or  more 
enriched  in  14C  than  the  EX-RES  and  significant  variation  in  S13C  is  not  observed 
between  fractions.  Where  observed,  this  difference  between  TOC  and  EX-RES  represents 
the  extraction  of  OC  that  is  relatively  enriched  in  14C  as  reflected  in  TLEs  that  are 
dominated  by  short  chain  n-alkanoic  acids  (30-60%,  Figure  4.3e  i)).  The  latter  have  been 
shown  to  be  significantly  enriched  in  14C  relative  to  the  TOC,  exhibiting  AI4C  values 
corresponding  to  that  of  post-bomb  surface  ocean  dissolved  inorganic  carbon  in  this 
region  (Ohkouchi  et  al.,  2003).  The  Emperor  site,  however  displays  a  significantly  more 
depleted  A,4C  of  TOC  (-557%o)  relative  to  EX-RES  (-413%0;  Table  4.2),  implying 
removal  of  14C-depleted  OC.  Analysis  of  the  TLE  (Figure  4.3e,  where  other  represents 
the  contribution  of  sterols)  does  not  reveal  significant  differences  in  biomarker 
contributions  in  Emperor  compared  to  the  other  sites,  but  the  fraction  of  OC  that  is 
extractable  is  less  (17%)  than  that  observed  in  the  other  Ross  Sea  surface  sediments  (28- 
51%;  Table  4.3).  Long  chain  n-alkanoic  acids  at  the  Emperor  and  Fairy  sites  are 
relatively  depleted  in  14C  (e.g.  /2-C28  FA  =  -375%o;  Ohkouchi,  et  al.,  2003)  compared  to 
their  short  chain  counterparts,  and  these  depleted  values  may  have  a  larger  impact  on  the 
isotopic  signature  at  the  Emperor  site  where  the  %TOC  extractable  is  less.  Fairy 
sediments  and  residues  exhibit  the  most  depleted  A14C  values  (-736  to  -713%o)  of  all  of 
the  Ross  Sea  sediments  examined  here,  and  differences  in  contributions  of  biomarkers  to 
the  TLE  and  SAE  are  observed  at  the  Fairy  site  compared  to  the  other  Ross  Sea 
sediments.  The  contribution  of  sterols  to  the  TLE  is  much  less  than  at  other  sites  (~10% 
compared  to  26-46%;  Figure  4.3e)  and  the  remainder  of  the  GC-amenable  components  of 
the  TLE  are  comprised  of  equal  contributions  of  short  and  long  chain  n-alkanoic  acids 
whereas  at  other  sites,  short  chain  n-alkanoic  acids  contribute  twice  as  much  as  their  long 
chain  counterparts.  This  lesser  contribution  of  sterols  and  short  chain  n-alkanoic  acids  at 
the  Fairy  site  may  result  from  either  a  deficit  in  the  delivery  of  fresh  marine  organic 
carbon  from  the  photic  zone  (due  to  the  absence  of  large  algal  blooms  at  the  Fairy  site), 
preferential  degradation  of  sedimentary  biomarkers,  dilution  of  the  signal  by  OM  from  a 
different  source  with  a  different  composition  from  lateral  transport  (Villinski  and  Hayes, 
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2005)  or  fossil  OC  inputs  (Ohkouchi  and  Eglinton,  2005).  Either  way,  these  inputs  result 
in  an  overall  decrease  in  the  proportion  of  easily  extractable,  fresh  14C  enriched  material 
in  the  sediments,  and  thus  give  rise  to  only  a  small  offset  in  A14C  TOC  and  EX-RES. 

Differences  in  A14C  and  813C  for  EX-RES  and  SA-RES  from  Fairy,  Gentoo  and 
Chinstrap  indicated  by  similar  vectors  observed  (symbolized  by  grey  dashed  arrows 
between  grey  filled  diamonds  and  open  diamonds)  in  Figure  4.17.  The  SA-RES  are 
significantly  more  depleted  resulting  from  the  removal  of  relatively  14C  enriched 
biomarkers  dominated  by  short  chain  «-alkanoic  acids  and  sterols  (Figure  4.3  e  ii).  The 
Fairy  site  is  the  anomaly  and  this  offset  is  much  smaller  and  in  the  opposite  direction. 
The  SAE  indicates  less  of  a  contribution  of  short  chain  n-alkanoic  acids  (-20%) 
compared  to  the  other  Ross  Sea  surface  sediments  (-60-80%;  Figure  4.3  e  i)),  which  may 
cause  this  difference. 

4.3.7  Diagenetic  Considerations 

The  effects  of  diagenesis  on  the  nature  and  distribution  of  OM  in  marine 
sediments  was  examined  for  select  samples.  In  general,  the  significant  offsets  in  ,4C  and 
13C  between  bulk  sediments  and  residues  observed  in  surface  sediments  diminish  at  depth 
(e.g.  ERM  and  WM4;  Figure  4.18).  In  some  samples,  however,  small  isotopic  offsets 
observed  in  the  surface  sediments  are  preserved  with  depth  (SMB  St.  2  and  Chinstrap). 
Changes  in  the  relative  contributions  of  biomarkers  with  depth  are  also  observed  in  some 
samples  and  help  explain  the  isotopic  offsets  observed. 

For  organic-rich  marine  sediments  from  BU,  offsets  in  A14C  values  between  bulk 
sediment  and  residues  observed  in  the  surface  sediment  (Figure  4.13),  are  not  observed  at 
depth  (13-15  cm;  Figure  4.18)  and  instead  the  EX-RES  and  SA-RES  are  similarly 
depleted  in  14C  and  13C  from  the  TOC.  Analysis  of  the  TLE  indicates  smaller  relative 
contributions  of  short  chain  n-alkanoic  acids  at  13-15  cm,  and  a  corresponding  increase 
of  long  chain  w-alkanoic  acids  (Figure  4.3a)  reflecting  the  increased  lability  of  the  former. 
The  quantity  of  TOC  that  is  extractable  and  saponifiable  is  also  less  at  depth  (Table  4.3), 
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and  thus  removal  of  these  pools  of  organic  carbon  (TLE  and  SAE)  has  less  of  an  impact 
on  the  isotopic  composition  of  the  bulk  organic  carbon  pool. 

Pre-bomb  sediments  deposited  in  the  depocenter  (10-1  lcm  interval,  185  yr  bp; 
Mollenhauer  and  Eglinton,  2004)  of  the  SMB  display  a  similar  offset  in  A14C  between 
EX-RES  and  SA-RES  to  that  observed  in  the  surface  sediment  (Figure  4.13  and  4.18), 
again  most  likely  due  to  the  absence  of  easily  extractable  younger  material  (Table  4.2). 


Figure  4.18  Differences  in  813C  of  the  extracted  residue  (EX-RES;  filled  grey  symbols) 
and  saponified  residue  (SA-RES;  open  symbols)  from  813C  of  bulk  sediment  (TOC;  filled 
black  symbol)  expressed  as  A813C  %o  vs.  the  differences  in  A14C  of  the  EX-RES  and  SA- 
RES  from  the  TOC  expressed  as  AA14C  %o  for  sediment  horizons  below  the  surface 
mixed  layer.  Error  bars  represent  0.1%o  error  for  813C  values  and  10%o  error  for  A14C. 
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Short  chain  n-alkanoic  acids  contribute  less  to  the  TLE  and  SAE  (~20%)  at  10-11  cm 
compared  to  0-1  cm  (~60%  and  ~70%  respectively)  and  instead  there  is  a  greater 
contribution  of  alkanes  to  the  TLE  (~20%)  and  sterols  to  the  SAE  (~70%;  Figure  4.3a  i)). 
This  decrease  in  contribution  of  the  short  chain  n-alkanoic  acids,  indicates  that  there  may 
be  some  post  depositional  degradation  occurring,  although  this  is  not  conclusive  without 
comparison  to  an  equivalent  periphery  sediment  horizon.  It  is  also  possible  that  changes 
in  source  function  may  give  rise  to  these  differences. 

Sediments  from  the  depocenter  deposited  under  oxic  conditions  prior  to  the 
switch  to  dysoxia  were  also  examined  (30-32  cm  interval,  corresponding  to  545  yr  bp; 
Mollenhauer  and  Eglinton,  2004)  and  show  significantly  lower  contributions  of  short 
chain  n-alkanoic  acids  (~3%)  and  higher  proportions  of  plant  wax  n-alkanes  (~70%; 
Figure  4.3a),  the  latter  having  a  CPI  of  2.9  (Figure  4.4).  This  decrease  in  n-alkanoic  acids 
may  be  due  to  either  diagenesis  and/or  the  higher  bottom  water  concentrations  of  oxygen 
when  these  sediments  were  deposited  that  would  result  in  enhanced  degradation  of  these 
compounds.  Compared  to  the  periphery  sediments,  however,  the  contribution  of  short 
chain  n-alkanoic  acids  is  still  very  low  suggesting  that  this  difference  has  to  be  due  to 
significantly  higher  bottom  water  oxygen  at  the  depocenter  relative  to  the  periphery.  The 
isotopic  offset  between  the  EX-RES  and  TOC  at  30-32  is  similar  to  that  observed  in  the 
surface  sediment  (Figure  4.13  and  4.18)  and  the  EX-RES  is  more  14C-enriched  relative  to 
the  TOC,  although  no  accompanying  enrichment  in  813C  is  observed.  There  is  also  no 
significant  offset  in  either  813C  or  A14C  between  SA-RES  and  TOC  at  this  depth. 

There  are  subtle  shifts  in  biomarker  contributions  in  the  extracts  of  WM4  0-1  cm 
compared  to  34-36  cm  sediments  (Figure  4.3b)  and  the  contribution  of  short  chain  n- 
alkanoic  acids  halves  from  ~40  to  ~20%  with  depth  whereas  long  chain  n-alkanoic  acids 
and  n-alkanols  contribute  more  significantly.  Sterol  contributions  also  decrease  with 
depth  from  ~35  to  ~15%  in  the  TLE,  but  increase  in  the  bound  pool  from  40  to  60%. 
Significant  offsets  in  isotopic  composition  are  not  observed  between  TOC  and  EX-RES, 
but  the  SA-RES  is  relatively  more  14C-depleted  relative  to  the  TOC  similar  to  general 
observations  for  the  surface  sediment  (Figure  4.14).  The  813C  is  invariant  ( — 22%o) 
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between  WM4  sediment  residues  and  with  depth  (Figure  4.18),  and  thus  we  can  infer  that 
diagenesis  does  not  appear  to  have  a  significant  effect  on  the  isotopic  composition  of  the 
sediments  at  this  site.  It  is  of  interest  to  note  that  WM4  is  positioned  at  the  depocenter  of 
the  slope  and  the  main  degradation  may  have  taken  place  during  lateral  transport  to  this 
site. 

The  large  offset  between  TOC  and  EX-RES  observed  in  the  surface  sediment  of 
ERM  is  not  observed  at  depths  of  100-101  cm  (Figure  4.18),  corresponding  to  circa  1830 
a.d.  (N.  Drenzek,  Pers.  Comm.).  The  distribution  of  n-alkanes  at  100-101  cm  suggests  a 
larger  terrestrial  and  smaller  contribution  of  petrogenic  hydrocarbons  relative  to  that 
contributing  to  the  surface  sediment,  and  this  is  reflected  in  the  CPI  (2.1;  Figure  4.5).  A 
relative  decrease  in  the  contribution  of  short  chain  n-alkanoic  acids  is  also  observed  with 
depth  (2%  compared  to  35%  in  the  surface  sediment),  but  significant  differences  are  not 
observed  for  the  SAE  with  depth  (Figure  4.3b  ii).  Less  material  is  solvent  extractable  at 
100-101  cm  compared  to  the  surface  sediment  (~5  and  50%  respectively;  Table  4.3)  and 
this,  coupled  with  differences  in  biomarker  contributions,  may  reflect  a  change  in  the 
type  of  OM  supplied  from  the  Eel  River  as  opposed  to  a  diagenetic  trend.  Increased 
sediment  accumulation  rates,  finer  sediments  and  increased  preservation  of  flood-layers 
have  been  observed  for  the  Eel  shelf  since  AD  1950  and  attributed  to  climate,  specifically 
increases  in  flood  magnitude  and  frequency  (Sommerfield  et  al.,  2002).  More  recently,  it 
has  been  suggested  that  these  differences  are  due  to  wet  conditions  coupled  with 
increased  sediment  supply  as  a  result  of  timber  harvesting  (Leithold  and  Blair,  2001). 
Considering  all  major  OC  sources  to  the  ERM  (terrestrial,  kerogen  and  marine),  only  the 
terrestrial  source  has  been  found  to  consistently  change  over  time  and  the  pre- 1950s 
contribution  of  terrestrial  material  relative  to  bedrock  was  found  to  be  greater  (Leithold 
and  Blair,  2001).  This  supports  observations  from  the  n-alkane  data  in  this  study.  It 
therefore  appears  likely  that  in  this  instance  the  absence  of  isotopic  offset  between  TOC 
and  EX-RES  at  100-101  cm  may  be  due  to  a  shift  in  OM  type. 

For  Antarctic  sediments  from  the  Chinstrap  site,  differences  in  §13C  and  A14C 
between  sediments  and  residues  for  sediment  from  12-15  cm  (Figure  4.18)  are  similar  to 
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those  observed  for  0-2  cm  (Figure  4.17).  Analysis  of  the  extracts  released  by  the 
sequential  treatments  indicates  a  significant  increase  in  the  relative  contribution  of  sterols 
at  12-15  cm  to  the  SAE  (Figure  4.3).  This  suggests  that  sediments  deposited  at  this  depth 
likely  have  a  more  significant  input  from  fresh  marine  organic  matter  than  those 
deposited  at  0-2  cm,  which  are  likely  more  diluted  by  allochthonous  inputs. 

4.5  CONCLUSIONS 

Differences  in  A14C  and  molecular  composition  of  different  OM  fractions  in 
sediments  are  evident  at  a  variety  of  sites  with  different  depositional  settings.  Only  a 
small  fraction  of  OM  present  in  sediments  is  GC-amenable,  and  thus  differences  in 
isotopic  compositions  between  bulk  sediments  and  residues  provides  information 
regarding  the  composition  of  pools  of  OC  that  are  undetectable  by  conventional 
geochemical  techniques. 

The  findings  of  this  study  compare  well  with  previous  work  that  has  examined  the 
A14C  of  different  OM  fractions.  Analysis  of  detrital  aggregates,  sediment  floe  and 
sediments,  revealed  14C-depleted  total  lipids  pools  relative  to  individual  amino  acid  and 
carbohydrate  fractions  (Wang  and  Druffel,  1996;  Wang  et  al.,  1998).  This  is  similar  to 
observations  in  organic-rich  marine  sediments  where  in  general,  the  EX-RES  is  more 
enriched  relative  to  the  TOC  as  a  result  of  the  removal  of  labile,  14C-depleted  lipids.  A 
comparative  study  between  sediments  from  the  Southern  Ocean  and  NE  Pacific  indicates 
that  a  ‘younger’  more  labile  organic  matter  pool  consisting  of  amino  acids  and 
carbohydrates  was  being  preserved  with  depth  at  the  NE  Pacific  site  where  sedimentation 
rates  and  TOC  contents  were  higher  suggesting  protection  of  these  more  labile 
compounds  by  sorption  onto  particle  surfaces  (Wang  and  Druffel,  2001).  This  is  similar 
to  the  findings  of  this  study  that  provides  evidence  for  the  protection  of  14C-enriched 
labile  marine  organic  carbon  (predominantly  short  chain  H-alkanoic  acids)  that  are 
associated  with  insoluble  sediment  residues. 
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It  is  important  to  reiterate  that  systematic  variations  in  isotopic  compositions  are 
not  only  observed  within  specific  depositional  settings,  such  as  organic-rich  marine 
sediments,  but  also  between  these  groupings.  The  two  most  obvious  examples  of  this  are 
the  similarities  observed  between  Arctic  and  Antarctic  samples  (excluding  Emperor)  and 
between  ERM  (0-1  cm),  Guaymas  basin  and  the  Emperor  site.  The  former  group  is 
characterized  by  the  presence  of  non-hydrolyzable  fossil  carbon  that  gives  rise  to 
extremely  depleted  A14C  values  for  the  bulk  unextracted  sediment.  These  A14C  values 
become  successively  more  depleted  as  relatively  14C-  enriched  material  is  removed  by  the 
sequential  chemical  treatments  at  each  step.  The  latter  group  are  also  controlled  by  the 
presence  of  fossil  carbon,  but  in  these  environments  it  is  solvent-extractable  and  the 
depleted  A14C  of  the  bulk  unextracted  sediment  becomes  more  enriched  as  this  fossil 
carbon  is  removed  by  the  sequential  chemical  treatments.  Comparisons  such  as  this 
across  different  depositional  settings  are  extremely  valuable  and  increase  our 
understanding  of  the  nature  and  distribution  of  organic  matter  in  marine  sediments. 

To  summarize  this  broad  array  of  data,  it  is  important  to  consider  the  global 
significance  of  this  work.  The  majority  of  OC  burial  in  the  contemporary  marine 
environment  occurs  on  open  shelf  and  slope  margins  (-45%)  as  well  as  deltaic  sediments 
(-44%;  Hedges  and  Keil  1995).  No  significant  offset  in  A14C  between  the  TOC  and  the 
EX-RES  is  observed  for  the  majority  of  surface  sediments  from  these  environments 
examined  in  this  study  (n=7).  The  EX-RES  from  other  surface  sediments  from  these 
environments  (n=6)  can,  however,  be  relatively  14C-depleted  compared  to  the  TOC  due  to 
the  removal  of  relatively  ,4C-enriched  material  that  has  been  more  recently  synthesized 
and  may  be  terrestrial  in  origin,  as  in  WM  sediments.  These  observations  in  A14C  offsets 
between  TOC  and  EX-RES  are  by  far  the  most  common  observed  in  this  study,  and  can 
be  expected  to  apply  for  -90%  of  organic  carbon  buried  in  contemporary  marine 
sediments.  The  main  exception  arises  from  inputs  of  fossil  carbon  and  in  sediments  from 
ERM,  GB  and  RS  Emperor,  the  EX-RES  is  more  14C-enriched  relative  to  the  TOC  due  to 
the  removal  of  fossil-derived  carbon  that  contains  no  14C.  Organic-rich  marine  sediments 
either  display  no  offset  in  ,4C  between  the  TOC  and  EX-RES  (n=2)  or  the  EX-RES  is 
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slightly  enriched  relative  to  the  TOC  (n=2)  due  to  the  removal  of  material  slightly 
depleted  in  14C.  These  highly  productive  environments  are,  however,  not  as  significant  as 
shelf,  slope  and  deltaic  systems  and  are  only  responsible  for  ~5%  of  global  organic 
carbon  burial  (Hedges  and  Keil,  1995).  For  the  majority  (n=13)  of  surface  sediments 
examined  in  this  study,  independent  of  depositional  setting,  the  SA-RES  is  more  14C- 
depleted  compared  to  the  EX-RES.  This  indicates  that  the  protection  of  relatively  ,4C- 
enriched  material  of  predominantly  marine  algal  origin,  is  a  widespread  phenomenon.  It 
is  clear  that  14C  analysis  provides  a  powerful  tool  for  determining  both  the  source  and 
fate  of  marine,  terrestrial  and  fossil  organic  carbon  in  the  marine  environment  as  well  as 
enhancing  our  understanding  of  the  general  distribution  of  organic  carbon  in 
contemporary  marine  sediments. 
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Santa  Monica  Basin,  St  2  (0-1  cm)  TOC  1  mg  OS-47278  -22.3  0.9277  0.0032  -78.4 

Santa  Monica  Basin,  St.  2  (0-1  cm)  EX-RES  1  mg  OS-47470  -22.3  0.9394  0.0044  -66.8 

Santa  Monica  Basin,  St.  2  (0-1  cm)  SA-RES  1  mg  OS-47321  -22.0  0.9204  0.0048  -85.6 
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Eel  River  Margin,  St  1(0-2  cm)  TOC  1  mg  OS-45557  -25.2  0.6847  0.0027  -319.7 


Eel  River  Margin,  St.  1  (0-2  cm)  EX-RES  1  mg  OS-45510  -25.6  0.7884  0.0032  -216.7 

Eel  River  Margin,  St.  1  (0-2  cm)  SA-RES  1  mg  OS-45511  -26.1  0.7737  0.0028  -231.4 

Eel  River  Margin,  St.  1(100-101  cm)  EX-RES  1  mg  OS-47478  -25.0  0.5376  0.0029  -466 
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0  Ai4C  values  displayed  in  table  are  raw  data  and  are  not  corrected  for  in-situ  14C  decay, b  EX-RES  is  the  extracted  sediment 
residue, c  S  A-RES  is  the  saponified  sediment  residue, d  TOC  is  the  untreated  sediment  prior  to  solvent  extraction. 
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CHAPTER  5 


RADIOCARBON  AS  AN  INVERSE  TRACER  OF  FOSSIL  FUEL-DERIVED 
CONTAMINANTS  IN  SEDIMENTS 


ABSTRACT 


We  have  exploited  natural  abundance  variations  in  radiocarbon  (14C)  as  an  inverse 
tracer  of  fossil  fuel-derived  contaminants  (14C  free)  within  different  chemically  defined 
fractions  of  contemporary  organic  materials  (OM)  (modem  14C  content)  from  Select 
marine  and  freshwater  sediments.  Hydrophobic  organic  contaminants  (HOCs)  are  known 
to  interact  with  OM,  but  these  processes  are  difficult  to  follow  and  quantitatively 
constrain.  Using  natural  abundance  14C  measurements  it  is  possible  to  make  quantitative 
assessments  through  the  principle  of  isotope  mass  balance.  This  study  focuses  on 
sediments  impacted  by  various  fossil  fuel-derived  HOCs  including  petroleum,  creosote, 
coal  tar,  polychlorinated  biphenyls  (PCBs),  benzotriazoles  (BZTs),  2,2- 
bis(chlorophenyl)-l,l,l-trichlorethane  (DDT)  and  other  chlorinated  pesticides.  A  mineral 
oil  contaminated  soil  and  an  urban  aerosol  sample  have  also  been  examined.  The  samples 
were  sequentially  chemically  treated,  first  by  solvent  extraction  and  then  saponification  to 
examine  the  quantity  of  HOC  and/or  its  derivatives  present  in  both  the  extractable  and 
unextractable  fractions.  Radiocarbon  analysis  of  the  original  bulk  sample  and  residues 
resulting  from  the  sequential  treatments  was  performed  to  trace  HOCs  and  their 
metabolites  in  the  different  pools  of  OM.  The  resulting  extracts  were  examined  using  gas 
chromatography/mass  spectrometry  (GC/MS)  to  identify  and  quantity  HOCs  and/or 
metabolites  released  from  the  sequential  treatments.  For  most  samples,  the  majority  of 
fossil  fuel-derived  carbon  is  solvent  extractable,  but  at  a  few  sites  either  chemical  or 
isotopic  evidence  indicates  varying  quantities  of  HOCs  and/or  derivatives  present  within 
the  insoluble  residues.  The  quantities  of  these  nonextractable  HOCs  are  for  the  most  part, 
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not  significant,  and  thus  incorporation  of  HOCs  into  insoluble  sediment  residues  cannot 
be  considered  an  important  sink  in  aquatic  sediments.  This  suggests  that  the  ability  of 
HOCs  to  persist  in  the  environment  is  likely  controlled  by  physical  associations  as 
opposed  to  chemical  linkages.  In  this  form,  HOCs  may  be  more  available  to  biota  and 
more  likely  to  be  subjected  to  transport  processes,  especially  if  sediment  conditions 
change  e.g.  if  the  sediment  is  resuspended  and  remobilized.  Overall,  the  lack  of 
significant  incorporation  of  HOCs  into  nonextractable  residues  suggests  that  these 
compounds  have  a  greater  potential  to  exert  longer-term  adverse  affects  on  the 
sedimentary  environment  in  which  they  persist. 

5.1  INTRODUCTION 

Elucidating  the  nature  of  interactions  between  hydrophobic  organic  contaminants 
(HOCs)  and  organic  matter  (OM)  is  of  great  importance  for  determining  the  fate  of 
xenobiotics  in  the  environment.  Much  of  the  work  to  date  has  focused  on  soils 
contaminated  with  polycyclic  aromatic  hydrocarbons  (PAHs)  (Richnow  et  al.,  1994, 
1995,  1997,  1998,  1999  and  2000)  and  other  HOCs  such  as  pesticides  (Gevao  et  al.,  2000 
and  references  therein),  fumigants  (Xu  et  al.,  2003),  the  fungicide  cyprodinil  (Dec  et  al., 
1997),  2,4,6-trinitrotoluene  (TNT)  (Eriksson  et  al.,  2004),  and  polychlorinated  biphenyls 
(PCBs)  (Kohl  and  Rice,  1998).  Experimental  and  analytical  techniques  for  determining 
the  presence  of  bound  HOCs  and  their  residues  vary  (Northcott  and  Jones,  2000  and 
references  therein),  but  in  most  instances,  specific  interactions  are  determined  by 
sequential  chemical  treatments.  Typically,  samples  are  first  solvent  extracted  to  isolate 
compounds  with  non-covalent  associations  to  OM,  such  as  hydrophobic  sorption,  charge 
transfer  complexes  and  hydrogen  bonding.  The  HOCs  and  residues  released  from  this 
process  are  referred  to  as  “free”  or  “unbound”  and  those  that  remain  associated  with  the 
OM  are  defined  as  “nonextractable”  or  “bound”.  Subsequent  chemical  treatments 
performed  on  the  remaining  extracted  soil  residues  or  humic  fractions  target  linkages  by 
specific  types  of  covalent  bond  between  HOCs  and  OM.  Alkaline  or  acid  hydrolysis 
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methods  are  specifically  aimed  at  cleaving  ester  bonds  (Richnow  et  al.,  1994,  1997, 1998 
and  2000),  and  harsher  degradation  techniques  such  as  boron  trichloride  degradation  and 
hydrogenolysis  with  rhodium-on-charcoal  target  ether  as  well  as  ester  linkages  (Richnow 
et  al.,  1994).  Analytical  pyrolysis  and  thermal  desorption  methods  have  also  been 
employed  to  release  compounds  associated  by  carbon-carbon  bonds  as  well  as  other 
bonds  (Richnow  et  al.,  1995;  Northcott  and  Jones  2000,  and  references  therein).  It  is 
important  to  note  that  not  all  nonextractable  residues  are  covalently  associated  to  OM  and 
some  may  be  physically  entrapped  in  OM  or  mineral  matrices  (Gevao  et  al.,  2000  and 
references  therein)  and  only  released  when  the  matrices  themselves  are  altered  by 
chemical  treatments. 

For  covalent  binding  to  occur,  the  contaminant  or  its  metabolic  intermediate  must 
possess  a  functional  group,  such  as  hydroxyl,  carboxyl,  nitro,  amino  or  phosphate  group 
or  other  reaction  center.  The  functional  group  of  a  contaminant  will  largely  determine  its 
reactivity  towards  OM,  and  it  has  been  demonstrated  in  soils  that  phenolic  and  aniline 
compounds  have  the  greatest  tendency  to  form  bound  residues  (Roberts,  1984).  Simple 
hydrocarbons,  however,  are  relatively  inert  and  undergo  the  same  fate  only  after  a 
preliminary  attack  and  oxidation  by  microorganisms  (Northcott  and  Jones,  2000). 

The  formation  of  nonextractable  residues  in  soils  can  be  quantitatively  significant, 
and  as  much  as  70%  of  anthracene  and/or  its  metabolites  have  been  demonstrated  to  bind 
to  soil  OM  in  aerobic  microcosm  experiments  (Richnow  et  al.,  1998).  This  quantity, 
however,  is  highly  variable  and  in  a  separate  incubation  of  phenanthrene,  only  ~6%  of  the 
corresponding  metabolites  became  bound  (Richnow  et  al.,  2000).  This  difference  is  likely 
due  to  variations  in  soil  microflora  and  the  efficiency  of  microorganisms  to  metabolize 
different  PAHs  into  more  reactive  metabolites  that  can  subsequently  bind  to  OM.  Bound 
residue  formation,  however,  has  been  proposed  as  a  major  sink  for  HOCs  in  soils  and  has 
been  considered  as  a  bioremediation  strategy  (Dec  and  Bollag,  1997).  Interactions 
between  HOCs  and  OM  have  not  been  examined  as  extensively  in  aquatic  sediments,  but 
initial  studies  suggest  that  binding  is  a  potentially  important  sink  (Beller  and  Simoneit, 
1988;  Richnow  et  al.,  1994;  Reddy  et  al.,  2000a;  Schwarzbauer  et  al.,  2003)  Aquatic 
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sediments,  however,  provide  a  very  different  chemical  and  microbial  environment  to  that 
found  in  soils  (White  et  al.,  2005)  and  it  is  therefore  of  interest  to  determine  whether  the 
significance  of  binding  in  soils  can  be  extrapolated  to  sediments. 

Natural  abundance  variations  in  radiocarbon  content  (AI4C)  have  previously  been 
used  to  trace  fossil  fuel  or  petrochemical-derived  HOCs  into  various  pools  of  organic 
matter  (Reddy  et  al.,  2002a;  White  et  al.,  2005)  as  well  as  enabling  quantitative 
assessments  of  fossil  carbon  contributions  through  the  principle  of  isotope  mass  balance 
(Pearson  and  Eglinton,  2000;  Petsch  et  al.,  2001).  The  absence  of  14C  in  fossil-fuel- 
derived  contaminants  serves  as  an  inverse  tracer  intrinsic  to  the  contaminant,  thereby 
negating  the  need  to  use  a  labeled  derivative  that  may  not  accurately  reflect  in  situ 
processes  (Steinberg  et  al.,  2003).  The  sensitivity  of  this  approach  is  an  important 
consideration  and  in  general,  with  a  conservative  estimate  of  analytical  error  of  10%o,  at 
least  1%  of  organic  carbon  in  a  particular  organic  matter  pool  of  interest  must  be  fossil 
fuel  derived  for  it  to  significantly  influence  the  overall  A14C  (White  et  al.,  2005).  The 
inverse  tracer  approach  has  been  previously  applied  to  anaerobic  salt  marsh  sediments  to 
determine  the  fate  of  petroleum  residues  persisting  in  sediments  for  over  30  years  (White 
et  al.,  2005).  At  this  site,  significant  binding  to  OM  was  not  observed,  and  this  was 
attributed  to  the  relatively  inert  nature  of  the  petroleum  residues  and  the  lack  of  more 
reactive  metabolites  resulting  from  the  absence  of  active  biodegradation  at  the  site  (Slater 
et  al.,  2005).  Understanding  whether  more  reactive  HOCs,  or  HOCs  in  more  reactive 
sediments  (e.g.  aerobic  as  opposed  to  anaerobic),  would  give  rise  to  more  significant 
binding  of  HOCs  to  OM  is  therefore  of  interest.  To  explore  this,  we  have  collected 
surface  sediments  from  a  variety  of  locations  contaminated  with  petroleum,  creosote,  coal 
tar,  polychlorinated  biphenyls  (PCBs),  benzotriazoles  (BZTs),  2,2-bis(chlorophenyl)- 
1,1,1-trichloroethane  (DDT)  and  other  chlorinated  pesticides. 

In  petroleum  contaminated  sediments,  PAHs  such  as  naphthalene  and 
phenanthrene  and  their  alkylated  derivatives  are  of  interest  due  to  their  toxicity  to 
organisms,  as  well  as  the  potential  of  their  corresponding  metabolic  intermediates  to  bind 
to  OM  (Richnow  et  al.,  1994, 1995, 1997,  1998, 1999,  2000).  The  structure  of  the  parent 
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Figure  5.1  Structures  of  hydrophobic  organic  contaminants  examined  in  this  study. 
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PAH  naphthalene  (Figure  5.1a  i))  and  its  metabolites  identified  from  aerobic  (Figure  5.1a 
ii;  Heitkamp  et  al.,  1997)  and  anaerobic  degradation  (Figure  5.1a  iii;  Zhang  and  Young 
1997)  in  sediments  are  shown  as  an  example.  Creosote  and  coal  tar  are  composed  of 
approximately  85%  PAHs,  10%  phenolic  and  5%  N-,  S-,  and  O-heterocyclic  compounds, 
and  examples  of  these  compounds  are  also  presented  (Figure  5.1b).  Polychlorinated 
biphenyls  are  mixtures  of  up  to  209  individual  chlorinated  compounds  known  as 
congeners  with  the  general  structure  Ci2Hio-xClx  where  x  =  1-10  (Figure  5.1c).  Structures 
of  the  most  abundant  BZTs  (Figure  5.1d),  and  DDT  and  its  common  metabolites  2,2- 
bis(4-chlorophenyl)- 1 , 1  -dichloroethane  (DDD),  2,2-bis(4-chlorophenyl)- 1 , 1  -dichloro- 
ethene  (DDE)  and  bis(4-chlorophenyl)acetic  acid  (DDA),  are  also  shown  (Figure  5.1e). 
All  of  these  HOCs  and/or  their  metabolites,  have  a  variety  of  functional  groups  and 
therefore  potential  reactive  sites  (“chemical  handles”)  to  bind  to  OM. 

5.2  MATERIALS  AND  METHODS 

5.2.1  Study  Area  and  Samples 

Geographic  locations,  dates,  depths  and  details  of  the  sites  examined  in  this  study 
are  provided  in  Table  5.1  and  Figure  5.2.  To  test  the  inverse  tracer  approach  on  a  sample 
known  to  contain  bound  HOCs,  soil  contaminated  with  mineral  oil  from  a  site  near 
Hamburg,  Germany  (SS),  previously  examined  by  Richnow  et  al.,  (1994)  was  obtained. 
Sediments  where  bound  HOCs  had  previously  been  reported  were  also  examined.  These 
included,  freshwater  sediments  contaminated  with  BZTs  from  the  Pawtuxet  River  (PR) 
(Reddy  et  al.,  2000a)  and  from  a  control  site  (Station  3),  upstream  of  the  chemical  plant 
that  released  the  BZTs  into  the  river  (Figure  5.2  insert),  and  freshwater  sediments 
contaminated  with  DDT  from  Teltow  Canal  (TC)  (Schwarzbauer  et  al.,  2003).  The  Palos 
Verdes  (PV)  shelf  is  a  coastal  marine  site  that  also  has  a  history  of  DDT  contamination 
(Eganhouse  et  al.,  2000)  and  thus  for  comparison,  sediments  from  this  site  were  analyzed. 
Creosote  contaminated  sediments  from  both  a  marine  site,  Eagle  Harbor  (EH),  and 
freshwater  site,  Chattanooga  Creek  (CC),  as  well  as  PCB  contaminated  sediments  from 
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Table  5.1  Locations,  dates,  depths  and  sedimentary  environments  of  samples  analyzed  in  this  study. 
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'nd  -  not  determined. 
PC  -  push-core. 

?BC  -  box-core. 

'  MC  -  multi-core, 
na  -  not  applicable. 
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Figure  5.2  Location  of  sediments  examined  in  this  study. 

New  Bedford  Harbor  (NBH)  and  Hudson  River  (HR)  were  examined.  Sites  contaminated 
with  petroleum  and  PAHs  were  sampled,  but  unlike  the  previous  investigation  of  anoxic 
sediments  from  West  Falmouth  (White  et  al.,  2005),  oxic  sediments  from  Paleta  Creek 
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(PC),  San  Francisco  Bay  (SF),  Boston  Harbor  (BH)  and  Winsor  Cove  (WC)  were  also 
analyzed.  Standard  Reference  Materials  1941b  “Organics  in  marine  sediment”  and  1649a 
“Urban  dust”  from  the  National  Institute  of  Standards  and  Technology  (NIST)  were 
selected  for  study  because  of  the  extensive  information  regarding  their  compositions  and 
broad  contamination  from  PAHs,  PCBs  and  pesticides,  as  well  as  providing  the 
opportunity  with  the  “Urban  dust”  sample  to  examine  the  significance  of  binding  in  a 
different  environmental  matrix  (atmospheric  particulates).  Data  from  West  Falmouth 
sediments  impacted  by  No.  2  fuel  oil  as  previously  described  (White  et  al.,  2005),  are 
also  included  in  this  study  for  comparison. 

5.2.2  Sequential  Chemical  Treatment 

Samples  were  air-dried  and  homogenized  with  a  mortar  and  pestle  before  passing 
through  a  1-mm  sieve.  Depending  on  the  quantity  available,  1  to  10  grams  of  each  sample 
was  extracted  with  a  mixture  of  dichloromethane  (DCM),  methanol  (MeOH)  and  water 
(H20)  (89:10:1)  by  pressurized  fluid  extraction  (100  °C,  1000  psi).  The  resulting  total 
lipid  extracts  (TLE)  were  acidified  with  4N  hydrochloric  acid  (HC1)  to  pH  4  to  ensure 
that  any  anionic  metabolites  would  be  protonated  before  back-extracting  into  DCM.  The 
TLE  was  then  reduced  in  volume  by  rotary  evaporation,  dried  over  Na2S04  and 
percolated  through  a  copper  column  plugged  with  glass  wool  to  remove  sulfur  and  any 
particulate  matter  prior  to  further  analysis.  The  remaining  solvent  extracted  residues  (EX- 
RES)  from  all  samples  were  then  refluxed  for  2  h  with  a  solution  of  0.5  N  potassium 
hydroxide  (KOH)  in  MeOH,  (100  mL)  and  H20  (20  mL).  After  cooling,  the  reaction 
mixture  was  separated  by  centrifugation  and  the  remaining  saponified  residue  (SA-RES) 
rinsed  with  MeOH,  followed  by  DCM,  and  then  hexane.  The  alkaline  solution  was  back 
extracted  with  hexane  (3  x  50  mL)  to  obtain  a  neutral  fraction  before  acidification  with 
4N  HC1  to  pH  2  and  extraction  with  DCM  (3  x  50  mL)  to  obtain  an  acidic  fraction.  Both 
fractions  were  reduced  in  volume  by  rotary  evaporation  and  dried  over  Na2S04. 
Combined,  these  fractions  comprise  the  saponification  extract  (SAE). 
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5.2.3  Demineralization  of  Pawtuxet  River  Sediments 

BZTs  had  previously  been  shown  to  be  present  in  the  nonextractable  fraction  of 
Pawtuxet  River  sediments,  and  were  hypothesized  to  be  linked  to  OM  via  an  ester  bond 
(Reddy  et  al.,  2000a).  To  examine  this  further  and  explore  the  role  of  minerals 
interactions  with  HOCs,  10  gram  aliquots  of  the  solvent-extracted  sediment  residues  were 
demineralized.  To  remove  carbonate  minerals,  the  extracted  sediments  were  placed  in 
60mL  acid-cleaned  Teflon  test  tubes  with  5  mL  of  12N  HC1  before  capping  the  tubes  and 
agitating  for  lhr.  Following  centrifugation,  the  supernatant  was  removed  and  this  was 
repeated  twice  more.  Silicate  minerals  were  then  removed  by  addition  of  10  mL  of 
concentrated  hydrofluoric  acid  (HF)  and  4  mL  of  HC1.  The  tubes  were  agitated  for  12  hrs, 
centrifuged  and  the  supernatant  removed.  This  was  repeated  twice  more  before  salting 
out  the  remaining  HF  with  boric  acid  and  several  water  rinses.  The  demineralized 
sediments  were  then  allowed  to  dry  before  adding  pyridine  (5  mL)  to  swell  the  sediments, 
followed  by  sonication  for  30  minutes.  The  tubes  were  then  centrifuged  and  the  pyridine 
extract,  which  was  now  a  dark  brown  color,  was  decanted  reduced  in  volume  and  dried 
over  Na2SC>4.  The  demineralized  sediment  was  then  saponified  with  KOH  in  MeOH  as 
previously  described  in  section  3.1. 

5.2.4  Bulk  Sediment  Analyses 

Aliquots  of  the  total  samples  prior  to  solvent  extraction  (TOC),  and  the  respective 
EX-RES  and  SA-RES  for  each  sample  were  treated  with  10%  HC1  to  remove  any 
carbonates.  These  bulk  samples  and  residues  were  then  analyzed  for  organic  carbon 
content  with  a  Carlo  Erba  1108  elemental  analyzer.  The  stable  carbon  isotope  ratio  (813C) 
and  radiocarbon  abundance  (A14C)  of  the  samples  and  respective  residues  were  measured 
on  purified  carbon  dioxide  (CO2)  after  their  combustion  in  the  presence  of  cupric  oxide 
(CuO).  Stable  carbon  values  were  determined  by  isotope  ratio  mass  spectrometry  (irMS) 
and  14C  content  by  accelerator  mass  spectrometry  (AMS)  at  the  National  Ocean  Sciences 
Accelerator  Mass  Spectrometry  (NOSAMS)  facility  at  Woods  Hole  Oceanographic 
Institution  (WHOI),  Woods  Hole,  MA.  Radiocarbon  analyses  were  performed  after 
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conversion  of  the  CO2  to  graphite  (McNichol  et  al.,  1994)  and  all  samples  were  run  as 
conventional  targets.  All  14C  measurements  are  normalized  to  8BC  values  of  -25%o  and 
expressed  as  A14C  values,  which  is  the  per  mille  (%o)  deviation  from  the  international 
standard  for  14C  dating,  Standard  Reference  Material  4990B  “Oxalic  Acid”.  Precision  for 
813C  and  A14C  measurements  are  ~0.1  and  l-5%o  (see  Appendix  for  errors  specific  to  each 
sample  analyzed),  respectively.  For  each  sample,  the  reported  A14C  =  \fme  (195°'x)X-l]1000 
(Stuiver  and  Polach,  1977),  where  X= 1/8267  (year'1),^n=  fraction  modem  14C,  corrected 
for  isotopic  fractionation  using  8I3C,  and  “x”  equals  the  year  of  collection  (see  Table  5.1). 

5.2.5  Gas  Chromatography-Mass  Spectrometry  (GC/MS)  Analysis 

Pawtuxet  River  sediment  extracts  were  acetylated  with  acetic  anhydride  and 
pyridine  and  the  resulting  product  was  solvent  exchanged  into  hexane  before  placing  on  a 
0.5cm  (i.d.)  x  15  cm  column  containing  activated  silica  gel.  The  first  fraction  was  eluted 
with  15  mL  of  a  mixture  of  hexane  and  DCM  (4:1)  and  the  second  fraction  containing 
BZT-acetates  was  eluted  with  15  mL  of  a  mixture  of  DCM  and  acetonitrile  (1:1).  These 
fractions  were  reduced  in  volume  by  rotary  evaporation  prior  to  further  analysis.  All  other 
extracts  were  derivatized  with  bis-(trimethylsilyl)-trifluoroacetamide  (BSTFA)  and 
pyridine  to  improve  chromatography.  Aliquots  of  extracts  from  all  samples  were  spiked 
with  an  internal  standard  of  W-C36  alkane  (4  pg).  For  compound  identification,  electron 
ionization  (El)  spectra  were  acquired  on  an  Agilent  6890  series  gas  chromatograph  (GC) 
interacted  to  an  Agilent  5973  mass  selective  detector  (MSD).  A  post  column  split  to  a 
flame  ionization  detector  (FID)  allowed  for  quantification  of  compounds  of  interest 
relative  to  the  internal  standard.  Extracts  were  separated  on  a  J&W  DB-5MS  capillary 
column  (60  m  x  0.32mm  ID,  0.25  pm  film  thickness)  with  helium  carrier  gas  at  a 
constant  flow  of  1  ml  min'1.  The  GC  oven  had  an  initial  temperature  of  40°C  (1  minute 
hold)  and  was  ramped  at  20°C  min'1  until  130°C  and  then  at  2°C  min'1  until  320°C  (30 
minute  hold).  Spectra  were  acquired  between  m/z  40-650  at  a  scan  rate  of  1  cycle  s'1. 
Compounds  present  in  the  extracts  were  identified  from  both  mass  spectra  and  gas 
chromatographic  retention  characteristics.  Contributions  of  an  unresolved  complex 
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mixture  (UCM)  to  some  of  the  extracts  were  quantified  by  integrating  the  total  FID  area 
of  the  UCM  and  using  response  factors  determined  from  No.  2  fuel  oil  standards  from  the 
Marine  Ecosystem  Research  Laboratory  (MERL)  as  in  Peacock  et  al.,  (2005). 
Saponification  extracts  from  HR  were  also  analyzed  by  electron  capture  negative 
ionization  (ECNI)  with  methane  as  a  reagent  gas  and  using  an  Agilent  6890  series  gas 
chromatograph  interfaced  to  an  Agilent  5973  mass  selective  detector  (MSD),  with  a  J&W 
Scientific  DB-XLB  column  (60  m  x  0.25mm  ID,  0.25  pm  film  thickness).  The  GC 
oven  had  an  initial  temperature  of  50°C  (1  minute  hold)  and  was  ramped  at  30°C  min'1 
until  120°C  and  then  at  8°C  min'1  to  300°C  (30  minute  hold).  Spectra  were  acquired 
between  m/z  50-750  at  a  scan  rate  of  1  cycle  s'1. 

5.3  RESULTS  AND  DISCUSSION 

5.3.1  Solvent-Extractable  Contaminants 

Variable  quantities  of  HOCs  were  extracted  from  the  samples  examined  in  this 
study  (Table  5.2).  Several  BZT  compounds  including  Cl-BZT,  C4-BZT,  C4-C1-BZT, 
C5-C1-BZT,  C8-BZT  and  C10-BZT  (Figure  5. Id)  were  isolated  from  Pawtuxet  River 
sediments  from  Stations  1  and  2.  Extracts  from  both  sediments  were  dominated  by  Cl- 
BZT  with  maximum  concentrations  of  6  pg  g'1  dry  weight  (dw)  sediment  (0.33  mg  g"1 
organic  carbon  (OC))  in  Station  1,  and  100  pg  g'1  dw  (1.9  mg  g'1  OC)  in  Station  2.  No 
BZT  compounds  were  identified  in  Station  3  sediments  due  to  it  being  upstream  from  the 
chemical  plant  that  originally  released  the  BZTs  into  the  river  (Figure  5.2),  and  it  was 
therefore  used  as  a  control  site.  Sediment  extracts  contaminated  with  PCBs  from  NBH 
were  dominated  by  di-,  tri-  and  tetrachlorobiphenyls,  but  also  contained  penta-  and 
hexachlorobiphenyls  with  a  summed  concentration  of  4200  pg  g"'dw  (37.5  mg  g"1  OC). 
Sediments  from  HR,  however  were  less  contaminated  and  mono-,  di-,  tri-  and 
tetrachlorobiphenyls  with  summed  concentrations  of  170  pg  g'1  dw  (7.7  mg  g'1  OC)  were 
present  in  the  TLE.  Sediments  from  TC  contained  a  variety  of  extractable  DDT 
metabolites  including  DDE,  DDD  and  DDA  (Figure  5.1e)  with  a  summed  concentration 
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of  740  (ig  g*'dw  (9.4  mg  g'1  OC),  whereas  PV  sediments  were  much  less  contaminated 
(0.002  jig  g"1  dw;  0.1  jig  g1  OC  of  DDE).  Creosote  contaminated  extracts  from  CC  and 
EH  were  dominated  by  PAHs  including  fluoranthene,  pyrene,  phenanthrene, 
benzo[6]fluoranthene,  benz [ci\ anthracene,  benzo[a]pyrene,  chrysene,  fluorene  and 
anthracene  with  individual  PAH  concentrations  between  10  to  50  jig  g'1  dw  (0.23  to  1.16 
mg  g"1  OC)  in  CC,  and  100  to  1400  jig  g1  dw  (9.1  to  127.3  mg  g'1  OC)  in  EH.  PAHs  also 
dominated  sediment  extracts  from  BH  with  fluoranthene,  pyrene  and  benzja] anthracene 
at  concentrations  of  1  to  2  jig  g'1  dw  (0.06  to  0.12  mg  g'1  OC)  as  well  as  SS,  which 
contained  phenanthrene,  anthracene,  pyrene,  fluoranthene,  chrysene  at  concentrations  of 
60  to  570  |ig  g_1dw  (2.1  to  19.7  mg  g'1  OC).  Details  of  solvent  extractable  contaminants 
and  their  concentrations  from  NIST  SRM  1941b  and  NIST  SRM  1649a  are  discussed  in 
the  certificates  of  analysis,  but  contribute  total  contaminant  concentrations  of  ~10  jig  g'1 
dw  (0.32  mg  g'1  OC)  and  ~60  jig  g'1  dw  (0.34  mg  g'1  OC)  respectively.  Varying 
quantities  of  UCM  with  different  carbon  ranges,  and  therefore  different  sources,  were 
also  detected  in  most  samples  except  TC,  SF  and  EH  (Table  5.2). 

Gas  chromatograms  of  TLE  from  select  sediments  including  NBH,  CC  and  EH 
indicate  the  dominance  of  HOCs  in  the  GC-resolvable  portion  of  the  TLE  (Figure  5.3a,  c 
and  d  respectively).  This  is  observed  for  the  majority  of  samples  and  in  general,  HOCs 
comprise  more  than  90%  of  the  resolvable  TLE  (Figure  5.4a).  Biogenic  hydrocarbons 
contribute  more  significantly  to  sediments  such  as  those  from  HR  (Figure  5.3b),  as  well 
as  PV  and  all  PR  sites  where  HOCs  constitute  50-90%  of  the  GC-amenable  TLE,  with  the 
remainder  consisting  of  biogenic  hydrocarbons  (Figure  5.4b).  In  some  samples,  such  as 
TC,  BH  and  SF,  biogenic  hydrocarbons  dominate,  and  HOCs  comprise  less  than  5% 
(Figure  5.4c).  Contributions  from  biogenic  hydrocarbons  to  the  TLE  consist  of  varying 
amounts  of  n-alkanoic  acids  from  n- Cm  to  H-C30  with  strong  even/odd  dominance  and  n- 
alkanols  from  n- C14  to  n-C 26  with  strong  even/odd  dominance.  The  short  chain  (n-Cn  to 
n- C22)  compounds  are  commonly  accepted  to  be  predominantly  derived  from  marine 
algae  or  bacteria  whereas  the  longer  chain  compounds  (n- C24  to  /1-C34)  have  terrestrial 
origins  from  higher  plant  epicuticular  waxes.  A  strong  odd/even  dominance  is  observed 
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Figure  5.3  Gas  chromatograms  of  total  lipid  extracts  from  a)  New  Bedford  Harbor,  b) 
Hudson  River,  c)  Chattanooga  Creek  and  d)  Eagle  Harbor.  The  dashed  line  is  the  GC 
baseline.  In  a)  PCBs  dominate  and  in  d)  PAHs  dominate.  In  b)  and  c)  contaminants  are 
designated  with  and  biogenic  hydrocarbons  w-alkanoic  acids,  n-alkanols,  alkanes 
and  sterols  are  labeled  with  “F”,  “O”,  “A”  and  “S”  respectively.  The  internal  standard 
hexatriacontane  is  designated  with  “/s”  and  an  external  contaminant  is  labeled  with  “e.” 
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Figure  5.4  Percent  contribution  of  biogenic  hydrocarbons  and  contaminants  (specified  in 
the  text)  to  the  total  lipid  extract. 
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for  n-alkanes  from  n- C23  to  n-C3i,  which  is  also  typical  of  higher  plants.  Sterol 
contributions  were  dominated  by  C27  to  C29  compounds  including  cholesterol,  a  general 
biomarker  for  plants  and  consumer  organisms  (Grice  et  al.,  1998),  and  (3-sitosterol  and 
stigmasterol,  which  are  the  major  sterols  in  higher  plants  (Figures  5.3  and  5.4). 

5.3.2  Nonextractable  Contaminants 

Chemical  evidence  for  the  incorporation  of  HOCs  into  nonextractable  residues  is 
observed  in  PR  Station  2,  TC,  CC9,  EH,  and  SS,  and  comprises  between  0.02  -  41%  of 
the  total  pool  of  HOCs  at  each  site.  At  PR  Station  2,  1  |xg  g_1of  Cl-BZT  (Figure  5.5a)  is 
released  by  alkaline  hydrolysis,  representing  0.4%  of  the  total  pool  of  BZTs  in  the 
sediment,  which  is  similar  to  the  quantity  observed  by  Reddy  et  al.,  2000a.  Known 
metabolites  of  DDT,  4,4’-dichlorobenzophenone  (4,4’-DBP)  and  2,2-bis(4- 
chlorophenyl)-ethanol  (4,4’-DOH)  (Figure  5.5b)  are  present  in  the  SAE  from  TC  at 
summed  concentrations  of  300  pg  g'1;  29%  of  the  total  HOC  pool.  Creosote-derived 
compounds,  fluoranthene,  pyrene  and  anthracene  carbonitrile  (Figure  5.5c)  were 
identified  in  the  SAE  from  CC  at  total  concentrations  of  7  pg  g1;  1%  of  the  total  HOC 
pool.  Analysis  of  the  SAE  from  EH  indicated  the  presence  of  9H-fluoren-9-one,  a  known 
metabolite  of  fluorene  (Grifoll  et  al.,  1992)  and  9,10-anthracenedione,  a  known 
constituent  of  creosote  (Figure  5.5d),  at  summed  concentrations  of  2  |ig  g'1;  0.02%  of  the 
total  pool  of  HOCs  at  the  site.  In  SS,  fluoranthene,  pyrene,  anthracene  carbonitrile 
(Figure  5.5c),  acridine,  and  aminofluoranthene  (Figure  5.5e)  were  identified  with 
summed  concentrations  of  180  pg  g'1;  0.02%  of  the  total  pool  of  HOCs  at  the  site. 

Significant  quantities  of  biogenic  hydrocarbons  are  also  released  as  a  result  of 
alkaline  hydrolysis  as  evidenced  in  chromatograms  of  the  acidic  (Figure  5.6)  and  neutral 
(Figure  5.7)  extracts  from  saponification  of  NBH,  HR,  CC  and  EH.  The  acidic  extracts 
are  dominated  by  H-alkanoic  acids  from  «-Ci2  to  M-C28  with  strong  even/odd  dominance 
and  the  neutral  extracts  are  dominated  by  n-alkanols  from  n- C12  to  n- C22  with  strong 
even/odd  dominance  and  sterols  as  before.  Nonextractable  HOC  residues  are  only 
observed  in  the  neutral  extracts  (Figure  5.7c  and  d)  and  the  presence  of  alkanes  in  the 
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acidic  extract  from  HR  (Figure  5.6d)  is  likely  an  external  contaminant  because  if  alkanes 
were  truly  present,  they  would  partition  into  the  neutral  fraction.  Overall  the  contribution 
of  biogenic  hydrocarbons  to  the  GC-amenable  portion  of  the  SAE  is  between  93  and 
100%  and  is  dominated  by  short  chain  «-alkanoic  acids  (Figure  5.8).  Where  detected, 
HOCs  and/or  their  metabolites  contribute  a  maximum  of  1-7%  of  the  resolvable  SAE 
(Figure  5.8a),  but  for  the  majority  of  samples,  nonextractable  residues  of  HOCs  are  not 
observed  by  GC-FDD  (Figure  5.8b). 


Figure  5.5  Structures  of  hydrophobic  organic  contaminants  released  by  alkaline 
hydrolysis  from  various  samples  examined  in  this  study.  Compounds  in  d),  e)  were 
tentatively  identified  by  GC/MS. 
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Figure  5.6  Gas  chromatograms  of  saponification  acidic  extracts  from  a)  New  Bedford 
Harbor,  b)  Hudson  River,  c)  Chattanooga  Creek  and  d)  Eagle  Harbor.  The  dashed  line 
is  the  GC  baseline.  In  all  samples,  the  w-alkanoic  acids  are  numbered.  The  n-alkanoic 
acids  in  d)  have  different  retention  times  as  they  were  methylated  in  the  saponification 
process.  The  internal  standard  hexatriacontane  is  designated  with  “/s”. 
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Figure  5.7  Gas  chromatograms  of  saponification  neutral  extracts  from  a)  New  Bedford 
Harbor,  b)  Hudson  River,  c)  Chattanooga  Creek  and  d)  Eagle  Harbor.  The  dashed  line  is 
the  GC  baseline.  In  all  samples,  the  n-alkanols  acids  are  numbered.  Contaminants  are 
designated  in  c)  and  d)  with  The  internal  standard  hexatriacontane  is  designated  with 
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Figure  5.8  Percent  contribution  of  biogenic  hydrocarbons  and  contaminants  (specified 
the  text)  to  the  saponification  extract. 


5.3.3  Nature  of  Nonextractable  Contaminants 

Possible  interactions  between  HOCs  and  nonextractable  sediment  residues  can  be 
inferred  from  the  chemical  structure  of  the  HOCs  present  in  the  SAE  (Figure  5.5),  as  well 
as  the  chemical  treatment  by  which  they  are  released.  Alkaline  hydrolysis  is  used  to 
target  ester  linkages  liberating  products  with  hydroxyl  or  carboxyl  functional  groups,  and 
thus  for  the  HOCs  released  in  this  study  (Figure  5.5),  only  Cl-BZT  and  4,4’DDOH  are 
likely  to  have  been  associated  via  an  ester  bond.  Demineralization  of  the  sediment  from 
PR  Station  2,  however,  liberated  the  same  quantity  (1  Jig  g"1  dw)  of  Cl-BZT  as  that 
released  by  alkaline  hydrolysis,  and  subsequent  alkaline  hydrolysis  of  the  demineralized 
sediment  did  not  release  more.  This  indicates  that  Cl-BZT  is  unlikely  to  be  bound  via  an 
ester  bond  to  OM  as  previously  suggested  (Reddy  et  al.,  2000a)  and  is  instead  likely  to  be 
associated  with  minerals.  The  previous  suggestion  that  Cl-BZT  was  bound  to  OM  via  an 
ester  linkage  was  based  on  the  finding  that  it  was  the  only  BZT  detected  in  significant 
quantities  in  the  bound  fraction  (as  much  as  8%  of  the  total  BZT  concentration)  compared 
to  C8-C1-BZT  (0.03%  of  the  total  BZT  concentration)  and  C10-BZT,  which  was  never 
detected  in  the  bound  fraction  (Reddy  et  al.,  2000a).  This  was  thought  to  result  from 
steric  hindrance  by  the  presence  of  alkyl  groups  at  the  3’  position  for  C8-C1-BZT  and 
C10-BZT  (see  figure  5.1  for  structures),  which  would  inhibit  binding  to  OM.  If  it  is  to  be 
concluded  that  Cl-BZT  is  instead  associated  with  minerals,  then  the  preferential 
association  of  this  compound  may  arise  from  the  fact  that  its  smaller  size  allows  it  to  be 
more  easily  incorporated  into  mineral  structures.  For  example,  naturally  occurring 
minerals  such  as  anionic  clays,  also  known  as  layered  double  hydroxides,  may 
incorporate  organic  compounds  between  the  layers  of  the  clay,  and  these  layers  can 
impose  a  restricted  geometry  on  their  interlayer  guests  (Wei  et  al.,  2005).  Alternatively, 
the  influence  of  the  3’  alkyl  group,  present  in  C8-C1-BZT  and  C10-BZT  may  have  an 
affect  on  the  acid  dissociation  of  the  adjacent  hydroxy  group  (as  suggested  by  Reddy  et 
al.  2000a),  which  may  then  decrease  the  potential  for  anionic  interactions  with  minerals. 

It  is  possible  that  aminofluoranthene  may  have  been  bound  to  OM  via  a  carbon- 
nitrogen  (C-N)  bond  as  oxidative  coupling  may  occur  between  OM  and  compounds  with 
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amino  functional  groups  (Sjobald  and  Bollag,  1977).  The  carbon-nitrogen  bond  formed 
between  the  aminofluoranthene  and  OM  in  sediment  would  be  cleaved  by  base  hydrolysis 
to  liberate  the  amino  compound  and  an  aldehyde  or  ketone  group  on  the  OM  (March, 
1992).  It  is  possible  that  aminofluoranthene  may  be  physically  entrapped  within  the 
nonextractable  residues  in  sediment,  and  in  general  this  appears  to  be  the  most  likely  fate 
for  the  other  HOCs  (Figure  5.5)  that  do  not  possess  functional  groups  through  which  they 
could  be  bound.  Alkaline  hydrolysis  releases  entrapped  HOCs  by  hydrolysis  and 
sufficient  breakdown  of  the  sedimentary  macromolecular  OM  structure  (Northcott  and 
Jones,  2000). 

5.3.4  Radiocarbon  as  an  Inverse  Tracer  of  Contaminants 

Assuming  that  all  contaminants  are  solvent  extractable,  the  relative  contribution 
of  solvent-extractable  fossil  fuel-derived  HOCs  to  the  samples  examined  in  this  study  can 
be  calculated  based  on  an  isotope  mass  balance  (see  White  et  al.,  2005)  where  in  terms  of 
radiocarbon  abundances  (A), 

f  =  (At  —  Ah)/(Ax  —  An)  (I) 

where  f  is  the  fraction  of  the  mixture  that  is  fossil  fuel-derived,  At  is  the  decay-corrected 
A14C  value  for  the  untreated  sediment,  An  is  the  corrected  A14C  value  for  the  extracted 
sediment  residue  (Table  5.3)  and  Ax  is  the  A14C  of  fossil  fuel-derived  HOCs  and  reflects 
an  absence  of  14C  (i.e.,  A14C  =  -1000%o).  Solutions  of  /  were  obtained  for  the  samples 
examined  in  this  study,  and  are  expressed  as  the  percentage  of  TOC  that  is  fossil  fuel- 
derived  (Table  5.2).  These  values  are  compared  to  those  calculated  from  resolvable 
contaminant  and  UCM  concentrations  (measured  by  GC-FID;  Table  5.2)  relative  to  TOC 
content  (determined  by  elemental  analysis;  Table  5.3).  The  mass  of  UCM  was  converted 
to  carbon  equivalents  considering  the  mass  contribution  of  hydrogen  to  the  hydrocarbon, 
which  in  an  alkane  typical  of  the  UCM  («-C10  to  n-C38  as  before)  is  approximately  15% 
of  the  total.  Similar  conversions  of  contaminant  quantities  to  carbon  equivalents  were 
performed  for  the  other  contaminants  examined  considering  non-carbon  elements  such  as 
oxygen,  nitrogen,  chlorine  and  hydrogen.  On  average,  95%  of  PAHs,  70%  of  BZTs,  60% 
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of  PCBs  and  50%  of  DDT  and  its  metabolites  are  derived  from  carbon.  To  evaluate  the 
quantity  of  bulk  OC  that  these  compounds  represent,  the  percent  of  the  TOC  that  is 
extractable  was  calculated  by  weighing  a  portion  of  the  TLE  and  converting  this  weight 
to  carbon  equivalents.  The  percent  of  the  TLE  that  is  GC  amenable  was  also  calculated 
(Table  5.4). 

For  the  majority  of  sites  examined,  there  is  reasonable  agreement  between  the 
percent  TOC  that  is  fossil  fuel-derived  determined  by  GC-FID  compared  to  that 
calculated  by  isotope  mass  balance  (Table  5.2).  To  directly  compare  these  approaches, 
the  difference  in  A14C  between  the  untreated  sediment  (A14Ctoc)  and  extracted  sediment 
residue  (A14Cex-res)  is  plotted  against  percent  TOC  that  is  fossil  fuel-derived,  as 
measured  by  GC-FID  (Figure  5.9).  The  difference  between  A14Ctoc  and  A,4Cex-res 
predominantly  represents  the  removal  of  fossil  organic  carbon,  assuming  that  the  A14C  of 
biogenic  hydrocarbons  contributing  to  the  TLE,  is  the  same  as  the  uncontaminated 
A14Ctoc  from  which  it  was  extracted.  In  general,  as  the  difference  between  A14Ctoc  and 
A14Cex-res  becomes  more  negative,  reflecting  greater  removal  of  fossil  fuel-derived 
carbon,  the  contribution  of  fossil  fuel-derived  carbon  measured  by  GC-FID  also 
increases.  This  inverse  relationship  (first,  considering  EH  as  an  outlier)  is  indicated  by  a 
solid  line  in  Figure  5.9.  The  quantity  of  fossil  fuel-derived  carbon  calculated  by  isotope 
mass  balance  and  GC-FID  for  EH  is  similar  (77  and  72%  respectively;  Table  5.2)  and 
thus  it  may  be  incorrect  to  consider  it  as  an  outlier.  If  we  instead  accept  the  EH  data,  then 
significant  differences  emerge  between  these  two  methods  for  samples  from  PY,  PR 
Station  2,  NBH,  PC,  HR,  SS  (Table  5.2;  circled  in  Figure  5.9).  Specifically,  the  quantity 
of  fossil  fuel-derived  carbon  measured  by  GC-FID  at  these  sites  is  underestimated 
relative  to  that  determined  by  isotope  mass  balance.  Excluding  these  samples  reveals  an 
inverse  relationship  with  a  strong  correlation  (r2^  0.96,  p=9.28xl0'8). 

For  sites  above  the  regression  line  such  as  TC  and  BH,  the  %TOC  that  is  fossil 
fuel-derived  as  calculated  by  mass  balance  is  less  than  that  measured  by  GC-FID.  In 
these  sites,  the  quantity  of  contaminant  is  small  (0.5  and  0.02%,  respectively)  and  the 
lowest  observed  in  the  range  of  samples  examined  (Table  5.2).  The  negative  values  of 
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Table  5.3  Percent  organic  carbon,  stable  catbon  and  radiocarbon  abundances0  for  bulk  samples  and  their  respective  residues. 
"Sample  %  OC  %  OC  %  OC  Wc  5^C  jj^C  A^C  A^C  A**c" 
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TOC  is  the  untreated  sample  prior  to  solvent  extraction. 

EX-RES  is  the  extracted  residue  remaining  after  solvent  extraction. 
SA-RES  is  the  saponified  residue  remaining  after  saponification. 
Data  for  this  sample  taken  from  White  et  al,  (2005). 


%TOC  that  is  fossil  fuel-derived  as  calculated  by  mass  balance  at  these  sites  and  SF, 
arises  from  the  A14Ctoc  being  more  enriched  than  the  A14Cex-res,  which  again  is  only 
observed  at  these  sites  (Table  5.3).  This  likely  occurs  as  a  result  of  the  removal  of 
biogenic,  solvent  extractable  14C  enriched  material  which  leaves  behind  a  relatively 
depleted  EX-RES.  This  is  supported  by  GC  analyses  of  the  TLE  for  these  sites,  which  are 
dominated  by  biogenic  hydrocarbons  (Figure  5.4c). 

Understanding  why  the  quantity  of  fossil  fuel-derived  carbon  is  underestimated  by 
GC-FID,  however,  is  more  complex  as  these  samples  are  contaminated  with  different 
HOCs  of  varying  concentrations  and  include  soil  as  well  as  marine  and  freshwater 
environments.  The  commonality  between  them  is  that  they  contain  significantly  more 
fossil  fuel-derived  carbon  than  the  sites  that  fit  the  trend  (19-52%  compared  to  0-10% 
respectively  by  mass  balance;  Table  5.2).  Good  agreement  between  the  percent 
contribution  of  petroleum  to  the  TOC  calculated  by  GC-FID  compared  to  isotope  mass 
balance  was  observed  in  the  study  of  West  Falmouth  sediments  (White  et  al.,  2005;  a 
representative  contaminated  sediment  horizon  is  shown  in  Figure  5.9),  where  the 
contribution  of  petroleum  was  always  less  than  10%.  The  exception  to  this,  however,  is 
EH,  which  is  the  most  heavily  contaminated  sample,  but  displays  good  agreement 
between  the  percent  contribution  of  fossil  fuel-derived  contaminants  determined  by  GC- 
FID  and  mass  balance  (72  and  77%  respectively;  Table  5.2).  Compared  to  other  samples, 
EH  is  composed  of  coarse  sandy  material  and  has  a  low  OC  content  when  the  creosote  is 
removed  (0.3%),  whereas  other  sites  have  OC  contents  of  0.5  -  13.4%  (Table  5.3). 
Organic  carbon  content  may  affect  the  ability  to  correctly  determine  the  %TOC  that  is 
fossil  fuel-derived  by  GC-FID  if  sediments  are  organic  carbon  rich,  and  large  quantities 
of  extractable  biogenic  hydrocarbons  lessen  the  impact  of  14C-dead  fossil  carbon  or  mask 
the  contribution  of  HOCs  observable  by  GC.  The  outliers  in  Figure  5.5,  however,  have 
varying  OC  contents  (1.6  -  1 1.2%;  Table  5.3)  and  a  range  in  quantities  of  %TOC  that  is 
extractable  (3  -  67%;  Table  5.4)  and  thus  this  does  not  appear  to  be  a  controlling  factor. 
The  quantity  of  the  TLE  that  is  GC  amenable,  however,  plays  a  more  significant  role.  In 
the  case  of  EH,  77%  of  the  TLE  is  GC-amenable,  but  for  the  samples  that  do  not  fit  the 
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%  TOC  fossil  fuel-derived 

Figure  5.9  Difference  between  A14C  of  unextracted  sediment  (TOC)  and  solvent- 
extracted  sediment  residue  (EX-RES)  as  a  function  of  percent  TOC  that  is  comprised  of 
solvent-extractable  fossil-fuel-derived  contaminants  measured  by  gas  chromatography. 


trend,  this  value  is  much  less  (32  -  61%;  Table  5.4).  To  examine  this  further,  the 
deviation  from  the  linear  regression  in  Figure  5.5  is  plotted  against  %TLE  that  is  non-GC 
amenable  (Figure  5.10).  An  inverse  relationship  with  strong  correlation  (i*=0.82, 
p=0.00004)  is  observed.  This  regression  excludes  BH,  CC,  and  WC  samples  that  exhibit 
a  high  percentage  of  non-GC-amenable  TLE  (67  -  93%),  but  display  reasonable 
agreement  between  the  %TOC  fossil  fuel-derived  calculated  by  the  two  methods.  Fossil 
fuel-derived  HOCs  and/or  metabolites  that  are  extractable  may  not  be  GC-amenable  if 
they  are  bound  within  the  TLE.  This  solvent  extractable,  but  chemically  bound  pool  may 
be  quantitatively  significant  and  requires  further  investigation. 
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Deviation  from  relationship  between  %  TOC  fossil  fuel- 
derived  calculated  by  GC-FID  and  isotope  mass  balance 


Figure  5.10  Difference  between  the  percent  of  total  lipid  extract  (TLE)  that  is  non-GC 
amenable  and  the  deviation  from  the  relationship  between  the  percent  total  organic 
carbon  (%TOC)  that  is  fossil-fuel  derived  calculated  by  GC-FID  and  isotope  mass 
balance. 


Additional  fossil  inputs  to  the  sediments  that  are  not  contaminant  derived,  but 
instead  come  from  “pre-aged”  natural  organic  carbon  sources  such  as  remobilization  of 
soil  carbon  and  weathering  of  sedimentary  rocks  (Masiello  and  Druffel,  2001;  Blair  et  al., 
2004;  Gofii  et  al.,  2005),  may  also  play  a  role  in  the  underestimation  of  %TOC  that  is 
fossil  fuel-derived  as  calculated  by  GC-FID.  There  is  no  obvious  molecular  evidence  for 
contributions  of  these  sources  in  the  SAE,  but  “pre-aged”  carbon  may  be  metabolized  by 
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microbes  and  incorporated  into  biogenic  lipids  (Petsch  et  al.,  2001)  that  would  not  be 
incorporated  into  the  GC-FID  quantification  of  fossil  carbon  derived  from  HOCs. 
Because  of  this  potential  uncontaminated  control  sites  of  similar  depositional  setting 
provide  a  useful  comparison.  This  was  possible  for  the  PR  sites  where  an  upstream  site 
(Station  3)  was  compared  to  a  contaminated  site  downstream  (Station  2).  The  A14Ctoc 
was  different  at  the  two  sites  owing  to  the  more  significant  contamination  at  Station  2 
(8.41  mg  g"1,  A14Ctoc  =  -354%o)  than  at  Station  3  (0.79  mg  g"1,  A14CToc  =  -292%o),  but 
once  the  majority  of  contaminant  was  removed  the  remaining  A14Cex-res  of  Stations  2 
and  3  are  very  similar  (-187%o  and  -192%o  respectively)  indicating  that  most  of  the  fossil 
fuel-derived  carbon  has  been  removed  by  solvent  extraction. 

5.3.5  Isotopic  Evidence  of  Nonextractable  Contaminants 

To  examine  the  quantity  of  fossil  fuel-derived  carbon  removed  by  alkaline 
hydrolysis,  equation  (I)  as  described  before  is  applied,  but  instead,  A t  is  the  corrected 
A14C  values  for  EX-RES  and  An  is  the  corrected  A14C  values  for  the  SA-RES  (Table  5.3). 
Although  bound  contaminants  are  observed  in  PR  St.  2,  TC  and  CC9,  as  previously 
described,  this  is  not  observable  by  mass  balance  as  they  contribute  only  0.002%,  0.4%, 
0.02%,  respectively,  of  OC  (calculated  by  GC-FID)  to  the  EX-RES.  This  is  insufficient 
to  be  detected  by  isotope  mass  balance  (see  earlier  discussion  on  sensitivity  of  approach). 
Samples  from  SS  and  EH,  do,  however,  exhibit  isotopic  evidence  for  the  presence  of 
nonextractable  (i.e.  ester-bound)  contaminant  residues.  For  SS,  isotope  mass  balance 
indicates  that  1.2%  of  OC  in  the  EX-RES  is  fossil  fuel-derived,  which  agrees  well  with 
that  obtained  by  comparing  the  amount  quantified  by  GC-FID  to  the  OC  content  of  EX- 
RES  (determined  by  elemental  analysis;  Table  5.3).  At  EH,  however,  isotope  mass 
balance  indicates  greater  formation  of  nonextractable  residues  than  that  measured  by  GC- 
FID,  5%  compared  to  0.1%  respectively.  This  discrepancy  may  be  due  to  the  reasons 
discussed  previously  as  only  12%  of  the  SAE  is  GC  amenable  at  this  site  (Table  5.4). 

Analysis  of  the  HR  saponification  extracts  by  GC-FID  did  not  reveal  any 
detectable  contaminants  (Figure  5.6b  and  5.7b)  or  potential  metabolites,  but  isotope  mass 
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balance  indicates  that  25%  of  OC  from  the  EX-RES  is  fossil  fuel-derived.  Because  this 
site  is  contaminated  with  PCBs,  electron  capture  negative  ionization  gas  chromatography 
was  used  to  probe  for  chlorinated  compounds.  The  major  chlorinated  compound  detected 
in  the  SAE  was  octachlorodibenzofuran  (OCDF;  Figure  5.5),  at  a  concentration  of  ~1  |ig 
g"1,  which  does  not  account  for  all  of  the  fossil  carbon  detected  by  mass  balance,  likely 
due  to  the  reasons  previously  discussed  for  solvent  extractable  HOCs.  Furans  are 
generally  unwanted  byproducts  from  the  manufacturing  of  other  chemicals  such  as  some 
disinfectants,  wood  preservatives  and  herbicides,  as  well  as  being  emitted  during 
combustion  processes  such  as  the  incineration  of  municipal  and  industrial  waste,  wood 
and  gasoline  burning.  Polychlorinated  dibenzo-p-dioxins  (PCDD)  have  previously  been 
observed  in  geological  formations  that  are  free  of  human  impact  (Ferrario  et  al.,  2000; 
Rappe  et  al.,  2000),  suggesting  that  these  compounds,  that  are  quite  similar  in  structure  to 
dibenzofuans  (PCDF),  may  be  formed  in  nature.  Whilst  PCDDs  and  PCDFs  are  observed 
in  surface  sediments,  corresponding  pre-industrial  coastal  sediments  (Hashimoto  et  al., 
1990)  and  lake  sediments  (Isosaari  et  al.,  2002)  detect  only  PCDDs  in  pre-industrial 
horizons  and  note  that  OCDF  is  below  detection  limits.  This  suggests  therefore  that 
OCDF  is  most  likely  to  be  anthropogenically-derived.  The  structure  of  OCDF  does  not 
suggest  that  it  was  chemically  bound  to  OM,  but  instead  was  likely  to  be  physically 
associated. 

At  sites  where  significant  quantities  of  nonextractable  HOCs  are  detected  (HR, 
EH  and  SS)  the  A14Csa-res  is  more  enriched  than  A14Cex-res  (Table  5.3),  due  to  the 
removal  of  this  fossil  (14C-dead)  carbon.  At  all  other  sites,  the  A14Csa-res  is  more 
depleted  than  the  A14Cex-res  (Table  5.3),  indicating  the  removal  of  relatively  14C  enriched 
material.  These  observations  are  summarized  in  Figure  5.11  where  the  difference  in  A14C 
isotopic  composition  between  the  two  residues  (AA14Cex-res  -  A14Csa-res)  is  plotted 
against  the  difference  in  813C  isotopic  composition  between  the  two  residues  (A813Cex- 
res  -  613Csa-res)-  The  majority  of  samples  examined  in  this  study  exhibit  AA14CEx-res  - 
A14Csa-res  differences  >1  indicating  the  removal  of  modem  organic  carbon.  This  has 
previously  been  observed  (White  et  al.,  2005)  and  similarly,  the  saponified  extracts  are 
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mainly  comprised  of  intrinsically  labile,  but  chemically  protected  biogenic  compounds 
dominated  by  short  chain  «-alkanoic  acids  (Figure  5.8).  Differences  in  AA14Cex-res  - 
A14Csa-res  <1  are  only  observed  for  SS,  EH  and  HR  samples  and  reflect  the  removal  of 
varying  degrees  of  fossil  fuel-derived  HOCs. 


removal  of  modern 
organic  carbon 


removal  of  fossil 
organic  carbon 


Figure  5.11  Difference  between  A14C  of  solvent-extracted  (EX-RES)  and  saponified 
(SA-RES)  sediment  residues  plotted  against  the  difference  between  813C  of  solvent- 
extracted  (EX-RES)  and  saponified  (SA-RES)  sediment  residues. 


By  comparison,  the  A813Cex-res  -  813Csa-res  is  much  less  and  only  a  ±l%o 
difference  is  observed.  In  general  differences  in  813C  of  the  bulk  sediments  and  residues 
are  not  observed  as  a  result  of  the  sequential  treatments.  This  is  likely  due  to  the  fact  that 
there  is  no  significant  isotopic  difference  between  the  biogenic  hydrocarbons  and  HOCs. 
The  813C  of  many  of  the  HOCs  examined  in  this  study  has  been  measured,  for  example, 
petroleum  derived  UCM  is  typically  ~-27%o  (Reddy  et  al.,  2002 b),  Aroclor  mixtures  of 
PCBs  range  from  -26.6  -  25.7%o  and  DDT  is  -28.3%o  (Drenzek  et  al.,  2002).  For  samples 
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from  NBH,  PC,  BH,  CC,  EH,  WC  and  1649,  the  813Ctoc  is  more  depleted  than  the 
513Cex-res  reflecting  the  contribution  of  these  HOCs  that  are  likely  more  depleted  than 
the  uncontaminated  TOC  represented  by  the  EX-RES.  At  PV,  PR  all  stations,  1941,  TC 
and  SF  the  813Cxoc  and  813Cex-res  are  the  same  within  error  (0.1%o)  indicating  that  unlike 
14C,  shifts  in  813C  are  too  insensitive  to  provide  evidence  for  the  presence  of  these  HOCs. 

5.4  CONCLUSIONS 

In  general  the  contribution  of  HOCs  to  the  TOC  in  sediments  as  calculated  by 
GC-FID  shows  good  agreement  with  that  calculated  by  isotope  mass  balance.  Offsets 
between  the  two  approaches  mostly  arise  from  the  fact  that  not  all  of  the  fossil  fuel- 
derived  carbon  in  the  extracts  is  GC-amenable.  Overall,  alkaline  hydrolysis  does  not 
release  significant  quantities  of  HOCs,  but  exceptions  exist  and  up  to  40%  of  the  total 
pool  of  DDT  and  its  metabolites  are  associated  with  nonextractable  residues  in  TC 
sediments.  This  quantity,  however,  only  comprises  a  small  fraction  of  the  TOC  and  thus 
the  presence  of  these  nonextractable  HOCs  is  not  detectable  by  isotope  mass  balance. 
Conversely,  in  HR  sediments,  25%  of  the  extracted  sediment  residue  is  comprised  of 
fossil  carbon  as  calculated  by  mass  balance,  but  there  is  little  evidence  by  GC/MS  of  the 
HOCs  that  constitute  this  pool.  It  is  therefore  unclear  whether  this  fossil  carbon  is  derived 
from  HOCs  or  from  natural  pre-aged  carbon  inputs. 

The  lack  of  significant  chemical  binding  of  HOCs  to  OM  in  aquatic  sediments 
may  arise  from  the  absence  of  a  suitable  chemical  environment  for  this  to  occur. 
Sediments  are  water  saturated  and  thus  condensation  reactions  that  form  ester  linkages 
are  energetically  unfavorable  and  reactive  species  such  as  hydroxyl  radicals  are  less 
present  than  they  are  in  a  soil  environment  for  example.  Binding  of  HOCs  to  sedimentary 
OM  is  therefore  more  likely  to  be  dependent  on  mediation  by  microorganisms  and  the 
lack  of  binding  in  these  sites  reflects  an  inefficiency  of  microbial  populations  to  mediate 
these  processes  that  occur  in  soil  environments. 


166 


It  is  important  to  note  that  the  sensitivity  of  the  inverse  14C  approach  is  linked  to 
the  TOC  content  of  the  sediment.  For  example,  in  sediments  where  TOC  and  extractable 
OM  are  high,  even  if  absolute  concentrations  of  HOC  are  also  high,  they  may  be  diluted 
by  abundant  and  relatively  modem  carbon.  Conversely  in  sediments  where  TOC  contents 
are  low,  HOCs  may  be  more  easily  detected.  This  study,  however,  has  examined  a  variety 
of  sediments  with  wide  ranging  TOC  contents  and  binding  of  HOCs  has  not  been  found 
to  be  significant.  Future  work  examining  contaminated  sediments  with  low  TOC 
contents,  or  those  where  binding  has  previously  been  observed  could,  however,  prove  this 
more  definitively. 

Chemical  binding  between  HOCs  and  OM  enables  HOCs  to  persist  in  sediments, 
but  in  this  form,  they  lose  their  chemical  identity  making  them  less  bioavailable  and  less 
toxic  (Gevao  et  al.,  2000).  HOCs  are  clearly  persistent  in  the  sites  examined  in  this 
study,  but  this  is  not  due  to  the  presence  of  covalent  bonds  between  HOCs  and  OM  and 
therefore  the  majority  of  HOCs  are  freely  associated  in  their  original  and  potentially  toxic 
forms. 
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Teltow  Canal  (TC)  SA-RES  1  mg  OS-44178  -27.7  0.4715  0.0026  -531.6 
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Boston  Harbor  (BH)  SA-RES  1  mg  OS-45701  -22.8  0.4661  0.0022  -536.9 

Winsor  Cove  (WC)  TOC  1  mg  OS-49326  -17.7  1.0822  0.0046  75.0 


Winsor  Cove  (WC)  EX-RES  1  mg  OS-49327  -17.1  1.1000  0.0049  92.7 

Winsor  Cove  (WC)  SA-RES  1  mg  OS-49189  -17.1  1.0968  0.0046  89.5 

Eagle  Harbor  (EH)  TOC  1  mg  OS-41947  -23.6  0.1316  0.0015  -869.2 

Eagle  Harbor  (EH)  EX-RES  1  mg  OS-41948  -22.4  0.5788  0.0028  -424.9 
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a  A  C  values  displayed  in  table  are  raw  data  and  are  not  corrected  for  in-situ  C  decay,  TOC  is  the  untreated  sediment 
prior  to  solvent  extraction, c  EX-RES  is  the  extracted  sediment  residue, d  SA-RES  is  the  saponified  sediment  residue. 
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CHAPTER  6 


ISOTOPIC  CONSTRAINTS  ON  THE  PROVENANCE  OF  INDIVIDUAL 
BIOMARKERS  IN  SALT  MARSH  SEDIMENTS 


ABSTRACT 


Natural  abundance  radiocarbon  contents  (A14C)  of  individual  short  (m-Cm-is)  and 
long-chain  («-C24-28)  fatty  acids  and  long-chain  (H-C26-C32)  alkanes  from  salt  marsh 
sediments  impacted  by  an  oil  spill  in  1969  were  measured  to  determine  their  provenance. 
The  oil  spill  represents  a  pulsed  introduction  of  14C-free  organic  matter,  however,  the 
timing  of  the  oil  spill  coincided  with  a  period  of  intense  above-ground  testing  of 
thermonuclear  weapons,  and  thus  salt  marsh  sediments  deposited  during  this  period  are 
also  influenced  by  an  excess  quantity  of  14C.  Compound  specific  isotope  analysis  (CSIA) 
of  free  and  bound  short-chain  fatty  acids  displays  14C  abundances  that  reflect  marine 
algal  and  bacterial  inputs  to  the  marsh  sediments,  while  isotopic  characteristics  (13C,  14C) 
of  long-chain  fatty  acids  are  consistent  with  an  origin  from  the  C4  salt  marsh  cord  grass 
Spartina  alterniflora.  Molecular  distributions  and  A14C  and  813C  values  of  n-alkanes 
indicate  that  they  are  predominantly  biogenic  with  lesser  contributions  from  pyrogenic 
sources  (late  20th  century  horizons).  Even  within  the  oil  contaminated  sediment  horizon, 
14C  data  unequivocally  show  that  petroleum-derived  w-alkanes  are  absent,  demonstrating 
that  these  compounds  have  been  completely  removed  through  processes  such  as 
evaporation,  dissolution  and/or  biodegradation.  In  contrast,  the  persistence  of  «-alkanes 
equivalent  in  carbon  number,  but  instead  derived  from  biogenic  and  pyrogenic  sources, 
indicates  that  their  resistance  to  degradation  is  most  likely  due  to  a  form  of  physical 
protection  or  association  with  the  sedimentary  matrix. 
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6.1  INTRODUCTION 


To  further  our  understanding  of  the  cycling  of  organic  carbon  in  salt  marsh 
sediments,  and  to  complement  previous  studies  examining  the  persistence  of  petroleum 
residues  in  these  environments  (Reddy  et  al.,  2002;  White  et  al.,  2005a;  Slater  et  al., 
2005),  compound  specific  isotope  analysis  (CSIA)  of  individual  solvent-extractable  and 
saponifiable  biomarker  lipids  including  long-chain  (71-C26-32)  alkanes  and  individual  short 
(w-Ch-is)  and  long-chain  (/1-C24-28)  fatty  acids  has  been  performed.  These  compounds 
were  isolated  from  different  sediment  intervals,  including  a  horizon  that  contains 
degraded  oil  residues  from  a  spill  of  No.  2  fuel  oil  that  occurred  and  inundated  the  salt 
marsh  in  1969  (Sanders  et  al.,  1980  and  references  therein).  Persistent  oil  residues  are 
present  in  the  form  of  an  unresolved  complex  mixture  (UCM)  of  hydrocarbons  in  the 
boiling  range  of  «-Ci3  to  «-C24,  consistent  with  a  biodegraded  No.  2  fuel  oil  source 
(Reddy  et  al.,  2002).  Individual  biomarkers  isolated  from  uncontaminated  sediment 
horizons  at  depths  above  and  below  the  oiled  sediments  were  examined  to  compare  and 
contrast  the  sources  and  fates  of  both  natural  and  anthropogenically  derived  carbon  inputs 
to  a  salt  marsh.  We  aim  to  not  only  improve  our  understanding  of  both  the  long-term 
impact  of  an  oil  spill  on  carbon  cycling  in  a  salt  marsh  environment,  but  also  gain  insight 
into  processes  controlling  the  cycling  of  natural  organic  matter  in  these  complex 
heterogeneous  systems. 

The  salt  marsh  examined  in  this  study  is  situated  in  West  Falmouth,  MA  and 
whilst  the  concentration  of  total  petroleum  hydrocarbons  (TPHs)  has  not  markedly 
changed  over  the  past  30  years  (Reddy  et  al.,  2002),  analysis  by  comprehensive  two- 
dimensional  gas  chromatography  (GCxGC)  indicates  that  petroleum-derived  n-alkanes, 
which  are  typically  abundant  in  No.  2  fuel  oil,  are  no  longer  detectable  (Frysinger  et  al., 
2003).  Petroleum-derived  n-alkanes,  therefore,  are  more  susceptible  to  loss  processes, 
such  as  evaporation,  water  washing  and  biodegradation,  than  the  other  oil-derived 
compounds,  such  as  branched  alkanes,  cycloalkanes  and  polycyclic  aromatic 
hydrocarbons  (PAHs)  that  comprise  the  persistent  UCM.  This  preferential  loss  is  in  direct 
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contrast  to  the  refractory  nature  of  higher  plant  n-alkanes  that  are  known  to  persist  in 
sediments  (Prahl  et  al.,  1997)  and  manifest  themselves  in  West  Falmouth  sediments 
including  the  oil  contaminated  horizons. 

Long-chain  (>«-C23)  alkanes  with  an  odd  over  even  predominance  are  derived 
from  higher  plant  epicuticular  waxes,  and  act  as  protective  coatings  on  leaf  surfaces 
(Eglinton  and  Hamilton,  1967).  They  are  of  particular  interest  because  their  abundance 
and  chemical  composition  varies  under  different  environmental  conditions  as  part  of  a 
physiological  response  by  the  plant.  In  addition,  their  ubiquity  and  refractory  nature 
relative  to  other  organic  compounds  renders  them  well  preserved  in  the  sediment  record 
(Prahl  et  al.,  1997).  These  properties  have  enabled  them  to  be  used  to  reconstruct 
palaeoenvironmental  changes  (Cranwell  1973;  Boom  et  al.,  2002).  Long-chain  (W-C24  to 
H-C34)  fatty  acids  with  even-over-odd  dominance  are  also  derived  from  plant  waxes 
whereas  their  short-chain  (n-Cn  to  n-Cri)  counterparts  are  derived  from  a  broader  range 
of  sources  including  marine  algae  and  bacteria.  Microbial  biomass  may  also  contribute 
significantly  to  fatty  acid  contents  in  marine  sediments  (Hollander  and  Smith,  2001). 

Compound  specific  isotope  analysis  is  a  powerful  tool  used  to  determine  sources 
of  organic  matter  in  marine  sediments  (for  example,  Freeman  et  al.,  1990;  Eglinton  et  al., 
1997;  Pearson  and  Eglinton,  2000).  Both  natural  abundance  radiocarbon  (A14C)  and 
stable  carbon  (813C)  of  individual  compounds  may  be  examined,  but  A14C  is  especially 
useful  as  it  has  a  greater  dynamic  range  (-1000  to  ~^-200%o)  than  813C  (~-32  to  -12%o) 
(Raymond  and  Bauer,  2001).  Due  to  their  geologic  age,  petroleum  residues  contain  no 
14C,  and  this  characteristic  has  previously  been  exploited  to  determine  the  fate  of 
petroleum  residues  in  these  West  Falmouth  salt  marsh  sediments  by  examining  the  A14C 
values  of  different  bulk  sedimentary  organic  matter  pools  (White  et  al.,  2005a).  The  latter 
study  concluded  that  the  majority  of  petroleum  residues  are  solvent-extractable  and  have 
not  been  incorporated  in  the  insoluble  organic  matter  of  the  sediment.  Natural  abundance 
14C  measurements  on  individual  phospholipids  have  also  been  examined  at  this  site  and 
these  analyses  indicate  that  the  oil  is  not  currently  being  degraded  by  the  resident 
microbial  population  (Slater  et  al.,  2005).  These  previous  studies  provide  an  isotopic 
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framework  to  not  only  enhance  our  understanding  of  the  fate  of  petroleum  residues,  but 
also  lipid  biomarkers  from  various  sources  in  marine  sediments. 

6.2  MATERIALS  AND  METHODS 

6.2.1  Study  Area  and  Samples 

A  15  cm  diameter  sediment  core  was  collected  in  August  2000  from  the  intertidal 
zone  at  the  M-l  station  in  Wild  Harbor,  West  Falmouth,  MA  (41°38.0’N,  70°38.5’W)  as 
previously  described  (Reddy  et  al.,  2002).  The  core  was  extruded  and  sectioned  at  2  cm 
intervals  from  0  to  36  cm  and  the  sections  placed  in  pre-combusted  glass  jars  and  frozen 
(-20°C)  until  further  analysis.  The  sediment  was  later  thawed  and  wet-sieved  through  a 
1.4  mm  sieve  to  remove  the  majority  of  intact  roots  and  other  discrete  vascular  plant 
debris  of  the  marsh  grass  Spartina  alterniflora  that  grows  at  this  intertidal  site.  The 
sieved  sediment  was  left  to  air-dry  at  room  temperature  before  homogenizing  with  a 
mortar  and  pestle. 

6.2.2  Lipid  Extraction  and  Purification 

Approximately  20  g  of  air-dried  and  homogenized  sediment  from  0-2,  2-4,  12-14, 
30-32  and  32-34  cm  sediment  horizons  was  extracted  with  dichloromethane  (DCM)  and 
methanol  (MeOH)  (9:1)  by  pressurized-fluid  extraction  (100  °C,  1000  psi).  The  resulting 
total  lipid  extracts  (TLE)  were  reduced  in  volume  by  rotary  evaporation,  dried  over 
Na2SC>4  and  percolated  through  a  copper  column  to  remove  sulfur.  The  TLEs  were  then 
separated  into  four  fractions  using  silica  gel  column  chromatography  (50g,  100-200 
mesh,  fully  activated),  eluting  with  5%  DCM  in  hexane  (Fj),  DCM/hexane  (1:1;  F2),  2% 
formic  acid  in  DCM  (F3)  and  2%  formic  acid  in  MeOH  (F4).  After  reduction  in  volume 
by  rotary  evaporation,  the  hydrocarbon  fractions  (Fi)  were  eluted  through  silver  nitrate 
(~10%  w/w)  impregnated  silica  gel  columns  with  hexane,  to  provide  saturated 
hydrocarbon  fraction  (Fm)  containing  n-alkanes. 
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Further  separation  was  necessary  for  the  12-14  cm  TLE  as  it  also  contained  an 
unresolved  complex  mixture  (UCM)  of  saturated  hydrocarbons  in  the  boiling  range  of  n- 
C13  to  n-C-24,  consistent  with  a  biodegraded  No.  2  fuel  oil  source  (Reddy  et  al.,  2002).  The 
Fm  of  the  UCM  consists  of  20%  branched  alkanes,  50%  1-ring  cycloalkanes  and  30%  2- 
ring  and  greater  cycloalkanes  (Frysinger  et  al.,  2003).  To  separate  these  compounds  from 
the  M-alkanes,  Fm  was  placed  into  a  pre-combusted  glass  vial  containing  a  zeolite 
molecular  sieve  (Geokleen,  GHGeochemical  Services,  3g,  activated  at  450  °C  for  4  hr) 
and  heated  at  80  °C  for  8  hr  (procedure  modified  from  McDuffee  et  al.,  2004).  The 
mixture  was  then  allowed  to  cool  before  passing  through  a  glass  pipette  plugged  with 
glass  wool.  The  sieve  was  then  washed  with  hexane  (20  ml)  to  remove  the  branched  and 
cyclic  hydrocarbons  that  comprised  the  UCM.  The  zeolite  was  allowed  to  dry  overnight 
and  was  then  transferred  to  a  Teflon®  vial  that  had  been  pre-cleaned  with  rinses  of 
MeOH,  DCM  and  hexane.  Concentrated  hydrofluoric  acid  (48%  HF)  was  added  drop 
wise  to  dissolve  the  zeolite  sieve.  The  sample  was  then  extracted  with  hexane  to  yield  a 
clean  M-alkane  fraction,  and  dried  over  Na2SC>4. 

The  F3  acidic  fraction  of  the  TLE  was  reduced  in  volume  and  back  extracted  into 
DCM  before  separating  into  two  fractions  by  aminopropyl-bonded  silica  gel  column 
chromatography  (20  x  0.5cm),  eluting  with  DCM/acetone  (9:1;  F3.1)  and  2%  formic  acid 
in  DCM  (F3.2).  The  acid  fraction  F3.2  containing  “free”  fatty  acids  (FFA)  was 
transesterified  using  acidic  methanol  of  known  carbon  isotopic  composition  (S13C  =  - 
39.59  %<>;  A14C  =  -989.5  %o)  and  purified  by  silica  gel  column  chromatography  (20  x 
0.5cm)  to  yield  a  fatty  acid  methyl  ester  (FAME)  fraction  in  5%  ethyl  acetate  in  hexane. 
This  was  finally  eluted  through  silica  gel  columns  impregnated  with  silver  nitrate  (~10% 
w/w)  using  DCMrhexane  (1 : 1)  to  provide  only  the  saturated  homologues. 

To  obtain  “bound”  fatty  acids  (BFA),  approximately  20  g  of  the  solvent  extracted 
sediment  from  0-2,  2-4,  12-14,  30-32  and  32-34  cm  sediment  horizons  was  taken  and 
refluxed  at  70  °C  for  2  hr  with  a  solution  of  0.5  N  potassium  hydroxide  (KOH)  in  MeOH 
(100  mL)  and  water  (H2O)  (20  mL).  After  cooling,  the  reaction  mixture  was  separated  by 
centrifugation  (1500  rpm)  and  the  remaining  sediment  rinsed  with  MeOH,  followed  by 
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DCM,  and  then  hexane.  The  alkaline  solution  was  back  extracted  with  hexane  (3  x  50 
mL)  to  obtain  a  neutral  fraction  before  acidification  with  4N  hydrochloric  acid  (HC1)  to 
pH  2,  and  extraction  with  DCM  (3  x  50  mL)  to  obtain  an  acidic  fraction.  Both  fractions 
were  reduced  in  volume  by  rotary  evaporation  and  dried  over  Na2S04.  The  acidic  fraction 
was  taken  and  transesterified  and  purified  by  silica  gel  column  chromatography  (20  x 
0.5cm)  and  silver  nitrate  (~10%  w/w)  impregnated  silica  gel  columns  as  previously 
described. 

Purified  n-alkanes,  and  bound  and  free  FAMEs  for  each  sediment  horizon  were 
quantified  by  GC/FID. 

6.2.3  Preparative  Capillary  Gas  Chromatography  (PCGC) 

Select  compounds  were  isolated  from  these  fractions  for  radiocarbon  (14C) 
analysis  by  preparative  capillary  gas  chromatography  (PCGC)  using  an  HP  5890  Series  II 
GC  system  coupled  to  a  Gerstel  preparative  fraction  collector  (see  Eglinton  et  al.,  1996 
for  detailed  description)  fitted  with  a  DB-XLB  (60m  x  0.53mm  i.d.  x  0.5fxm  film 
thickness)  capillary  column.  GC  temperature  programs  for  all  samples  were  40°C  (1 
min),  20°C/min  to  160°C  and  4°C/min  to  320°C  (30  min).  The  U-tube  traps  were  kept  at 
room  temperature  for  w-alkane  samples  and  all  FAMEs  greater  than  H-C20.  FAMEs 
smaller  than  n-C 20,  were  collected  with  the  U-tube  traps  chilled  to  0°C  to  prevent  volatile 
losses.  In  order  to  provide  sufficient  mass  (>30pg  carbon)  for  AMS  14C  analysis  C14,  Ci6 
and  Cj8  free  FAMEs  were  combined  in  the  12-14  cm  and  32-34  cm  sediment  horizons. 
C]4  and  Ci8  free  FAMEs  were  combined  and  Ci6  was  analyzed  separately  in  the  0-2  cm 
interval.  For  bound  FAMEs  C14,  Ci6  and  Cis  homologues  were  combined  in  the  0-2  cm 
and  32-34  cm  horizons  and  C14  and  Cis  homologues  were  combined  and  Ci6  analyzed 
separately  in  the  12-14  cm  interval.  For  w-alkanes,  C27,  C29  and  C31  were  analyzed 
separately  and  the  even  numbered  chains  C26  through  C32  were  combined. 

Column  bleed  was  removed  from  pure  «-alkanes  by  eluting  the  recovered  material 
with  hexane  through  silica  gel  columns  (prepared  in  Pasteur  pipettes)  into  pre-combusted 
quartz  tubes.  For  the  FAME  fractions,  the  column  bleed  was  removed  by  eluting  the 
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recovered  material  with  10%  ethyl  acetate/hexane  through  silica  gel  columns  into  pre¬ 
combusted  quartz  tubes.  The  solvent  was  evaporated  from  all  quartz  tubes  under  ultra- 
high-purity  nitrogen  gas. 

6.2.4  Radiocarbon  Analysis 

Pre-combusted  copper  oxide  was  added  to  the  quartz  tubes  containing  the  samples 
and  the  tubes  were  then  attached  to  a  vacuum  line,  evacuated,  flame  sealed  and 
combusted  to  CO2  (850°C  for  5  hrs).  The  CO2  was  then  purified  on  a  vacuum  line.  The 
stable  carbon  isotopic  ratio  (8I3C)  of  the  samples  were  determined  by  isotope  ratio  mass 
spectrometry  (irMS),  and  radiocarbon  (14C)  content  measured  by  accelerator  mass 
spectrometry  (AMS)  at  the  National  Ocean  Sciences  Accelerator  Mass  Spectrometry 
(NOSAMS)  facility  at  Woods  Hole  Oceanographic  Institution  (WHOI),  Woods  Hole, 
MA.  Radiocarbon  analyses  were  performed  after  conversion  of  the  CO2  to  graphite 
(McNichol  et  al.,  1994).  All  14C  measurements  are  normalized  to  8*3C  values  of  -25%o 
and  expressed  as  AI4C  values,  which  is  the  per  mille  (%o)  deviation  from  the  international 
standard  for  14C  dating,  Standard  Reference  Material  4990B  "Oxalic  Acid".  Precision  for 
813C  and  A14C  measurements  are  ~0.1  and  4-2  l%o  (see  Appendix  for  errors  specific  to 
each  sample  analyzed),  respectively.  For  each  sample,  the  reported  A14C  =  \fme  (1950’9\. 
1]  1000  (Stuiver  and  Polach,  1977),  where  1/8267  (year '),  fm=  fraction  modem  14C, 
corrected  for  isotopic  fractionation  using  813C,  and  “x”  equals  the  year  of  collection.  The 
813C  and  A14C  values  of  “free”  and  “bound”  FA  were  corrected  for  the  addition  of  the 
derivative  carbon  using  a  two  end-member  isotopic  mass  balance. 

6.3  RESULTS  AND  DISCUSSION 

6.3.1  Fatty  Acid  and  Alkane  Distributions 

Free  fatty  acids  from  n-Cu  to  M-C32  (Figure  6.1)  and  bound  fatty  acids  from  n-Cn 
to  n- C26  (Figure  6.2)  are  observed  in  all  horizons  and  exhibit  a  strong  even-over-odd 
dominance.  Short-chain  (n-Cn  to  «-Cig)  homologues  are  the  dominant  FFAs  in  the 


183 


2.0  -| 


O 

o 


14  15  16  17  18  19  20 


■ 

0-2  cm 

□ 

12-14  cm 

□ 

32-34  cm 

22  23  24  25  26  27  28  29  30  31  32 


Chain  length 


Figure  6.1  Abundance  of  individual  free  fatty  acids  from  sediment  horizons  0-2,  12-14 
and  32-34  cm. 
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Figure  6.2  Abundance  of  individual  bound  fatty  acids  from  sediment  horizons  0-2,  12-14 
and  32-34  cm. 
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surface  (0-2  cm)  sediment  whereas  their  long-chain  counterparts  are  more  significant  at 
depth.  Summed  concentrations  of  FFAs  decrease  significantly  from  5.0  mg  g'1  of  organic 
carbon  (OC)  in  the  surface  sediment  to  0.8  mg  g'!  OC  at  12-14  and  32-34  cm.  In  contrast, 
BFA  contents  do  not  change  much  between  depths  of  0-4  and  12-14  cm  with  summed 
concentrations  of  2.5  mg  g"1  OC.  At  30-34  cm,  however,  the  content  decreases  to  0.5  mg 
g'1  OC.  It  is  not  clear  whether  this  change  in  abundance  of  BFAs  that  occurs  in  the 
deepest  sediment  horizon  results  from  diagenetic  processes  or  a  change  in  the  source  of 
organic  matter. 

The  abundance  of  n-alkanes  in  0-2,  12-14  and  32-34  cm  sediment  horizons  are 
shown  in  Figure  6.2.  The  most  abundant  homologue  in  all  horizons  is  n- C29  and 
concentrations  of  individual  n-alkanes  decrease  slightly  with  increasing  depth  within  the 
sediment  from  0.22  mg  g'1  OC  (0-2  cm),  to  0.19  mg  g'1  OC  (12-14  cm)  and  0.16  mg  g'1 
OC  (32-34  cm).  There  is  a  strong  odd-over-even  carbon  number  predominance  in  all 
sediment  horizons  which  can  be  numerically  represented  by  carbon  preference  index 
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Figure  6.3  Abundance  of  individual  n-alkanes  (C25  through  C33)  from  sediment  horizons 
0-2, 12-14  and  32-34  cm. 
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(CPI,  =  £[odd  numbered  alkane  abundances  from  n-C 23  to  «-C35]/Z[even  nmnbered 
alkane  abundances  from  H-C22  to  n-C34]).  The  CPI  increases  slightly  with  depth,  from  4.4 
at  0-2  cm  to  7.6  in  12-14  and  8.5  at  32-34  cm.  These  values  are  within  ranges  previously 
exhibited  for  plant  epicuticular  waxes  (CPI  =  4-10;  Tulloch,  1976). 

Rates  of  decrease  in  these  biomarkers  over  the  depth  examined  (0-34  cm)  are 
~0.12  mg  g'1  OC  cm'1  for  free  fatty  acids  and  ~0.02  mg  g'1  OC  cm"1  for  n-alkanes, 
highlighting  the  preferential  degradation  of  free  fatty  acids. 

6.3.2  Compound  Specific  Isotope  Analysis  of  Fatty  Acids 

Compound  specific  A14C  and  813C  values  for  free  and  bound  fatty  acids  (FFA  and 
BFA  respectively)  and  bulk  A14C  and  813C  values  for  the  solvent-extracted  sediment 
residue  (EX-RES)  are  shown  in  Figure  6.3.  The  isotopic  composition  of  EX-RES 
represents  the  bulk  OC  uncontaminated  by  petroleum  (White  et  al.,  2005a).  The  813C 
values  in  this  study  agree  well  with  those  for  organic  carbon  in  salt  marsh  sediments 
where  Spartina  is  the  dominant  macrophyte.  Values  of  813C  between  -24  and  -17%o 
reflect  a  mixture  of  phytoplankton,  benthic  algae  and  vascular  plant  material,  including 
the  13C  enriched  C4  macrophyte  Spartina  (Boshker  et  al.,  1999).  In  this  study,  the  813C 
for  short-chain  (h-Ci4_]8)  FFA  (filled  triangle  and  diamond  symbols)  are  relatively  similar 
throughout  the  core  (-23.6  to  -22.7%o),  and  are  consistent  with  a  predominantly  algal 
and/or  bacterial  source  and  could  also  be  Spartina  derived.  The  H-C14-18  BFA  are 
somewhat  more  13C-depleted  (-25.6  to  -22.9%o)  and  likely  have  a  more  significant 
contribution  from  marine  sources.  The  long-chain  (H-C24-28)  FFA,  however,  are  much 
more  13C-enriched  (-20.9  to  -18.7%o)  consistent  with  an  origin  from  Spartina. 

A14C  values  of  short  and  long-chain  FFA  and  short-chain  BFA  from  the  surface 
sediment  (0-2  cm)  are  very  close  to  that  of  the  solvent  extracted  sediment  residue  (EX- 
RES),  a  finding  similar  to  that  for  A14C  of  phospholipids  fatty  acids  (PLFAs)  previously 
examined  (Slater  et  al.,  2005).  The  S13C  and  A14C  values  of  the  EX-RES  are  driven  by  the 
dominant  input  of  Spartina  (White  et  al.,  2005a),  but  also  reflect  a  contribution  from 
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marine  inputs.  With  depth,  PLFAs  examined  by  Slater  et  al.,  (2005),  and  the  long-chain 
FFA,  and  short-chain  BFA  (examined  in  this  study)  continue  to  track  the  A14C  EX-RES, 
with  slight  enrichments  in  14C  for  the  BFAs  and  PLFAs.  The  short-chain  FFA,  however, 
become  significantly  more  depleted  in  14C  with  depth,  accompanied  by  a  significant 
decrease  in  their  concentration,  as  previously  noted.  The  decrease  in  concentration  may 
result  from  biodegradation  of  the  fresher  and  labile  material,  which  is  likely  the  most 
enriched  in  14C,  and  hence  the  A14C  of  the  short-chain  FFA  becomes  more  depleted  with 
depth.  The  short-chain  FFA  that  persist  are  significantly  older  and  may  have  a  different 
source  such  as  resuspended  marine  carbon  that  is  washed  in  by  the  tide,  and  would  likely 
be  older  (Megens  et  al.,  2001),  and  potentially  more  resistant  to  degradation  due  to 
association  with  minerals  (Middleburg  et  al.,  1997). 

The  peak  in  above-ground  testing  of  thermonuclear  weapons  reached  a  maximum 
at  1963-1964.  This  corresponds  to  a  depth  of  approximately  13  cm  in  the  salt  marsh 
sediment  assuming  a  sedimentation  rate  of  0.35  cm  yr'1,  and  is  consistent  with  a 
maximum  of  137Cs  at  this  depth  (White  et  al.,  20056).  Whilst  long-chain  FFA  and  short- 
chain  BFA  are  relatively  enriched  at  this  depth  (142%o  and  187-210%o),  a  pronounced 
bomb-spike  is  not  observed  in  the  A14C  of  the  long-chain  FA  as  might  be  expected  for 
vascular  plant  tissue  such  as  that  predominantly  derived  from  Spartina.  This  would  be 
expected  as  a  result  of  photosynthetic  fixation  of  atmospheric  carbon  dioxide  (CO2)  that 
would  have  been  greatly  enriched  in  14C  at  this  time.  The  amount  of  14C  in  the 
atmosphere  nearly  doubled  in  the  early  1960s,  and  the  A14C  of  tropospheric  CO2 
increased  to  a  maximum  value  of  >900%o  (Levin  et  al.,  1985). 

The  absence  of  a  pronounced  bomb-spike  in  the  A14C  of  the  long-chain  FA  may 
arise  from  a  smearing  of  bomb-produced  14C  as  a  result  of  sedimentary  mixing,  which 
would  in  turn,  dampen  the  amplitude  of  the  bomb-spike.  This  would  also  be  evident  in 
the  activity  profile  of  137Cs  at  the  site,  which  is  present  throughout  depths  of  9-15  cm 
(White  et  al.,  20056)  corresponding  to  the  years  1951-1974.  The  first  detection  of  137Cs  in 
these  sediments  shows  good  correlation  with  the  early  1950s  increase  in  total  fission 
yields  (Carter  and  Moghissi,  1977)  and  the  general  time  span  agrees  well  with  the 
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Figure  6.4  A14C  and  813C  compositions  of  individual  free  (filled  symbols)  and  bound 
(open  symbols)  short  (triangles  and  diamonds)  and  long-chain  (squares)  fatty  acids  from 
0-2, 12-14  and  32-34  cm  compared  to  the  A14C  and  813C  composition  of  the  solvent 
extracted  sediment  residue  (EX-RES)  from  White  et  al.,  2005a.  Errors  for  A14C  and  813C 
are  smaller  than  the  symbols  representing  the  data  and  thus  are  not  observed. 
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atmospheric  history  of  137Cs  fallout  from  nuclear  weapons  testing  previously  documented 
from  1952-1970  in  Miami,  Florida  USA  (Robbins  et  al.,  2000).  It  therefore  does  not 
appear  that  the  suppression  of  14C  in  these  sediments  is  due  to  mixing.  This  is  further 
supported  by  comparison  to  the  137Cs  activity  profile  in  anoxic  and  laminated  sediments 
from  the  Pettaquamscutt  River  where  sediment  mixing  is  absent  (Lima  et  al.,  2005).  The 
peak  in  137Cs  activity  corresponding  to  nuclear  weapons  testing  at  this  site,  spans  a 
similar  time  frame  (1953-1974)  and  peak  width  (8  cm;  14-22  cm)  to  that  observed  in 
West  Falmouth  sediments.  It  is  worth  noting  that  the  137Cs  peak  in  West  Falmouth 
sediments  is  not  as  smoothed  as  that  observed  in  the  laminated  sediments  of  the 
Pettaquamscutt  River  and  analysis  at  higher  resolution  would  more  definitively  determine 
the  significance  of  sedimentary  mixing  at  the  site. 

It  is  unlikely  that  there  is  a  lag  between  the  incorporation  of  bomb  14C  into  the 
Spartina  and  its  deposition  in  the  sediments  as  it  is  a  perennial  plant  and  incorporates 
CO2  year-round.  A  further  alternative,  therefore,  is  that  a  pronounced  bomb-spike  is 
simply  not  recorded  in  these  sediments.  This  absence  could  also  occur  if  the  Spartina 
fixes  inorganic  carbon  from  sources  other  than  the  atmosphere.  Carbon  dioxide  derived 
from  sediments  via  the  rhizome  system  and  transported  by  mass  flow  can  be  responsible 
for  more  than  90%  of  the  total  carbon  uptake  in  emergent  plants  (Madsen  and  San- 
Jensen,  1991).  Such  a  mechanism  would  easily  account  for  the  observations  in  this  study 
if  marine  porewater  dissolved  inorganic  carbon  (DIC)  was  the  source  of  inorganic 
carbon.  However,  in  large  emergent  plants  such  as  Spartina  altemiflora,  porewater  DIC 
is  not  thought  to  be  a  major  contribution  because  the  distance  from  root  to  shoot  is  too 
long  to  support  high  rates  of  photosynthesis  in  the  leaves  (Bowes,  1993).  For  Spartina 
altemiflora,  only  1-10%  of  carbon  for  photosynthesis  has  been  observed  to  be  sediment 
derived  (Hwang  and  Morris,  1992)  and  whilst  this  may  lead  to  a  slightly  muted  bomb 
signal  in  the  long-chain  FA  it  would  not  dilute  it  to  the  value  observed.  The  study  by 
Hwang  and  Morris,  however  was  performed  on  Spartina  grown  in  the  laboratory  and 
therefore  without  direct  14C  measurements  of  Spartina  biomass,  and  analysis  of  uptake  of 
DIC  in  the  field,  this  cause  cannot  be  discounted. 
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6.3.3  Compound  Specific  Isotope  Analysis  of  Alkanes 

Compound  specific  A14C  and  813C  for  H-alkanes  are  shown  in  Figure  6.4.  As  for 
the  fatty  acids,  the  813C  of  H-alkanes  becomes  more  enriched  with  depth,  from  -26.1  to 
-24.6%o  (0-2  cm)  to  -24.0  to  -22.8%o  (12-14  cm)  to  -22.4  to  -20.8%o  (32-34  cm)  indicating 
selective  preservation  of  Spartina-derived  carbon.  In  general  813C  values  of  H-alkanes  are 
~4-6%o  more  depleted  than  fatty  acids,  and  this  is  attributed  to  additional  fractionation 
against  13C  during  their  biosynthesis  comparing  well  with  previous  work  (Collister  et  al., 
1994). 

Consistent  with  observations  made  for  fatty  acids,  significant  enrichment  resulting 
from  the  bomb-spike  is  not  observed  in  the  plant  wax  derived  H-alkanes  (C29  and  C31)  at 
12-14  cm  although  they  are  more  enriched  in  14C  than  the  EX-RES  (206-26  l%o  vs. 
164%o)  in  this  interval.  Instead,  these  H-alkanes  have  similar  A14C  values  to  short-chain 
BFA  and  long-chain  FFA  previously  discussed.  It  is  evident,  therefore  that  all  of  these 
compounds  are  predominantly  derived  from  Spartina.  No  evidence  is  found  for  a 
contribution  from  petroleum  hydrocarbons  despite  the  depth  corresponding  to  the  peak 
contamination.  Indeed,  the  odd-chain  H-alkanes  (C29  and  C31)  are  more  enriched  in  14C 
than  the  long-chain  FAs  in  the  12-14  cm  interval,  likely  due  to  the  enhanced  preservation 
of  these  compounds,  and  a  greater  legacy  of  compounds  produced  during  a  bomb  spike. 

In  the  surface  sediment  and  12-14  cm  interval,  the  even-chain  H-alkanes  are 
significantly  more  depleted  in  14C  than  their  corresponding  odd-chain  counterparts.  For 
example  at  0-2  cm,  the  even  homologues  exhibit  a  A14C  value  of  -476%o  whereas  n- C27  is 
more  enriched  (-224%o),  and  H-C29  is  further  enriched  (0%o).  This  clearly  indicates  the 
input  of  fossil  carbon  that  has  a  larger  impact  on  the  even-chain  H-alkanes  than  it  does  on 
the  odd  homologues.  This  is  due  to  the  lack  of  carbon  number  preference  in  thermally 
mature  fossil  sources  (Brassell  and  Eglinton,  1980)  relative  to  the  odd-over-even 
predominance  from  vascular  plant  sources.  The  fossil-fuel  carbon  inputs  therefore  have  a 
disproportionate  impact  on  the  A14C  value  of  the  even  homologues  because  they  are  not 
abundantly  produced  by  terrestrial  or  marine  biomass.  The  fossil  carbon  input  impacts 
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Figure  6.5  A14C  and  813C  compositions  of  individual  n-alkanes  from  0-2, 12-14  and  32- 
34  cm  compared  to  the  A14C  and  813C  composition  of  the  solvent  extracted  sediment 
residue  (EX-RES)  from  White  et  al.,  2005a.  Errors  for  A14C  and  813C  are  smaller  than  the 
symbols  representing  the  data  and  thus  are  not  observed. 
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n- C29  less  than  /1-G27  as  a  result  of  the  increased  abundance  of  M-C29  in  plant  waxes. 
Similar  observations  are  also  made  at  32-34  cm,  although  the  offset  is  more  subtle.  Mixed 
sources  of  n-alkanes  to  marine  sediments  have  previously  been  observed  in  sediments 
from  the  Santa  Monica  Basin  (Pearson  and  Egliriton  2000).  In  this  study,  A14C  and  813C 
measurements  of  long-chain  alkanes  (H-C24-33)  were  measured  and  a  three-  component 
mixing  model  used  to  assess  the  contribution  of  petroleum  (12%),  modem  plant  wax 
(80%)  and  shale-derived  alkanes  (8%). 

The  fraction  of  n-alkanes  that  are  fossil-derived  in  these  sediment  horizons  can  be 
calculated  based  on  isotope  mass  balance.  Assuming  that  the  total  n-alkanes  in  the 
sediment  are  a  binary  mixture  (?)  comprised  of  fossil  carbon  (x)  and  Spartina  plant  waxes 
(n)  that  can  be  represented  by 

(alkane)/  =  (alkane)*  +  (alkane)/?  (I) 

and  the  isotopic  composition  of  the  mixture  is  given  by 

Ft  =  fFx  +  (l  -  f)Fn  (II) 

where  the  Fs  represent  isotopic  abundances  and  /  is  the  fraction  of  the  mixture  that  is 
fossil-derived.  In  terms  of  radiocarbon  abundances  (A),  we  can  write 

f  =  {At -An) /{Ax -An)  (III) 

The  corrected  A14C  values  for  the  total  n-alkanes  were  used  as  At  and  the  corrected  A14C 
values  for  the  long-chain  fatty  acids  were  used  to  represent  An.  Assuming  Ax  to  reflect  an 
absence  of  14C  (i.e.,  A,4C  =  -1000%o),  solutions  of f  were  obtained  for  the  n-alkanes  in  the 
three  sediment  horizons  examined  (Table  6.1).  The  isotope  mass  balance  calculation  for 
sediments  from  0-2  and  12-14  cm  reiterates  that  a  greater  fraction  of  even-chain  n- 
alkanes  are  fossil  derived  compared  to  the  odd-chain  lengths  (0.24  to  0.51  vs.  -0.10  to 
0.27).  This  offset  is  not  observed  in  sediments  from  32-34  cm  and  in  general  the  fraction 
of  fossil  carbon  is  much  less  for  all  n-alkanes  (0  to  0.10)  in  this  horizon.  Isotope  mass 
balance  also  indicates  that  the  fraction  of  fossil  carbon  is  less  for  n-C29  less  than  n-C27  in 
sediments  from  0-2cm  and  32-34  cm  as  previously  described.  Negative  contributions  of 
fossil  carbon  in  sediments  from  12-14  cm  arise  from  the  fact  that  A14C  values  for  the 
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long-chain  fatty  acids  used  to  represent  the  Spartina  plant  waxes  are  more  enriched  in  14C 
than  the  w-alkanes  due  to  the  influence  of  the  bomb-spike  at  this  depth. 

The  source  of  this  fossil  carbon  to  the  surface  and  12-14  cm  sediment  horizon  is 
likely  fossil  «-alkanes  entrained  on  combustion-derived  particles  derived  from  local 
and/or  distal  sources.  The  chain  length  of  the  n-alkanes  up  to  C32  indicates  a  source  from 
motor  boat  fuel  such  as  No.  6  fuel  oil.  These  H-alkanes  may  be  physically  protected  from 
degradation  by  strong  association  with  combustion  particles  (McGroddy  et  al.,  1995),  or 
it  may  be  that  their  long-chain  lengths  impart  a  chemical  recalcitrance.  Plant-wax 
«-alkanes  may  also  be  more  likely  to  persist  in  sediments  because  of  the  way  in  which 
they  are  packaged.  Leaf  waxes  occur  as  microcrystallites  in  the  pm  range  varying  greatly 
in  size,  shape  and  cuticular  distribution  (Martin  and  Juniper,  1970)  and  in  this  form  are 
less  likely  to  be  bioavailable. 

Table  6.1  Fraction  of  fossil-derived  w-alkanes  in  West  Falmouth  sediments 


depth 

(cm) 

approx. 

age 

#i-alkane 

fraction 

fossil-derived 

0-2 

1997 

n-  C27 

0.27 

AJ-C29 

0.07 

H-C26+28+30+32 

0.51 

12-14 

1963 

«-C29 

-0.06 

«-c3, 

-0.10 

H-C26+28+30+32 

0.24 

32-34 

1906 

n- C27 

0.10 

n-  C29 

0.08 

W-C31 

0.00 

H-C26+28+30+32 

0.07 
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6.4  CONCLUSIONS 


Natural  abundance  14C  measurements  of  specific  biomarkers  provide  a  powerful 
means  to  determine  the  source  of  organic  compounds  in  contemporary  environmental 
matrices.  In  this  study,  contributions  from  petrogenic,  pyrogenic  and  biogenic  sources 
were  assessed.  Salt  marsh  derived  n-hydrocarbons  and  fatty  acids  are  predominantly 
derived  from  the  dominant  macrophyte  Spartina  alterniflora,  although  inputs  from 
marine  as  well  as  pyrogenic  sources  are  important.  Compound-specific  isotope  analysis 
of  long-chain  fatty  acids  and  alkanes  not  only  provides  information  about  sources  of 
organic  matter,  but  also  about  post-depositional  processes  that  occur  in  sediments.  For 
example,  petroleum  derived  n-alkanes  from  the  oil  spill  do  not  persist  in  sediments, 
whereas  pyrogenic  and  biogenic  n-alkanes  do,  likely  as  a  result  of  protective  matrix 
association.  It  is  evident  from  this  study  that  the  mode  of  association  or  “packaging”  of 
biomarkers  plays  a  crucial  role  in  determining  their  long-term  fate.  Here,  we  observe  that 
the  same  compound  can  undergo  very  different  fates  depending  on  how  it  is  introduced 
to,  and  physically  disposed,  in  the  sediment. 
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CHAPTER  7 


CONCLUSIONS  AND  FUTURE  WORK 


Natural  abundance  radiocarbon  (A14C)  measurements  provide  a  powerful  means 
to  determine  the  source  of  organic  compounds  in  contemporary  environmental  matrices. 
Analysis  of  the  14C  content  of  different  fractions  of  sedimentary  organic  matter  has 
enhanced  our  understanding  of  the  associations  between  specific  organic  compounds  and 
macromolecular  organic  matter  (MOM),  and  we  can  conclude  that  these  associations 
influence  and  in  some  instances,  determine  the  long-term  fate  of  organic  compounds. 

Surface  sediments  examined  from  a  variety  of  different  depositional  settings 
display  A14C  values  from  -737%o  to  +44%o,  encompassing  most  of  the  dynamic  range  of 
A14C.  This  highlights  the  significant  inputs  of  both  recently  photosynthesized  carbon,  in 
the  case  of  the  West  Falmouth  saltmarsh,  and  pre-aged  carbon  to  marine  sediments  from 
a  variety  of  depositional  settings  including  the  Eel  River  Margin,  Guaymas  Basin,  Ross 
Sea,  Mackenzie  Shelf  and  Beaufort  Slope. 

In  terms  of  associations  between  natural  organic  carbon  and  MOM,  we  conclude 
that  for  the  majority  of  sediments  studied,  marine  algal-derived  carbon  is  preferentially 
preserved  by  association  with  MOM  via  ester  linkages.  Differences  in  A14C  values  and 
molecular  compositions  between  sediments  and  their  respective  residues  also  highlight 
the  effects  of  diagenesis  (e.g.  Benguela  Upwelling  and  Washington  Margin  sediments) 
and  the  presence  or  absence  of  oxygenated  bottom  waters  (e.g.  Santa  Monica  Basin 
sediments). 

Hydrophobic  organic  contaminants  (HOCs)  persisting  in  marine  sediments,  do  not 
have  to  be  “bound”  or  associated  with  nonextractable  sediment  residues,  to  persist.  This 
is  evident  in  sediments  from  Wild  Harbor,  West  Falmouth,  MA  where  the  majority  of 
No.  2  fuel  oil  residues  are  found  to  be  solvent-extractable  and  not  bound  to  insoluble 
MOM.  Good  agreement  is  observed  between  isotope  mass  balance  calculations  and 
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values  calculated  when  total  oil  concentrations  (measured  by  gas  chromatography)  are 
compared  to  total  organic  carbon  (TOC)  content  (determined  by  elemental  analysis).  This 
highlights  the  utility  of  14C  content  as  an  inverse  tracer  of  fossil  fuel-derived  HOCs  in 
marine  sediments. 

Studies  of  various  other  HOCs  including  creosote,  coal  tar,  polychlorinated 
biphenyls  (PCBs),  benzotriazoles  (BZTs),  2,2-bis(chlorophenyl)-l,l,l-trichloroethane 
(DDT)  provide  similar  conclusions  to  those  made  for  West  Falmouth  sediments  and  the 
majority  of  fossil  fuel-derived  carbon  is  solvent-extractable  and  not  bound  to  MOM.  In 
instances  where  HOCs  and/or  their  respective  metabolites  are  found  to  be  associated  with 
nonextractable  sediment  residues,  their  structures  suggest  that  they  are  not  covalently 
linked  and  are  instead  likely  entrapped  in  mineral  or  organic  matter  matrices.  Overall, 
these  findings  indicate  that  HOCs  may  be  more  available  to  biota,  and  more  likely  to  be 
subjected  to  transport  processes  and  thus  potentially  have  a  greater  impact  on  the 
environments  in  which  they  persist.  Where  evidence  of  covalent  bonding  between  HOCs 
and  MOM  is  observed,  it  is  of  interest  to  determine  whether  binding  occurred  prior  to 
entering  the  aquatic  environment.  It  is  possible  that  binding  may  have  occurred  in  soils  or 
as  a  by-product  of  their  synthesis  and  determination  of  this  will  be  the  focus  of  future 
work. 

Quantities  of  fossil  fuel-derived  carbon  released  by  solvent  extraction  and/or 
saponification  as  determined  by  gas  chromatography  (GC)  do  not  always  agree  with 
quantities  derived  from  isotope  mass  balance  calculations.  Discrepancies  between  GC 
and  isotope  mass  balance  point  to  aspects  that  we  do  not  currently  understand.  Future 
work  will  be  aimed  at  understanding  these  offsets  so  that  the  fate  of  HOCs  in  the  marine 
environment  can  be  more  conclusively  determined. 

Sediment  from  the  hydrothermal  environment  of  the  Guaymas  basin,  where 
petroleum  is  actively  produced  in  the  basin  sediments,  presents  an  alternate  view  of  the 
fate  of  fossil  carbon  in  sediments.  Here,  the  solvent-extractable  fossil  carbon  pool 
consists  primarily  of  n-alkanoic  acids  (supported  by  Pearson  et  al.,  2005),  as  opposed  to 
those  derived  directly  from  petroleum.  Petroleum-derived  compounds  such  as  polycyclic 
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aromatic  hydrocarbons  (PAHs)  are  instead  associated  with  non-extractable  sediment 
residues.  We  therefore  conclude  that  PAHs  are  protected  from  degradation  or  other  loss 
processes  by  their  association  with  MOM  at  this  site.  The  PAHs  are,  however,  non- 
functionalized  and  thus  do  not  appear  to  be  covalently  linked  to  MOM  and  are  instead 
likely  associated  with  minerals  or  encapsulated  in  organic  matter  domains  as  observed  for 
the  HOCs  previously  described. 

Finally,  the  importance  of  the  mode  of  association  between  organic  compounds 
and  MOM  is  emphasized  by  molecular-level  14C  analysis  of  individual  biomarkers  from 
West  Falmouth  sediments.  Here,  compound  specific  isotope  analysis  of  short  chain  («- 
C  14-Cis)  and  long  chain  (M-C24-C28)  fatty  acids  and  long  chain  alkanes  (n- C26-C32)  not 
only  provides  information  regarding  the  sources  of  organic  matter,  but  also  about 
sedimentary  processes  and  interactions  with  the  dominant  macrophyte  at  the  site, 
Spartina  altemiflora.  Petroleum-derived  n-alkanes  from  the  oil  spill  do  not  persist  in 
these  sediments,  but  pyrogenic  and  biogenic  n-alkanes  do.  The  persistence  of  plant-wax- 
derived  n-alkanes  likely  results  from  a  form  of  protective  matrix  association  thus 
providing  evidence  that  mode  of  association  can  play  a  crucial  role  in  determining  the 
long-term  fate  of  biomarkers. 

The  absence  of  a  pronounced  bomb  spike  in  Spartina-derived  biomarkers 
indicates  that  either  sediment  mixing  dampens  the  amplitude  of  the  peak  or  that  there  are 
processes  concerning  the  uptake  of  CO2  by  Spartina  that  are  not  fully  understood  and  are 
more  complex  than  the  fixation  of  atmospheric  CO2  alone.  Future  work  will  be  aimed  at 
performing  Cs  analysis  at  a  higher  resolution  to  constrain  sediment  mixing  in  the  salt 
marsh.  Radiocarbon  analysis  of  individual  Spartina  biomarkers  at  a  higher  resolution  in 
the  sediment  as  well  as  14C  analysis  of  the  Spartina  itself  (including  roots  and  leaves)  and 
sedimentary  dissolved  inorganic  carbon  should  also  be  performed.  These  results  will 
allow  us  to  more  fully  constrain  the  inputs  of  carbon  to  the  saltmarsh  so  that  we  can 
improve  our  understanding  of  organic  carbon  cycling  in  this  environment.  Coastal 
environments  and  shelf  and  slope  sediments  play  an  important  role  in  the  global  organic 
carbon  cycle  at  the  interface  between  land  and  ocean  (Hedges  et  al.,  2000).  Developing  a 
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more  complete  understanding  of  the  processes  that  control  these  systems  and  the  sources 
and  sinks  of  organic  carbon  has  implications  for  understanding  the  global  organic  carbon 
cycle  as  a  whole. 
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